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to  arrive  at  any  operational  conclusion  or  new  discoveries. 
This  Space  Shuttle  Simulation,  which  is  aimed  at  training 
and  contractual  acceptance  of  procured  hardware,  is  different 
from  the  AFFTC  Office  of  Advanced  Manned  Vehicles'  Space 
Shuttle  Simulation,  which  is  used  to  support  Shuttle  develop- 
ment and  flight  testing. 
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PREFACE 


The  following  document  was  intended  to  serve  as  a reference  manual 
which  will  enable  the  reader  to  gain  a quick  understanding  into  the  design 
and  working  arrangements  of  the  Space  Shuttle  Entry  and  Approach  Test 
simulation  as  used  on  the  Simulator  for  Aircraft  Flight  Test  and  Development 
Simulation  Benchmark.  By  reading  through  this  report  he  is  expected  to 
know  how  to  make  use  of  the  program  and  to  find  his  way  around  in  case  any 
modification /addition  is  comtemplated.  In  brief,  the  main  purpose  for  this 
report  is  to  insure  a smooth  transferability  from  the  originators  to  the 
users  with  a minimum  period  of  learning  and  associated  delay. 

The  program  was  written  in  Fortran  for  the  SEL  32/55  digital  computer 
with  48K  600  nsec  memory  equipped  with  floating-point  hardware  and  tied  to 
a functioning  mock-up  of  an  aircraft  cockpit.  Some  part  of  this  program 
was  coded  in  Assembly  language  to  improve  the  execution  speed.  Therefore 
the  package  cannot  be  transferred  to  another  type  of  machine  without  some 
re-programming  efforts. 

Significant  contributions  to  all  phases  of  the  simulation  and  report-writing 
were  made  by  the  following  personnel: 

Mr.  Richard  Hansen,  Supervisor 
Mr.  Steve  E.  Louton,  Engineer 
Mr.  Larry  V.  LeDuc,  Engineer 
Capt.  Charles  L.  Bozeman,  USAF 
Lt.  Douglas  V.  Palmer,  USAF 
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1.  INTRODUCTION 


The  Shuttle  software  package  is  a general-purpose  real-time 
digital  computer  program  that  simulates  all  aspects  of  the  flight: 
aerodynamics,  flight  control  system  connected  to  a cockpit  for 
man-in- the-loop  operation.  The  Space  Shuttle  mission  profile  is 
composed  of  many  modes.  They  are  very  briefly:  ignition  6 lift- 
off, pitch  initialization  after  about  6 seconds,  staging  122 
seconds  after  lift-off  (in  the  neighborhood  of  2 35-nautical-mile 
altitude),  external  tank  separation  6 disposal  then  orbit  inser- 
tion. On  the  way  back  to  earth  the  mission  can  be  subdivided 
into:  deorbit  phase  at  160-nautical-mile  altitude,  entry  mode 
starting  at  400,000ft,  terminal  phase  starting  at  50,000ft  and 
finally  landing  and  touchdown.  The  simulation  contained  in  this 
report  is  devoted  to  the  last  two  phases:  Entry  (beginning  at 
10*  ft)  and  Terminal  Area  Energy  Management  (TAEM) 

Handling  quality  constitutes  the  main  objective  of  this 
study.  It  does  not  address  the  Shuttle's  full  mission  nor  does 
it  simulate  the  navigation  to  pinpoint  a landing  site.  Actual 
touchdown  will  not  be  simulated  because  the  model  does  not 
include  ground  effects,  the  lowering  of  landing  gear  assembly 
and  nose  wheel  steering.  The  Space  Shuttle  simulation  described 
in  this  report  is  designed  for  training  and  contractual  accept- 
ance of  procured  simulator  hardware.  It  is  different  from  the 
Shuttle  simulation  used  by  the  AFFTC  Office  of  Advanced  Manned 
Vehicles  which  is  used  for  development  and  flight  testing  of 
the  Shuttle. 


2 . MATHEMATICAL  MODEL 

By  writing  the  aerodynamic  equations,  three  explicit  coordinate 
systems  are  used.  They  are  H-frame  (flight-path  axes  for  describing 
orbital  mechanics),  body-axes  and  earth  axes.  The  choice  of 
various  coordinate  systems  is  made  to  enhance  the  relative 
computational  accuracy  of  each  segment  of  the  model.  For  example, 
it  makes  no  sense  writing  the  vehicle  rotational  equations  in  any 
coordinate  system  except  the  ones  with  an  origin  residing  within 
the  vehicle  itself  such  as  body-axes,  stability  axes,  or  wind  axes. 

The  transformation  from  body-axes  to  earth-axes  is  accomplished 
with  a matrix  composed  of  direction-cosine  terms.  These  direction 
cosines  are  listed  below  but  they  are  not  explicitly  evaluated  in 
the  program. 

= cos  0 cos  ip 


Z2  = cos  0 sin  ip 


£g  = -sin  0 


( 2 -1) 


™1  = -cos  <J>  sin  \p  + sin  <J>  sin  0 cos  ip 


cos  ip  + sin  <J>  sin  0 cos  iJj 


(2-2) 


<(>  cos  0 


*1  = sin  <p  sin  ip  + cos  <p  sin  0 cos 


2 = -sin  ip  cos  ip  + cos  <p  sin  0 sin  ip 


3 = cos  0 cos  4> 


( 2 -3) 


The  angles  0,  p and  ip  represent  the  pitch,  roll  and  yaw  angles. 
Together  they  are  commonly  referred  to  as  the  Euler  angles. 
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Tor  clarity,  the  model  equations  are  arranged  in  groups  , each 
group  describing  an  aspect  of  the  flight.  Subscripts  II,  D,  and  E 
stand  for  the  three  coordinate  systems  mentioned  earlier.  What- 
ever symbol  not  explicitly  named  can  be  found  in  alphabetical 
order  under  Chapter  8 . 

2.1  Aerodynamic  equations;  Acceleration  terms  in  body-axis 

reference  frame  are  the  results  of  thrust  vectors  T f t , and  . 

x y ^ 

The  subscript  "G"  stands  for  landing  gear  and  even  though  there 
is  no  engine,  the  thrust  terms  are  included  in  the  equations  for 
completeness . 


XR  = =& 
m 


'c0  + 


!l_  + SdF  + C°6SF  <{SIi> 


C ( G)  ! T 
+ G I + — 

J m 


Yu  = { [Cy  ( 6SB  - 6SBnom)  + Cy  1 6r  + Cy  (6a)  (2-4) 


+ Cy  (8)  } + — — 
\3  m 


ZB  = tCNo  + CN6c  (6o)  + CNxnv  (6BF)  + CNxsn  ( 6SD  - 


6S\-0J!) 


* C«5G  («0)  ♦ AC  (CN  >1  + IT 

(h) 


Note  that  the  subscript  "N"  refers  to  the  normal  force  as  opposed 
to  subscript  "n"  which  is  associated  with  yawing-moment  coefficients. 


Resolved  into  the  earth-axes  reference  frame,  these  acceleration 
terns  become,  after  compensation  for  the  earth’s  oblateness 


r ’ 

E = XBS,1  + VL,m1  + ZDn1  - 3J2g  — 


cos  p sin  p 


E ” XD  l‘A  + YEm2  + ZDn2 


( 2 -5) 


E = xb£n  + Yurnn  + zBnn 


Y.lirrc  u and  J0  am  tlie  colatitude  (^-  - I.)  and  the  neasure  of  cart!, 
oblateness  respectively.  The  dynamic  pressure  (j  in  the  above 
equations  and  other  related  terms  such  as  Mach  number  and  the  mass 
n are  defined  below. 

— 1 n 

q 4py^ 


Mach 


_ V 


— ; a = sonic  speed 
a ' 


( 2 -6) 


m = 


... 

•*o  + r 


Here,  'n  and  l.j.  stand  for  the  vehicle  empty  weight  and  fuel  weight 
respectively.  Also  the  airspeed  V,  which  is  normally  obtained  by 


WUj,"  + Vj.^  + can  simply  be  established  by 

V = 


UD 


cos  a • cos  h 


(2-7) 


Latitude  I.  and  longitude  X are  derived  from  the  launching-point 
coordinate  (Ln,  X0)  and  tlie  earth’s  rotation  ue 


X = X. 


2 . 2 


V 

+ 

1*1 

r cos  L e 

Jo 

J 

Equations:  hirst 

dt 


( 2 -8 ) 


phase.  Lot  \J>j,  be  the  angle  between  the  direction  of  travel  and 
true  nortli  and  r the  vehicle  altitude  measured  from  the  earth  center, 
Let  (Xjj,  yjj)  be  the  local  plane  perpendicular  to  the  line  of  r 
(or  zT.)  such  that  x.,  is  pointing  in  the  general  direction  of  travel 

.1  i» 

and  yI(  perpendicular  to  (x^,,  z,,)  to  form  a r ight -handed  triad  (see 

t 


ref.  9.1.  3)  Then 


rIiI  = CrlTlf>0  + 


r::.,dt 

i i 


+ 

0 Jo 


rti+\,T  ' w y 

~ tan  L + ~ 


dt 


( 2 -<>) 


8 


Where 


“h  ■ ”iu  * f 'a* 


■o  Jo 


6R  = 6RQ  + 


un  = (rlIi) 
r 


(-vrlt  )dt 
0 11 


r = Tq  + 6R  ( 6R  measured  from  sea  level) 


Xjj  = XE  cos  i|)TT  + Yr  sin  i|r 


II  E 


( 2 -10) 


Yjj  = -Xj,  sin  + Yj,  cos  if) 


z.i  ’ zr. 


And  the  velocity  and  the  gravity  terms  are  obtained  as  follows 


Re  = c„  cos  ^T( 


Vr  = Ejj  sin  ipj 


( 2 -11) 


V'E  V‘II 


C R0  _£  2 ~ |jog  (r q 4 (3  cos2  y - 1) 

/ Ir 


(un)2 


Next,  consider  the  suborbital  phase 


L’  =11  f + v f.  + C. 
Ln  ^EA*!  EA*2  E 3 


VD  = UEAml  + VEAm2  + V,Em3 


( 2-12) 


UEAnl  + VEAn2  + VEn3 


where  the  effective  airspeeds  l'p^  and  Vp^  arc  defined  in  terms  of 
the  nortli  (V.'u)  and  east  (b'v)  component  of  the  wind  shears.  No 
updraf t/downdraf t are  included  in  this  model.  Values  of  V.  and 
VT  are  stored  at  5,000  foot  increments,  from  0 to  100,000  feet 
altitude  where  this  simulation  starts. 


TA 


= iv  - v: 


v 


Ext 


= Vp  - (r  cos  L)we  - W 


( 2 -13) 


2 .3  Rotational  Equations:  The  moments  actine  on  the  vehicle 

are  estimated  about  the  body  axes  as  shown  below.  Let  e be  the 

r 

reference  length  of  the  vehicle.  Let  Jp,  Jp  and  be  the  dimension- 
less numbers  of  roll,  pitch  and  yaw  jet  firings  and  L and  Ly  be 
the  rolling  moments  due  to  one  roll  and  one  yaw  reaction  jet 
firing,  respectively . Then 


r - npv 

L1J  " 2V~ 


[C.  (P„)  + c 


(R  )]  + qCb  { [C 

n ^ 


+ c. 


( 6e ) 


oe 


+ (C:, 


+ C, 


) ( 63L  - 25°)  + C. 


( 6EL  - 00°) 


’6SD., 


'6SD, 


'fiStt. 


( 2 -14) 


+ C 


(6G) 


1 (P)  + C.  (6a)  + C.  (fir)  } 
'6a  *6r 


+ ia-£r  1 (znoii  - z)  * Yn  - y * Zb)  + (V  JI 


Mn  = <crn,  > «*>  + nsc  [C 


C.  c 

( 60 Bright  + (fio)  li.it 


nu 


( 2 -15) 


+ c: 


n6SR  (6SIJ  ' 6SBKOM>  + CnrT  t,  <*Br>  + Cru.  (6G) 


6LF 


6 G 


+ ACm,]  + m**rt  - <zNOM  - z)  * *‘L  + (xno«:  - x>  • Si;]  + (lp)jp 


10 


NB  " a|^l  fcnp  + CnR  (V]  + ^Sb{  <C n (6e) 


+ (C  + C ) (6SD  - 25° ) + C (6SB  - 60°) 

W <W  °(W 


+ Cn„  (6G)  + CnJ<®>  + Cn . <6a>  + lcn 


( 2 -16) 


+ Cn6r  (5r>>  + (fCB  ' *YB)  (m)  (£r)  + (LY)jY 

6 


Appropriate  react  ion-jet  terms  can  be  added  to  the  MR  and  NR 

equations  if  later  vehicles  are  so  equipped.  The  mean  aerodynamic 
chord  length  is  represented  by  c in  the  above  equation.  The 
symbol  x represents  the  x-axis  displacement  of  center-of-gravity 
from  an  arbitrarily  picked  reference  point.  Similar  definitions  apply  to  y and  z. 


* = (cgv  - cgv  ) + _F  Ax„  + 


>Y  Vbv  / ^ _ 

0 “NOf!  V,’ 


sin  6 


r (Axa  ) Vsin  o 


' max  / 


| sin  0rj 


y = (cgv  - cgY  ) 

* r\  * * 


( 2 -17) 


Z = (ceZn  * cgz  ) + / W F zF 

0 ZN0M  W„  F 


Here,  0p  the  angle  made  by  the  surface  of  the  fuel  with  respect  to 


the  axis 


-1  ZB 

0p  - 90°  + tan  ^ 

D 


( 2-18) 


— . JL.  ... 


Using  the  moments  Lj,,  Mp,and  we  can  obtain  the  vehicle  angular 
accelerations  about  its  own  axes  - - Pp,  Qp  and  Rp 


Pn  = 1 lLn  - lvv(PnRn  - Qn)  + ^Y7.(PnQr,  + 


B yi-  AXYViBIVB  ^B ' XZV  B Uj  “E; 

*XX 


+ IYZ(%  “ PlB  ^ + %RB^TYY  " *ZZ^ 


' T“  <md  * hr  <%nD  * V - Ixz(pu2  - nB2) 


( 2 -19) 


~ IYZ^PBf‘B  " PB^  + ^B^B^ZZ  “ *XX^ 


B j [Np>  - IXY^%  PI5  ) IXZ(fJBnB  “ PB^ 

Z/j 


+ IYZ(PBRB  + QB^  + PBC*B ^ 1 XX  ~ *YY^ 


where  the  inertia  terms  are  defined  as  follows 


!X  = XXX  XXX, 


"®*F  \ I, 


lY  = XW  XYY0  + U— M IyYj 

\ Fmax / 


ZZ0  ZZF 

■ max  , 


Ij!z0  + 

i max  / 


O W, 


ixz_  + L_F_  i, 


( 2 -20) 


..JL 


+ 


Stability-axis  angles,  namely  angle  of  attack  a and  angle  of 
sideslip  3 can  be  computed  as 


tan-1  W- 


iU. 


B 


( 2-24) 


B - tan"1  /VB  cos  »' 


u 


B 


Flight  path  angle 


y “ 0 - a cos  <{>  - 6 sin 


( 2-25) 


Let  be  the  given  angle  between  true  north  and  the  vector; 

then  down  range  and  cross  range  can  be  computed  as: 


FvD  = r0  l(L  " COS  + (x  - xnl  cos  L sin  'M 


KR 


R 


(2-26) 


Rx  = r0  [(L  - L„)  sin  i|ip  + (X  - Xn)  cos  L cos  (|/R1 


'R 


R 


Equivalent  airspeed  REAS  and  calibrated  airspeed  LCAS 


KEAS 


Ik 


o 


3600 

G08b" 


nautical  miles 


KCAS  = V1  . (KEAS) 


( 2 -27) 


The  quantities  v,  and  vn  are  defined  by  the  following  empirical  formulae: 

JL  4- 


»1=f(x)  ; x =.  Mach  number 


( 2-28) 


f ( x ) ; x « (.001511)  (KCAS) 


Where  f(x)  assumes  the  following  form: 
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Geometric  altitude  h and  geopotential  altitude: 


Air  density  can  be  approximated  a 
upon  standard  atmospheric  data: 


a function  of  altitude,  based 


3 . MAIN  PROGRAM 

No  efforts  are  spared  to  organize  the  main  program  in  a self- 
explanatory  manner.  Parameters  and  data  are  grouped  together 
according  to  their  roles  and  the  omnipresence  of  the  comment  cards 
should  guide  the  reader  through  the  program.  It  is  suggested 
that  subsequent  users  do  not  deviate  from  the  existing  logical 
organization  when  making  modifications/additions,  so  not  to 
defeat  the  purposes.  In  order  to  simplify  the  loading  process, 
all  parameters  and  data  (including  aerodynamic  coefficients)  are 
stored  in  the  main  program  , thus  avoiding  the  use  of  COMMON  blocks. 

As  indicated  by  figure  3 -1,  the  main  program  is  composed  of 
three  parts:  the  initialization,  the  Hi-speed  loop  (starting  from 
statement  47),  and  the  EASE  subroutine  designed  to  display  and 
modify  any  parameter  of  the  main  program.  EASE  was  written  to 
allow  interchange  of  parameters  without  the  use  of  COMMON  blocks. 

3.1  Initialization:  Most  data  parameters  and  initial  conditions 
are  on  punched  cards  to  allow  maximal  flexibility  on  flight  conditions 
and  vehicle  configurations.  All  aerodynamic  derivatives  are  stored 
in  the  form  of  lG-bit  normalized  fixed  point  fractions.  The  detailed 
discussion  of  their  handling  will  be  covered  under  the  heading  of 
"Function  Generation".  During  the  initialization  phase  the  program 
will  ask  the  user  if  he  is  through  with  changing  parameters  before 
entering  into  the  next  phase,  by  typing  the  question  "ARE  FURTHER 
CHANGES  DESIRED?"  A "YES"  answer  will  automatically  call  in  the 
EASE  program.  A "NO"  answer  will  set  in  motion  the  high-speed 
loop. 

Not  all  parameters  can  be  independently  initialized.  For 
instance, Vp  is  called  an  "intermediate"  variable  because  it  is  a 


function  of  ij>p  and  Up, 


Thus  the  initial  value  of  Vp  is  determined 


once  the  initial  conditions  of  and  Un  are  specified,  and  vice  versa. 


3.2  Hi-Specd  Loop:  This  loop  is  actuated  by  an  interrupt 
which  occurs  at  synchronous,  regular  intervals  At  whose  duration 
was  set  a-priori.  During  the  time  At,  the  computer  must  finish 
solving  one  frame  of  the  dynamic  equations  before  the  occurrence 
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f igure  3-1  (Continued) 


On  Next  l’age 

18 


Figure  3-1  (Continued) 


Integration  to  obtain  body  rate 
of  rotation  P , Q , R 


/ICN 

JlXor  HOLf 
mode> 


Xic  \ 

or  HOLET 
\mode 


Calculate  Euler-angle 


rates  9 , cp  and  ^ 


Want 


3 -DOF 
\onlyj> 


v|/  = cp  = 0 


Integration  to  obtain  9 , cp  and  (j/ 


'WanK 
5 -DOF 
.only?'' 


Set  accelerations  UB  , 
VB  and  WB  to  zero 


To  Be  Continued  On  Next  Page 
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Figure  3 - l (Continued) 


1)  Coordinate  transformation  to  obtain  earth- 

axis  accelerations  Xg  , Y^,  , 

2)  Coordinate  resolution  to  obtain  H- frame 

accelerations  Xjj  , , ZH 

3)  Computation  of  gravity 


• ' 

Calculate  H-frame  derivatives  rUE  , 


'W  ant\ 
5 -DOF 
vonly?/ 


/IC\ 
or  HO  LIT' 
\mode  l/ 


1) Integrate  to  obtain  rU^  , ^ , 

WH  , L , A and 

2)  Calculate  r , Ujj  separately 


1)  Compute  down  range,  cross  range  & altitude 

2)  Compute  vehicle  ground  speeds  UE  , VE  and  W£ 


WIND  = 0 


^atmosphere  — ■ 


Wind  model 
subroutine 


420' 


Figure  3-1  : Continued 


1)  Compute  vehicle  airspeed  components  U£A  : veA  ’ h i 

2)  Generate  body-axis  velocities  UB  ; Vg  ; WB 


Evaluate  a and  0 


/ Want 
5 DOF  only. 


-H  True  airspeed  V=constant 


1)  Calculate  p , q , equivalent  airspeed,  calibrated 
airspeed,  Mach. 

2)  Digital  FCS  processing 

3)  Examine  pulse  flags  (DEPULSE,  DAPULSE 
and  DRPULSE  ) 

4)  Apply  pulse  to  the  desired  surface 

5)  Digital  FCS  Processing  (bypassed  if  any  surface 
is  pulsed ) 

6)  Examine  surface  flags  ( FCSDE  , FCSDA  , FCSDR, 
FCSSB  , FCSBF  ) 

7)  Set  surfaces  that  are  Ragged  equal  to  their  IC  values 
( DELIC  , DERIC  , DRIC  , DSBIC  , DBFIC  ) 

8)  Preparation  of  output  


To  be  continued  on  next  page 
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of  the  next  interrupt.  In  order  to  achieve  the  real-time  (or  time- 
critical)  goal,  great  care  was  made  to  minimize  the  execution  time. 
Even  so,  a number  of  calculations  have  to  be  delegated  to  the  analog 
computer  such  as 

- Sin/cos  functions  of  angles  0,  <J>,  4>,  a for  instrumentation 

- generation  of  test  pulses 

In  fact,  if  more  equations  are  to  be  added  into  the  dynamic  model, 
some  other  existing  equations  must  be  deleted  from  the  high-speed 
loop.  An  alternative  would  be  solving  the  deleted  equations  on 
another  digital  or  analog  processor.  The  dynamic  model  is  given 
in  Chapter  II,  thus  not  repeated  here.  During  run  time, the  operator 
can  intervene  by  throwing  a number  of  switches  which  in  turn 
activiate  a number  of  interrupts  in  the  digital  computer.  The 
switches  are  situated  in  the  general  cockpit  area  and  they  consist  of 
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- IC  Switch:  sets  the  program  back  to  Initial  Condition 
mode.  All  information  concerning  the  preceding  run  is 
lost. 

- HOLD  Switch:  also  stops  execution  but  all  state  variable 
values  are  frozen  at  their  respective  values  reached 
prior  to  activation  of  the  switch.  This  mode  allows 

the  user  to  investigate  the  status  of  his  program  in 
post-mortem  fashion. 

• • 

- 3 DOF  Switch:  sets  the  yaw  rate  and  the  roll  rate  0 

to  zero,  the  vehicle  can  then  be  studied  entirely  in 
the  pitch  plane,  having  three  degrees-of-freedom:  x, 
z,  and  6.  This  mode  is  useful  in  planning  and  perfor- 
mance analysis  such  as  plotting  the  gliding  envelope. 

In  this  case,  we  are  interested  only  in  getting  from 
one  point  to  another  point. 

- 5 DOF  Switch:  imposes  a nearly  constant  altitude,  constant 
total-velocity  flight  regime  . This  mode  offers  a con- 
venient way  to  study  vehicle  handling  quality  and 
maneuverability.  To  different  people,  5-DOF  simulation 
means  slightly  different  things.  Thus  we  derive  the  5-DOF 
situation  used  in  this  simulation  for  reference  purposes. 

Consider  the  acceleration  equations  in  body-axes,  excluding  the 
thrust  and  aerodynamic  force  on  the  surfaces. 


g 
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WB 

“VB 
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+ Slp.Sfy) 
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- c4>.s<t>) 
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Inertial-axis  — *»  body-axis  transformation 


3 -1) 


Where  cip  and  sip  are  the  abbreviations  for  cos  ip  and  sin  ip  respectively . 
For  5 DOF,  assume  that  Ug  = Vg  = WQ  = 0 and  solve  for  the  acceleration 
component  XQ 


X 


B 


ZBWB  + YBVB 
U„ 


+ g 


sin0 


cos0  (VBsin<|>  - WgCost})) 


( 3 -2) 


Equation  ( 3 -2)  is  used  in  lieu  of  XQ  described  in  equations  ( 2 -4) . 

The  only  other  changes  in  the  5-DOF  case  are 

d ( Ar)  _ 

Ht 0 

( 3 -3) 

and  V = VIC  , a constant 


The  closed -loop  nature  of  the  hi-speed  loop  is  illustrated  by  the 
block  diagram  of  figure  3 -2.  This  block  diagram  depicts  the  main 
stream  of  the  calculations  and  how  different  groups  of  calculations 
are  related  to  each  other  functionally. 
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Figure  3-2:  Block  Diagram 
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H - Axis  Velocities 
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4.  FLIGHT  CONTROL  SYSTEM  (FCS) 

4.1  Description:  In  this  simulation  we  investigate  the  Orbiter 
handling  qualities  under  manual  control  of  a pilot.  There  is  no 
provision  for  an  autopilot  in  this  program.  From  the  time  the  Orbiter 
enters  the  atmosphere  at  a nominal  altitude  of  400,000ft  to  touchdown 

it  is  not  power-assisted  except  for  the  Reaction-jet  Control  System  (RCS). 
Thus  the  importance  of  the  FCS  can't  be  too  much  exagerated. 

The  Orbiter  is  equipped  with  an  all-digital  fly-by-wire  FCS  which 
insures  aerodynamic  stability  with  the  help  of  five  on-board  computers 
(four  redundant  primary  and  one  secondary).  The  computers,  in  turn,  feed 
signals  to  the  actuators  of  the  control  surfaces  . A picture  of  the 
control  surfaces  can  be  found  on  page  22  of  reference  9.2.2. 

a / Elevons  ( Elevator-aileron  combination):  They  are  full-span 
in  construction,  used  for  affecting  both  pitch  and  roll.  This  program  does 
not  exercise  the  inboard  and  outboard  elevons  separately. 

b / Rudder  / Speedbrakes:  The  conventional  rudder  is  made  up 
of  two  symmetrical  halves  which,  when  deployed,  constitute  the  speed 
brakes. 

c / Body  flap  is  found  on  the  very  end  of  the  fuselage  below  the 
main  engines  and  serves  mainly  as  augmentation  to  longitudinal  trim. 

Whenever  the  control  surfaces  prove  to  be  ineffective  (high  altitude,  low 
dynamic  pressures)  pitch  and  roll  reaction-jets  are  fired  to  guaranty  the 
orbiter' s aerodynamic  stability.  The  yaw  reaction  jets  are  used  more 
often,  including  the  cases  just  mentioned  above. 

4.2  FCS  Block  Diagrams:  The  pitch,  roll  and  yaw  modes  of  the  FCS 
are  depicted  by  the  block  diagrams  on  figures  4_1  , 4-2  and  4-3  in 

that  order.  Let  us  discuss  one  of  the  modes  in  detail,  say  the  pitch  mode 
of  figure  4-1  . A command  from  the  RHC  ( Rotational  Hand  Command  ) 
is  fed  through  the  deadband  and  pitch  shaping  networks  then  compared 
with  the  current  Orbiter  performance  ( a function  of  the  pitch  rate  Q ) 
to  generate  a pitch  rate  error  signal  called  DPJ.  Because  of  the 
coupling  effects  between  the  axes,  the  nose  tends  to  dip  during  a turn. 

To  compensate  for  this,  we  want  an  additional  negative  elevon  deflection 
( trailing  edge  up  , by  convention  ) and  the  term  R*tancpdoes  exactly 
this  by  contributing  a negative  influence  to  the  elevator  command  6 

CMD  . 


When  the  speedbrakes  are  deployed,they  force  the  nose  upward  by 
contributing  a moment  around  the  Orbiter's  center  of  gravity  ( CG  ). 

Note  how  the  speedbrake  increment  ( DSBPC  ) signal  is  used  to  negate 
the  unwanted  pitch-up  by  a direct  reduction  in  the  pitch-position  trim. 

The  FSC  roll  channel  is  represented  on  figure  4-2.  The  switches 
EARLY  and  LATE  correspond  to  MACH  >1.5  and  MACH  ^1.5  respectively. 
Notice  the  roll-rate  command  main  path.  Starting  with  the  roll  stick 
position,  the  signal  goes  through  deadband,  shaping  and  first-order  lag 
circuits  before  being  added  with  the  integral  of  the  roll-rate  trim  term 
at  SUM21.  From  there, the  control  signal  is  sent  through  SUM  24  where 
the  roll -rate  feedback  and  other  coupling  terms  are  subtracted.  All  the 
feedback  signals  taken  together  are  called  PSTAB  (roll-rate  stability). 

SUM24  is  directed  to  a switch  whose  poles  are  marked  EARLY  and  LATE. 

The  EARLY  side  of  the  switch  allows  only  reaction  jets  to  control  the  roll 
channel  because,  at  that  time,  the  Orbiter  is  still  at  high  altitude.  Only 
during  the  LATE  stage  can  the  control  surfaces  perform  effectively. 

The  FCS  yaw  channel  is  illustrated  on  figure  4-3.  The  rudder  pedal 
generates  a signal  which  is  fed  through  the  dead-band  and  shaping  networks. 
Then  it  is  summed  with  the  integrated  yaw-trim  signal.  Before  reaching 
the  actuators  proper, this  signal  is  further  modified  by  Mach  number  and 
dynamic  pressure  q.  An  arrangement  similar  to  the  roll-axis  channel  is 
implemented  here  to  allow  the  yaw  jets  UZCMD  to  take  over  at  high  altitude 
(EARLY).  Only  at  low  altitude  (LATE)  can  the  rudder  work  the  way  it  works 
on  an  airplane. 

Figure  4-4  describes  how  the  Orbiter's  control  surfaces  are  exercised 
by  different  actuators. 

4.3  FCS  Computer  Implementation:  In  the  block  diagrams  the  transfer 
functions  and  filters  are  expressed  in  the  Laplace  s-domain.  We  could 
have  used  the  s-plane  if  an  analog  computer  were  used.  Since  the  FCS 
program  is  done  digitally  it  is  advantageous  to  map  the  transfer  functions 
from  the  s-plane  to  another  complex  plane,  the  z-plane.  In  simple  terms, 
z can  be  defined  as  z=esT  where  T is  the  sampling  period.  A comprehensive 
treatment  of  infinite-impulse-response  digital  filters  can  be  found  in  reference 
9.1.6  , chapter  4 . 

It  is  a known  fact  that  a stable  analog  filter  might  or  might  not  map  into 
a stable  digital  filter  if  we  just  go  about  replacing  the  differentials  by  finite 
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FIGURE  4-1=  FCS  PITCH  CHANNEL 


FIGURE  4-2:  FCS  ROLL  CHANNEL 


DEADBAND 


< MACH  < I 0 


differences  (backward,  forwarder  central  differences).  A more  elaborate 
process  is  used  in  this  program  to  prevent  computational  instability  and 
aliasing  problems.  It  is  called  Bilinear  transformation  (simple  conformal 
mapping)  and  it  is  represented  by  the  following  transformation: 


- 2 ( 1 - z"1  ) 
S T ( 1 + z-1  ) 


(4-1) 


Consider  a first-order  transfer  function  as  an  example.  Given 


H(s) 


10 

s + 10 


Substituting  s for  z as  dictated  by  equation  (4-1) 


H(z) 


10_ 

2_  ( 1 - z~j) 
T ( 1 + z"1) 


+ 10 


/ 

H(z) 


10T  + lQTz~ 1 
2 - 2z_1  + 10T  + lOTz  1 


(4-2) 


(4-3) 


In  this  program  the  frame  time  is  T=.04  second.  Any  first-order  transfer 
function  can  be,  for  our  programming  purposes,  reduced  to  the  form  shown 
below: 


x 

X 


G1  + G2  * z~ 1 
1 + G3  * z'1 


(4-4) 


Similarly  we  put  all  second-order  transfer  functions  in  the  following  form 
prior  to  digital  coding: 


Gi  + G2  *z-1  + G3*Z~2 


1 + G4*z_I  + G5*z*'2 


(4-5) 


Thus  the  z -transform  of  equation  (4-2)  finally  becomes,  after  introducing 
the  proper  value  for  T into  equation  (4-3) 


H(z)  = 


.1667  + . 1 667  z ~ 1 
1 - .6667  z"1 


(4-6) 


Computer  mechanization  of  equations  (4-4)  and  (4-5)  are  done  in  the  form 
of  function  subroutines  called  FILT1  and  FILT2  respectively.  FILT1  and 
FILT2  stand  for  first-order  filter  and  second-order  filter.  For  FILT1 
the  expression  comes  directly  from  equation  (4-4) 


Gl*Xn  + XOLD 


XOLD  = G2*Xn-l  - G3*Yn-l 


For  FILT2,  equation  (4-5)  can  be  put  into  the  format 


(4-7) 


G-,*X  + X 

1 n nodel 


(4-8) 


-here  Xnodel  - G2*Xn.1  - G4*Vl  + Xnode2 

xnode2  = G3*Xn-2  “ G5*Yn-2 

Refer  to  Table  4-1  to  correlate  an  s-domain  transfer  function  to  its 
z-domain  counterpart.  Besides  the  filters,  the  hysteresis  transfer 
function  also  deserves  a brief  mention.  Figure  4-5  depicts  the  flowchart 
of  the  HYSTER  subroutine  and  it  is  self-explanatory. 

A detailed  flowchart  of  the  Shuttle  Flight  Control  System  (SHTLFCS 
subroutine)  is  contained  in  Figures  4-6.  The  actual  listing  of  SHTLFCS 
can  be  found  in  Appendix  D. 
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TABLE  4-1  : Transfer  Functions 


o 


- Generate  SIN(ALPHA)  Schedule  Breakout 

- Generate  COS( ALPHA)  , SIN (PHI) 


COS  (PHI)  , COS(THETA) 

- Generate  roll,  pitch  & yaw  trim  breakouts 

for  DETMPAN , DETMRHC , DATMPAN 
DATMRHC , DRT  and  DBFMAN 

- Generate  values  of  Mach-dependent  function 

schedule  for  GPS,  GJET,  ALPHMAX, 
ALPHMIN,  KPIT , GTRE , GDSB,  KDAMPAR 
KDAMLIN,  GRS , GDAC,  GP,  GNYDRM, 
GRAY,  KDRC , GTRR  " ' 
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To  be  continued  | on  next  page 
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To  be  continued  ^on  next  page 
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EASE  PACKAGE 


5 . 


EASE  is  an  interactive  package  that  allows  the  researcher  easy 
access  to  the  Simulation  program  in  the  format  that  is  natural  to 
him.  To  qualify  as  a general-purpose  package,  the  program  must 
readily  accept  different  vehicle  configurations  and  flight  condi- 
tions which  vary  from  case  to  case.  The  list  of  parameters  and 
variables  can  be  found  in  Chapter  VI  together  with  their  units 
and  computer  symbols.  Users  of  EASE  are  required  to  follow  two 
simple  formats  corresponding  to  an  interrogation  and  an  input. 

EASE  can  be  accessed  through  the  Main  program  or  by  setting  Control 
Switch  12  to  .TRUE. 

a/  Examples  of  Interrogation  Mode 

THETA?  What  is  the  current  value  of  0?  (in  F25.10  format) 

*12AF?  What  is  the  integer  content  of  address  12AF? 

(In  decimal  format  125).  The  asterisk  (*) 
specifies  the  integer  mode. 

43600?  What  is  the  floating-point  content  of  address 
43600? 

b / Examples  of  Input  Mode 

OMEGA  = .000072923  Replace  the  content  of  OMEGA  with 

. 72923E-4 . 

*17777  = 50  F ..;,:ce  the  content  of  cell  17777 

wi.n  50  decimal 

1FF2C  = .25  , Put  .25  (in  floating-point  format) 

into  cell  1FF2C. 

cl  To  exit  EASE,  type  an  exclamation  mark  "I"  followed  by 
a carriage  return  (C/R). 

For  unusally  large  numbers  or  unusually  small  numbers,  check  the 
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accompanying  listing  to  avoid  truncation  of  the  input  data.  For 
example,  the  I/O  floating-point  number  is  determined  at  format 
F25.10.  Any  digits  beyond  10  places  to  the  right  of  the  decimal 
point  will  be  ignored.  The  user  should  also  safeguard  against 
entering  a value  for  a parameter  that  is  larger  than  the  largest 
allowable  for  the  parameter.  For  example,  the  pitch  rate  Q0 
has  a range  of  +2  rad/sec  and  is  scaled  at  [Qb/2]  when  appearing 
to  the  Flight  Control  System.  If  the  user  of  EASE  inadvertently 
setSQB  = 3,  the  output  of  the  DAC  (digital-to-analog  converter) 
corresponding  to  Qq  is  completely  wrong. 

D 

When  a wrong  symbol  was  typed,  EASE  will  respond  with  "NO 
SUCH  SYMBOL  IS  DEFINED  TO  EASE".  There  are  other  interactive 
features  built  into  EASE  to  help  the  users  in  case  of  mistakes. 
Here  are  a few  common  typing  mistakes  to  serve  as  examples: 

a/  OMEGA  = 2.6F19 

EASE  will  come  up  with  the  message  "MEMORY  DATA  MUST  BE 
A DECIMAL  NUMBER"  because  there  is  no  provision  for  EASE 
to  accept  hexadecimal  numbers. 

b / *RHO?  (even  though  p was  declared  a floating-point  in 
the  main  program) . The  resulting  message  is 
"WRONG  INPUT  FORMAT". 

cl  BYTE  = 377  (even  though  BYTE  was  declared  as  8 bits) 

In  that  case  255  £ BYTE  < 0.  EASE  will 
alert  the  user  with  the  warning  "DATA 
OUTSIDE  RANGE". 

The  simplified  flowchart  for  the  EASE  package  is  contained  in 
figure  5 -1  for  reference  purposes. 
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1)  LETTER-. TRUE. 

2)  J=1 

3)  TEMP-"  " 


Figure  5-1  (continued),  Letter 

Processing 


SYMBOL  irTVW 
LIST(K) 


Display  "NO  SUCH  SYMROL  IS 
DEFINED  TO  EASE" 


r~ 


Figure  5-1  (continued),  Number  Address 


5 DIGITS" 


Figure  5-1  (continued).  Number  Address 


Output  Data  To  CRT  Display 


6.  FUNCTION  GENERATION  PACKAGE  (FGP) 


Fast  function  generation  constitute  the  main  objective  of  this 
software  package  so  it  can  successfully  support  a time-critical 
simulation.  The  package  is  Fortran  compatible  making  it  adaptable 
elsewhere  if  needed.  On  the  other  hand,  because  of  speed  consideration, 
some  part  of  this  program  were  written  in  Assemble  language.  Physi- 
cally the  FGP  is  embedded  into  the  Main  program  and  considered  an 
integral  part  thereof.  This  leature  also  helps  speed  up  the  execution 
time.  Aerodynamic  data  for  this  particular  simulation  are  composed  of 
about  40  functions  having  from  one  to  three  arguments.  Overall,  there 
exist  more  than  4000  data  points  to  be  handled  by  the  FGP.  If  there  is 
a shortage  of  memory  space  each  dal  a point  could  be  converted  into  a 
fixed-point  quantity,  normalized  then  made  to  fit  into  a 16-bit  half-word 
location.  But  in  that  case,  the  FGP  must  be  greatly  modified. 

V/e  are  provided  with  data  in  the  form  of  punched  cards. 

The  card  image  of  a representative  two-argument  function  is  shown 
below.  In  this  case  the  function  is  the  pitching-moment  coeffi- 
cient (or  derivative)  bias,  CM^  (a,  Mach),  which  is  a function  of 
the  angle  of  attack,  and  the  velocity. 


Ob  45 
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Note 

that  the  first  card 

(called 

the  Header  card)  gives 

the  deri- 

vative  name  (CMO)  and  the 

number 

of  data  points  (45), 

Derivative 

names 

must  start  with  the 

letter 

"C"  or  they  will  be  rejected  . 

There 

is  no  mention  of  how  many 

arguments  are  involved. 

Only 

when 

we  make  use  of  C"n  in  the  Main  program  that  the  type  of 

argument  and  its  breakpoints  appear  explicitly.  Tor  this  example, 
the  two  arguments  are  arranged  in  the  order  ) ; 


' - - — 
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CMq(  ci2  , Mj)  ; 

CMq  (a2 » M2)  ; CMj^o^jlV^)  ; . 


;CM0(a5,M1)  ; CM0(a1,M2)  ; 

. until  CMq  (a^.  Mg)  . 

The  FGP  is  made  up  of  three  separate  sub-programs,  each  doing  a 
specified  task.  They  are:  DATASTORE,  POINT  and  DERIVE4.  Their 
particular  roles  are  described  below. 

6.  1 DATASTORE  Program  : This  program  performs  the  following  tasks: 

- Read  aerodynamic  data  from  the  card  file  as  described  earlier. 

- Check  the  validity  of  data  including  the  name. 

- Rearrange  the  data  in  a predetermined  order. 

- When  an  end-of-file  (EOF  ) card  is  encountered,  write  the  whole 
data  array  on  a disc  file. 

The  flowchart  for  DATASTORE  can  be  found  on  figure  6-1.  We  can  analyze 
DATASTORE  program  by  tracing  through  the  handling  of  one  stability 
derivative,  say  CMO  again.  When  MAIN  program  calls  DATASTORE,  it 
reads  and  inspects  the  first  card,  also  called  the  Header  card.  That  card 
should  contain  an  8-character  derivative  name  followed  by  a blank  and 
followed  by  a 6-digit  integers  ( hlanks  $ included).  For  example 


Derivative  Columns  #16  thru  #80 

name  reserved  for  comments 

f \ t s 

CMOMMM£Mtf4  5,COMMENTS.  . . . 

Number  of 
data  points 

One  header  card  is  required  per  derivative  by  DATASTORE.  Anytime 
an  EOF  card  is  encountered  in  place  of  a header  card,  control  is  reverted 
to  the  MAIN  program.  Thus  the  very  last  card  must  be  an  EOF  or  the 
program  will  hang  up. 
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Since  all  derivative  names  stored  start  with  the  letter  C , 
a test  is  made  of  the  first  character  to  eliminate  any  unnecessary 
search  of  the  Data  List.  Only  those  names  previously  defined  in 
the  Main  program's  Data  Table  ( such  as  CMq,  C . , C , ...  ) 

9 6a 

are  considered  legitimate.  If  the  first  letter  of  the  header  card 
is  a "C  ",  control  is  passed  to  the  "name  list  search  portion  of 
DATASTORE  program  (beginning  with  statement  #50).  In  the 
event  the  first  character  of  the  header  card  is  not  a C , 
DATASTORE  allows  the  user  the  choice  of  aborting  the  run 
( returning  to  Main  program  ) or  the  choice  of  correcting  the  typo- 
graphical mistake  by  inputting  a new  name  on  the  CRT. 


After  the  derivative  name  and  address  have  been  found 
in  the  tables,  DATASTORE  returns  to  the  card  reader  and  reads 
the  information  on  the  cards  stacked  behind  the  header  card. 

The  newly  arrived  data  points  are  stored  in  the  temporary  array 
called  DATA(I)  , ( 1 = 1,  LENGTH  ).  In  our  example  LENGTH  = 45 
as  indicated  on  the  header  card  earlier.  When  all  the  expected  data 
points  are  in  ( an  EOF  card  is  read  ) DATASTORE  terminates 
the  reading  phase  and  all  the  aerodynamic  data  are  written  out  on 
a disc  file.  Notice  that  DATASTORE  is  excuted  only  once, 
during  program  initialization  phase  and  the  disc  file  will  be 
utilized  by  the  Derivative  subroutines  later  on. 

It  should  be  aware  that  DATASTORE  is  a separate  program 
that  must  be  executed  prior  to  all  other  programs  in  the  normal 
sequence  of  operation.  Also  note  that  DATASTORE  program  must 
be  modified  and  rerun  for  any  change  in  the  length  or  table  arrange- 
ment of  the  aerodynamic  data.  These  parameters  must  match  those 
of  the  MAIN  program  to  insure  that  the  data  are  stored  in  the  same 
relative  location  where  they  will  be  used.  Specifically  care  and  atten- 
tion should  be  paid  to  the  following: 

- The  order  and  dimension  of  the  data-table  declaration  statements 

- The  DATA  statements  for  variable  SIZE  and  NUMBER 
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- The  ASCII  variables  in  NAMELIST 


- The  order  of  the  data  table  ACW  (Address  Constant  Word)  array. 
These  parameters  effect  the  proper  ordering  of  the  data  table 
in  the  array  which  is  dumped  to  disc. 

6.2  POINT  Subroutine  : This  subroutine,  which  is  called  by  the  Main 

program,  performs  all  the  preliminary  tasks  prior  to  final  interpo- 
lation. One  call  to  POINT  subroutine  is  required  per  argument  list 
and  a representative  calling  sequence  applicable  to  our  present  example 
might  appear  as  follows: 


'It  lit 

CALL  POINT  (ALPHA  , ALPHAPT(2)  , ALPHAPT(l)  , ALPHAT(l)  , ALPH AT(3)  ) 


CALL  POINT  (MACH  , MACHPT(2)  , ALPHAPT(l) , ALP.HAT(l)  , ALPHAT(3) ) 

1 | t r i 


The  outgoing  arrows  are  added  to  visually  indicate  arguments  needed 
by  the  subroutine  for  processing  the  calls  and,  conversely,  the  in- 
coming arrows  show  that  information  is  being  returned  to  the  Main 
program.  Two  calls  are  made  to  POINT  subroutine  because  in  our 
current  example  , CMO  is  a function  of  two  arguments,  a and  Mach. 
Let  us  inspect  one  of  the  calls  above. 

ALPHA  = an  argument,  the  present  value  of  a 

ALPHAT(3)  = starting  point  of  a list  of  breakpoints  associated  with 
the  independent  variable  a.  This  list  consist  of 
floating-point  values. 

ALPHAPT(l)  = returns  to  Main  program  as  an  integer  value 

(called  "pointer")  indicating  the  smallest  breakpoint 
adjacent  to  the  argument  a,  say  . This  information 
will  be  used  during  subsequent  interpolation. 


w 


I T' 


/ ; _ 

^Subroutine  POINT 


( 


1 


Acquire  argument  value 
Load  current  POINTER 




* 


j Update  BREAKPOINT 
— j corresponding  to 
I current  POINTER 


^BRKPTl  > ARG^ompar>^BRKpT)<  A 
( 70;< <v(BRKPT)  to  > — * 30V 


POINTER  \.o>— — 
POINTER  >-LIM  \ _LIM  " 


POINTER < MAXBP.KPT 


late  BRKPT 


POINTER=POINTED+ 1 

•flompare 
(BRKPT+1)  tc 

(BRKPT+1)=ARG  \ ARG  / (BRKPT+1) < ARC 


(BRKPT+ 1)>(  ARC 


DIVI  = (BRKPT+ 1 ) - (BRKTT) 
FRAC  = ARC- (BRKPT)  /DIVL 


Notes : Abbreviations  used: 

BRKPT  = Breakpoint  number 
-LIM  = 1,  the  first  brkpt 
MAXBRKPT=  no.  of  brkpts  of 
an  argument. 

ARG=value  of  the  argument 
(BRKPT) =the  value  of  the  argument 
corresponding  to  BRKPT 
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RETURN  to  "\ 
ain  program/ 


Figure  6-2  : Subroutine  POINT 


I 


ALPHAT(l)  = returns  to  Main  program  representing  the 
distance,  FRAC  = ( a-  ) / ( a^+1  - Oj  ) 

Finally,  ALPHAPT(2)  contains  the  number  of  breakpoints 

associated  with  the  argument.  In  this  case  there 
are  9 breakpoints  as  declared  by  a DATA  state- 
ment in  the  Main  program. 

When  subroutine  POINT  is  entered,  a test  is  made  to  compare  the 
argument  ( a in  this  instance  ) against  the  extrema  of  the  argument 
values  (a-j  and  a 9 here).  Should  a<aj  then  we  set  the  pointer 

ALPHAPT(l)  to  1.  Conversely  whenever  a>a9  we  set  ALPHAPT(l) 
to  9.  A simplified  flowchart  of  the  POINT  subroutine  can  be  found 
in  figure  6-2. 

6.3  Derive  Subroutine:  This  section  was  written  using  the  FUNCTION 
approach  instead  of  the  standard  SUBROUTINE  approach  in  order 
to  reduce  program  length  and  minimize  the  transfer  of  arguments 
between  different  programs.  There  exist  four  DERIVE  functions  : 
DERIVE1  , DERIVE2  , DERIVE3  and  DERIVE4 . Together  they  can 
handle  any  function  having  from  one  to  four  arguments.  A sample 
use  of  DERIVE  1 is  as  follows: 

FUNCT  = DERIVE1  (DAT  AT  ABLE  , MACHT  ) 

where  DAT  AT  ABLE  is  the  starting  location  of  the  memory  table 
containing  the  data  which  were  were  placed  there  earlier  by  the 
DATASTORE  subroutine.  MACHT  is  an  argument  table  containing 
FRAC  in  the  first  word  (see  figure  6-2  for  the  significance  of  P'RAC), 
the  breakpoint  pointer  in  the  next  half-word  and  the  number  of  Mach 

breakpoints  in  the  Last  half-word. 


MACHT 

FRAC 

| POINTER 

MACHT  LENGTH 
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For  our  existing  example,  the  pitching  derivative  bias  CMO 
appears  at  last,  on  the  right-hand  side  of  an  assignment  statement 

CMO  = DERIVE2  ( CMOT  , ALPHAT  , MACHT  ) 


7.  REPRESENTATIVE  SOLUTIONS 


No  simulation  can  be  completely  trusted  unless  it  has  undergone 
preliminary  STATIC  and  DYNAMIC  CHECKS.  It  is  imperative  that 
such  checks  be  carried  out  before  the  very  first  production  run  and 
especially  if  the  simulation  is  just  started  up  after  some  lay-off 
time.  Check  conditions  need  not  be  unique;  any  reasonable  set  of 
conditions  will  suffice.  However,  some  considerations  should  be  given 
to  checking  conditions  commonly  encountered  in  the  flight  regime. 

7 . 1 General  Procedures 

The  STATIC  CHECK  for  the  rotational  equations  could  be 
performed  as  follows:  an  initial  flight  condition  is  established  by 
fixing  Mach  number,  dynamic  pressure  and  angle  of  attack  at  pre- 
selected values.  Input  variables,  angular  rates,  Euler  angles,  control 
surfaces,  etc  . , to  each  of  the  equations  are  then  perturbed  one  at  a 
time  or  in  pairs  to  excite  particular  terms  of  the  equations  to  be 
checked.  The  outputs  should  be  scrutinized  for  trends  and  for  possible 
wrong  signs. 

By  choice  of  certain  parameters  we  can  isolate  each  portion 
of  an  equation  to  see  if  it  gives  the  correct  answer.  A case  in  point 
is  illustrated  by  considering 


qSb  [" 

lx  <- 

(C,  PR+C.  R_)  — 
*p  B % B 2 V 

+ C*s 

rXZ 

flY  - lZ 

+ 

IJ  <rb+  pb%>  + 

{ >x 

If  we 

set  Pg 

= Qg  = 0 and  B = 6g  = 

6r  = XXZ 

+ C . S + C . 6a  + C,  6r 


'6a 


'6r 


(7-1) 


= 0 we  can  observe  the 


B 


HRIS>. 


By  resetting  different  parameters,  we  can  gradually  examine  the  influence 
of  pach  of  the  terms  in  the  PR  equation,  remembering  to  take  into 
account  the  sign  of  each  term. 
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A similar  procedure  is  used  to  check  the  aerodynamic  coef- 
ficient equations  for  CN  ' Cc  - Ct  ‘ Cm  > cn  and  Cy  as  functions 
of  a , 9 , Mach  number  and  control  surface  configuration. 

DYNAMIC  CHECKS  should  be  attempted  only  after  the  user 
is  fully  satisfied  with  the  static  tests.  Again  a few  checks  are  sug- 
gested herein  for  illustration  purposes.  They  are  neither  unique  nor 
exhaustive.  To  test  the  dynamic  behavior  of  the  rotational  equations, 
a stable  flight  condition  is  picked  and  initialized  by  setting  a,  3 , 

66  TRIM  ’ 6rTRIM  ' 6aTRIM  ’ V and  h ^ TRIM  conditions  are  selected 

so  that  the  initial  moments  are  zero,  i.e.,  a and  0 remain  unchanged). 

Each  surface  (&e  , 6r  , 6&  ) is  pulsed  separately  and  the  resulting 

• • • 

aerodynamic  responses  on  PB  • PB  ' ^B  * Qb  » RB  • RB  * a » 6 » 
acceleration  terms  NZ  & NY  , 9,  cp  and  ♦ are  recorded  on  a strip- 
chart  and  compared  to  known  correct  responses.  The  suggested  shape 
of  the  bidirectional  ramp  shown  below  could  be  used  for  forcing 
function  (deflection  varies  from  ±5°  for  subsonic  region  to  ±20°  for 
supersonic  region ) 


time  (sec) 


It  is  recommended  that  this  test  be  conducted  on  a daily  basis  to 
verify  the  operation  worthiness  of  the  equipment.  DYNAMIC  CHECKS 
should  be  carried  out  on  the  translational  (orbital)  equations  also. 

If  all  aerodynamic  effects  are  nullified  by  zeroing  Cc  , C , cN  , lb, 

• » i ^ A 

Ng  < rb  » Qb  * rB  at  t*ie  same  time  fixing  a , 9 , 9,  t and  cp  , the 
vehicle  can  then  be  considered  as  a point  mass.  Additional  checks  of  the 
orbital  equations  are  also  performed  by  initializing  C£  and  C^  such 
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that  the  aerodynamic  effects  of  these  equations  can  be  assessed. 

It  should  be  noted  that  dynamic  checks  for  handling  qualities  are 
normally  performed  in  five  degrees-of-freedom.  Dynamic  checks 
for  performance  (point-mass)  are  performed  in  three  and  six 
degrees-of-freedom. 

The  control  system  mechanization  is  checked  by  varying 
the  input  variables  (control  stick,  rudder  pedal,  angular  rates, 
angle  of  attack,  Mach  number,  dynamic  pressure,  etc)  and 
observing  the  control  surface  responses. 

The  cockpit  instruments,  switches,  control  stick  and  rudder 
pedals  must  be  checked  for  proper  mechanization  (magnitude  and 
sign) . 

7 . 2 Static  Checks 

The  values  being  established  for  Static  Check  can  be 
completely  arbitrary  but  they  must  be  reasonable.  For  instance, 
we  can't  choose  a test  value  for  a - -10°  because  the  Shuttle  will  not 
operate  in  that  range.  On  the  other  hand  one  can  pick  a test  value 
of  a=28°  or  a-30°  and  one  value  is  as  good  as  the  other.  In  this 
program  the  Orbiter  can  assume  either  the  ENTRY  phase  or  the 
TAEM  phase.  Each  phase  should  be  tested  individually. 

7.2.1  Equations-of-Motion  Checks:  By  setting  certain  flags 

the  FCS  can  be  effectively  bypassed  to  allow  us  to  check  only  the 
equations  of  motion.  Specifically,  to  disengage  the  FCS,  reset  the  flags 
concerning  aileron  and  rudder  and  trim  the  elevon  6e  = -16.33° 

Then  establish  the  following  initial  conditions 


Mach  - 3 

a = 20 

q = 150 

h = 99480 

V = 2970 

S = 2690 

b = 78.06 

c = 39.57 

WQ  = 155,000 

lx  = 600,000 

IY  *IZ  = 5.5*106 

rxz  = 105 

x = 6 5% 

z = 29.06% 

6BF  = ~5 

6SB  = 25° 

Gear  in  Up  position 

6e  = -16.33° 

With  the  above  IC's,  equation  (2-19)  become. 


Substituting  known  quantities  into  ( 8-2  ) 


P0=  .1667  (Rb+  PgQg)  + 52.495  [ (-.00168- . 00001  *(-16.33)  ) 8 

+ .000547  6a  +.00022  6rj  +.  68986  (-.  261  Pg  +.  073  (7-3) 

Next  we  manipulate  the  state  variables  to  obtain  a number  of 
conditions.  The  series  of  checks  appear  lengthy  at  first  glance, 
but  they  are  necessary  , because  each  set  of  conditions  allows  us 
to  verify  a different  part  of  the  set  of  equations  of  motion. 


a/ 

For 

PB 

= Q 

^B 

= 0 ; R 

B 

= .4 

rad/ sec  ; 8=  0° 

xxz 

= 0 

resulting  in 

PB  = 

. 0201  rad/ sec^ 

b / 

For 

qb 

= rb  = 

= 0 : P=  6a  : 

: «r 

= 0 ; 

PB  = 1 ' 

rxz 

= 0 

PB  = 

-.1801 

c / 

For 

PB 

= Qb 

= rb  = o ; 

6 = 

-5°  ; 

6a  = 6r  = 0 

*xz 

= 0 

PB  = 

.3990 

d / 

For 

p 

B 

qb 

= RB  = ° ; 

6a  = 

10°  ; 

8 = 6 r = 0 

*xz 

= 0 

PB  = 

.2871 

e/ 

For 

PB 

= Q 
^B 

= RB  - 0 ; 

6r  = 

; 20°  ; 

8 = 6a  = 0 

rxz 

= 0 

PB  = 

.2310 

f / 

For 

PB 

= 0 ; 

QB  = RB  = • 

4 ; 

8=  6a 

= = 0 

*xz 

= 0 

PB  = 

. 1534 

In  the  next  two  cases,  g)  and  h)  , Ixz  is  not  set  to  zero  but  to  105 
as  indicated  earlier. 

g / For  PB  = Qb  = rb  = 0 - 6a  = 10°  J 3=  6r  = 0 

RB  = -0058  Resulting  in  Pg  = .2881 

hi  For  PB  = 1 ; QB  = .4  ; RB  = 0 ; = 10°  ; 8 = 6r  = 0 

Rb  = -.3541  P = . 1148 

D B 

Still  under  the  IC's  imposed  earlier,  calculate  the  static  test  concerning 
the  pitch  acceleration.  According  to  equation  (2-19) 


A / lz  ' lx  \ P p + lxz 

(r2  - 

P2  \ 

qSc  r 

QB  ViT  J b b iy 

>vRB 

PB  ' 

+ C + AC  + 

Iy  L m0  me 

c + c 

m6BFl  m£SBl 

(6sb-  550) 

]+  2VIy(CraQQB) 

(7-4) 

Plugging  in  the  inertia  terms 

Qb  = .8909PBRg  + .0182  (r|  - p| ) + 2.903  . 023  + ACme  - .00051  5gF 

+ . 00023  ( 6sr  - 55°  )J  + .0193  (-2.3  Qg  ) (7-5) 

As  a consequence 

i/  For  Pg  = Qg  = Rg  = 0 ; 6e  = 0 ; 6gF  = -5°  ; <^g  = 25° 

Whereupon  the  check  value  for  the  pitch  acceleration  is 

Qg  = -.0838  rad/sec^ 

H For  PQ  - Qg  = Rg  = 0 ; 6e  = -16.33°  ; 6gp  = 16.3° 

6sb  = 25°  Qr  = -.0831 

k/  For  PB  = Qb  = Rg  = 0 ; 6e  = -16.33°  ; 6gp  = -5° 

%B  = 87.2°  Qg  = .072 
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For  PB  = 1.  ; QB  = HB  = 0 ; 0 = 6a  = 5r  = 0 
JXZ  = ° Rb  = -0008 


q / For  PB  ” = Rb  = 0 ' B = -5°  ; = 6r  = 0 

XXZ  = 0 Rb  = .0275 


■JL  -w*-  . - 
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r / 


s / 

t/ 


B 


For  P 

*xz 

For  P 


B 


I 


XZ 
For  P 


B 


I 


XZ 


= qb 

= 0 

= 0 
= 1. 

= 0 


= rb 


'•  qb 


= o ; 6a  = 20°  ; 6r  = 0=  0 

R0  = .0011 

= 0 ; 6r  = 20°  ; 0 = 6a  = 0 
R_  = -.0481 

D 

= .4  ; RB  = 0 ; 0 = 6a  = 6r  = 0 
R„  = -.3556 

±5 


Now  we  again  impose  a value  of  105  on  Ixz 

u/  For  pb  = ^b  = rB  = 0 5 6a  = 20°  J 3=^  = 0 

Pg  = .5761  rad/sec2  RB  = .0116 

v/  For  PB  = 0 ; Qg  = Rg  = . 4 ; 6a  = 20  ; 0 = 6r  = 0 

pB  = .5913  Rb  = -.0184 


7.2.2  FCS  Checks  : Set  ALTITUDE  = 92,800  ft  ; VELTRUE  = 2960  ft/sec 
1 URNCORD  = .TRUE.  ( to  allow  the  FCS  to  maintain  the  nose  level 
during  a turn)  and  Rudder  Pedal  Position  RPTASOP  = 0.  Then  proceed 
to  parts  a)  and  b)  below. 


a)  Pitch  Plane  Test  (ENTRY  Phase)  : Fix  the  Pitch  RHC  , 
PRHCSOP  = -10°  and  QB  = 2°/sec.  Then  the  eleven  deflection,  the 
final  object  of  our  Pitch-plane  test , should  be  6e  = -6.42°.  The  cal- 
culation of  intermediate  variables  leading  to  &e  is  left  to  the  reader. 

b)  Roll  and  Yaw  Test  (ENTRY  Phase)  : Establish  the  Roll  Stick 
Position  RRHCSOP  = 16.15°  and  also  the  following  conditions 


PRHCSOP = 

0° 

Pg  = 5°/sec 

CD  = 

CO 

o 

o 

Rg  = 1°/ sec 

0 = 

o 

o 

a = 15° 
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We  should  obtain  these  final  results 

6 = - 2.66° 

a 

6r  = -1.9° 

and  Yaw-jet  Command  DRJET  = .63  (dimensionless)  corresponding 
to  the  firing  of  two  reaction  jets  (UZCMD  = 2). 

Now  we  are  ready  to  test  the  FCS  in  TAEM  phase.  Assume 
ALTITUDE  = 48,600  ft  and  VELTRUE  = 1017  ft/sec . This  will  result  in 
Mach  =1.05  and  q = 200  lb /ft*2  . Proceed  to  parts  c)  and  d)  next. 

c)  Pitch-Plane  Test  (TAEM  Phase)  : Choose  PRHCSOP  = -10°  ; 
RRHCSOP  = 0°  ; RPTASOP  = 0°  ; QB  = -2°/sec  ; Rg  = -2°/sec  and 

cp  = 30°  . With  the  above  initial  conditions  we  should  command  the 
following  responses: 

Pitch -jet  command,  DPJET  = -.79 
Pitch  command  signal,  BCSL  = .365 

d)  Roll  & Yaw  Test  (TAEM  Phase)  : Select  the  conditions  listed 
below. 


PRHCSOP 

= 0° 

RRHCSOP  = 11.15' 

RPTASOP 

= 11.125° 

Pg  = 2°/sec 

cp 

= 30° 

R = 5°/sec 

e 

= 20° 

NY  = . 25  g 

a 

* 15° 

rb  = o 

Look  for  the  following  results  out  of  the  FCS 

6 = .771° 

a 

6r=  22.09° 

DRJET  = -10.68  (4  jets) 
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7 . 3 Dynamic  Checks 

Again,  two  principal  parts  of  the  program  are  tested  semi- 
independently  of  each  other.  We  choose  to  test  the  simulation  in 
5DOF's  only,  by  setting  FIVDEG  = . TRUE.  . In  5~DOF  mode  the 
Shuttle  flies  at  a constant  altitude  at  a constant  axial  velocity. 

To  disturb  the  equilibrium  (or  steady-state)  conditions,  we  introduce 
a forcing-function  of  short  duration  and  observe  the  system  response. 
The  triangular  -waveform  forcing-function  shown  in  paragraph  7.1  is 
generated  on  the  analog  computer  and  patched  into  an  A-D  convertor. 
It  has  its  peaks  at  ±1  . Once  converted  into  a digital  quantity,  this 
waveform  is  further  modulated  to  give  the  proper  pulse  amplitude 
(for  example  +16°  for  6e  pulse).  During  each  test  only  one  surface 
is  disturbed.  All  other  surface  deflections  are  returned  to  initial- 
condition  values  or  set  to  zero.  Both  the  equations -of- motion  test 
and  the  FCS  test  use  the  same  standard  conditions  for  altitude, 
true  airspeed,  weight,  surface  area,  inertial  coefficients,  Euler 
angles  and  trim  surface  deflections,  namely 


Ix  = 6*105 

x = . 65 

VT  = 2970 

Iy  = 5.5*106 

y = o 

q = 150 

Iz  = 5.5*106 

z = .2907 

Mach  = 3 

IXY  = 0 

h = 99,480 

S = 2690 

>xz  ■ 105 

WQ  = 1.55*10 

5 b = 78.03 

c = 39.56 

= 107.5 

The  following 

angular  conditions  are 

also  arbitrarily  imposed 

^5B  = 25° 

cp(  0 ) =0° 

L (0)  = 0° 

6BF  = '5° 

Y(0)  = 0° 

\(0)  = 0° 

a (0)  = 20° 

6(0)  = 20° 

6=0° 

a 

3(0)  = 0° 

6e  = -16.33° 

6r  = 0° 

7.3.1  Equations -of-Motion  Tests  : Three  runs  are  made  , all  with 
the  following  conditions: 

FCSDE  = FCSDA  = FCSDR  = .TRUE. 

Various  output  quantities  are  fed  to  the  DAC's  and  appear  on  the  strip 
charts  as  analog  signals.  Under  the  best  circumstances  the  traces  can 
be  read  at  an  accuracy  of  1%  and  thus  can  be  used  as  a qualitative 
check  only.  Distinguish  the  three  runs: 

- Case  a : 6e  is  pulsed  with  a triangular  wave  (see  page  65) 

having  ±16°  peaks  . See  figure  7-1. 

- Case  b:  6a  is  pulsed  with  the  same  waveform  resulting  in 

the  traces  on  figure  7-2. 

- Case  c:  6r  is  pulsed  with  a triangular  signal  having  +20° 

peaks.  Consult  figure  7-3  for  the  airframe  responses. 

7.3.2  Flight-Control -System  Tests:  The  FCS  tests  helps  to  verify  the 

end-to-end  performance  of  the  simulation.  Since  we  operate  at  Mach  3 
and  since  rudder  pedal  inputs  are  bypassed  by  the  FCS  at  Mach  >1.5, 
there  is  no  need  to  pulse  the  rudder  surface.  In  the  two  remaining 
test  cases,  a full  stick  input  is  applied  to  the  aileron  in  one  case  and 
to  the  elevon  in  the  other  case.  The  testing  is  intended  not  only  for 
dynamic  system  response,  but  also  in  assisting  in  the  discovery  of 
possible  FCS  discontinuities  or  errors.  Both  cases  are  made  with  the 
following  conditions: 

FCSDE  = FCSDA  = FCSDR  = .FALSE. 

As  before,  the  system  response  is  observed  through  the  use  of  strip- 
chart  recorders.  The  two  cases  are  perturbated  in  the  following  manner: 

- Case  d : is  pulsed  with  a triangular  wave  having  peaks  at 

±20°  . See  figure  7-4. 

- Case  e : 6a  is  pulsed  by  the  same  signal.  See  figure  7-5. 

RHC 

It  should  be  noted  that  in  both  cases  the  Flight  Control  System  must 
be  allowed  to  stabilize  the  airframe  for  approximately  3 to  5 seconds 
before  applying  the  disturbances.  In  this  manner  the  pitch  axis  will 
settle  to  a trim  condition,  reducing  the  inter-axis  interaction. 
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8.  TERMINOLOGY,  SYMBOLS  AND  DEFINITIONS 
8 . 1 List  of  Acronyms 


ASCII 

A-D 

CG 

American  Standard  Code  for  Information  Interchange 
Analog -to-Digital 

Center  of  Gravity 

CMD 

CRT 

C/R 

Command  ; Commanded 

Cathode-Ray  Tube  (Computer  terminal) 

Carriage  Return 

DAC 

DOF 

EOF 

Digital -to-Analog  Convertor 

Degree  Of  Freedom 

End-Of-File 

FCS 

FGP 

FILT 

Flight  Control  System 

Function-Generation  Package 

Filter 

Hex 

IC 

I/O 

Hexadecimal  number 

Initial  Condition  mode 

Input  / Output 

NOM 

RHC 

RCS 

Nominal  value 

Rotational  Hand  Control 

Reaction-jet  Control  System 

SEL 

Systems  Engineering  Laboratories 

\ 


8 . 2 Lists  of  Symbols  Used  in  Computer  Programs 

The  following  abbreviation  is  used  throughout  this  chapter:  DL  = dimen- 
sionless quantity.  The  subscripts  on  various  quantities  have  the 
following  significance: 

6 a = aileron  incremental  deflection 
B = written  in  Body  coordinate  system 
6BF  = Body  Flap  incremental  deflection 


l 
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6e  = elevon  incremental  deflection 

6r  = rudder  " " 

6SB  = Speed  Brake  incremental  deflection 

E or  e = Earth  axes  or  relating  to  earth 
F = Fuel 

G = Gear  extended 

H = written  in  H-frame  (orbital  reference  axes) 

L , R two  halves  ( Left  & Right)  of  the  elevon 
0 = bias  value  ; IC  value  ; sea-level  value 


Be  aware  that  stability  coefficients  with  various  final  subscripts  may 
have  different  units.  Consider  , C™  , ( bias  value  of 

nornal  force  coefficient  , normal  force  derivatives  due  to  elevon  and 

body  flap  incremental  deflection)  ; they  will  not  appear  as  three  separate 

entries  on  the  table . Note  that  Ctvt  is  dimensionless  while  CN  and 

1N0  6e 

N^p.v  have  units  of  sec-1  . Other  subscripts  used  are  a,  p,  p,  q and  r. 

Sometimes  a quantity  is  represented  by  two  computer  variables  which 
are  not  equivalent.  For  instance  ALPHA  stands  for  a in  radians  . At  the 
same  time  ALFA  also  symbolizes  a but  in  degrees.  Similarly  LAT  and 
LA  both  stand  for  latitude,  but  LAT  is  a single-precision  quantity  while 
LA  is  computed  in  double  precision.  The  computer  symbols  appear  on 
two  separate  tables.  The  second  table  lists  the  names  used  in  the  FCS 
program  and  the  first  table  covers  the  remaining  programs  ( MAIN 
plus  all  subroutines) 


Table  8-1  : MAIN  PROGRAM  SYMBOLS 


Symbol 

» 

Explanation 

Computer 

Representation 

Unit 

a 

b 

Sonic  speed 

Wing  span 

SOUND 

B 

ft/  sec 
ft 

cc 

Total  trim  chord  coefficient 

CC 

DL 

cgx 

Longitudinal  distance  of  CG  measured 
from  a reference  point  situated  238" 
in  front  of  the  nose  (expressed  in  % 
of  Orbiter's  lenght) 

(not  used) 

% 

Total  rolling- moment  derivative 

CL 

DL 

c 

Cm6e 

Mean  aerodynamic  chord 
Pitching-moment  derivative  due  to  6g 

C 

CMDE 

ft 

deg"1 

Cm6SB 

Pitching-moment  derivative  due  to 

CMSB 

deg'1 

Cn5r 

Yawing-moment  coeff.  due  to  rudder 

CNDR 

deg  1 

CN 

Total  normal  force  derivative 

CN 

DL 

CVB 

Side-force  coefficient  due  to  8 

CYB 

, -1 
deg 

ft 

Fuel  consumption  rate 

(not  used) 

1 / sec 

g 

h 

Gravitational  acceleration 

Geometric  altitude 

GRAVITY 
or  G 

ALTITUDE 

ft/ sec^ 

ft 

% 

Geopotential  altitude  (used  for  calcu- 
lating air  density  & temperature) 

GEOALT 

ft 

lxx 

Polar  moment  of  inertia  (along  the 
longitudinal  axis).  Or 

IXX 

slug,  ft^ 

lXY  ' 

rxz 

► 

Products  of  inertia 

IXY  ^ 

IXZ 

y 

slug,  ft" 

lYZ 

/ 

IYZ 

1 

Table  8-1  : Continued 


Symbol 

Explanation 

Computer 

Representation 

f 1 

Unit 

*YY;IZZ 

Transverse  moments  of  inertia 

IYY  ; IZZ 

slug,  ft ^ 

KCAS 

Calibrated  airspeed 

KCAS 

knot 

KEAS 

Equivalent  airspeed 

KEAS 

knot 

J2 

A constant  characteristic  of  the 

earth  mass  distribution 

J2 

DL 

Jp 

Constants  characteristic  of  the 

r UYCMD 

I 

j 

J"R 

thrusts  of  the  pitch,  roll  and 

- 

UXCMD 

DL 

J Y 

yaw  reaction  jets. 

v UZCMD 

L 

Latitude 

LAT  OR  LA 

radian 

lb 

Total  rolling-moment  in 

body  axes 

LBODY 

lb.  ft 

h ■ h-h 

Direction  cosines 

LI  ; L2  ; L3 

DL 

Length  of  the  Orbiter  (nose  to 

tail) 

LENGTH 

ft 

m 

Orbiter' s mass 

MASS 

slug 

mi  ;m2:m3 

Direction  cosines 

Ml  ; M2  ; M3 

DL 

! mb 

Total  pitching  moment 

MBODY 

lb.  ft 

Mach 

Mach  number 

MACH 

DL 

nl  ' n2  ’ 

n3 

Direction  cosines 

N1;N2 ;N3 

DL 

Total  rolling  moment 

NBODY 

lb.  ft 

PB 

Rolling  rate 

PBODY 

rad /sec 

q 

Dynamic  pressure 

QBAR 

lb/ ft 2 

qb 

Pitching  rate 

QBODY 

rad/sec 

Ji 


J.  - 


; 


Table  8 -1  : Continued 


Symbol 

Explanation 

Computer 

Representation 

Unit 

ro 

Equatorial  radius  of  earth 

RO 

ft 

r 

6R 

Distance  between  Ortiter  and 
earth  center 

Altitude  measured  from 
sea  level 

radius 

or 

RAD 

DELR 

ft 

ft 

rb 

Yawing  rate 

RBODY 

rad/ sec 

l 

rd 

Down  range 

DOWNRNG 

ft 

1 

Rx 

Cross  range 

CROSSRNG 

ft  j 

s 

Reference  wing  surface  area 

S 

ft2 

T 

Thrust 

(not  used) 

lb 

TX'TY 

TZ 

Thrust  vector  in  the 

3 body  axes 

(not  used) 

lb 

UEA*  V 

UE  ' 

EA 

Horizontal  components  of 
airspeed  in  earth  axes 

Orthogonal  velocity  components 

UAIR(or  UA) 

V AIR  (or  VA) 

f UEARTH 

ft/sec 

ft /sec 

VE 

WE 

/ 

> 

expressed  in  earth  axes 

< VEARTH 
l WEARTH 

VELTRUE 

V 

Orbiter's  total  velocity 

or 

TAS 

ft  / sec 

w0 

Orbiter's  weight 

WEIGHT 

lb 

WF 

Fuel  weight 

(not  used) 

lb 

Wu 

Wind  component  in 

North  direction 

NW1ND 

ft/sec 

Wv 


Wind  component  in 
East  direction 


EWIND 


ft/sec 


Table  8-1  : Continued 


Symbol 


Explanation 


Computer 

Representation 


Longitudinal  location  of 
Orbiter's  CG  (expressed  in 
percent  of  Orbiter's  length) 


XBAR 


Longitudinal  acceleration 
component 


XBODY , XB 
or 
AX 


Lateral  location  of  CG 


Lateral  acceleration  component 


YBAR 

YBODY , YB 
or 
AY 


Vertical  location  of  the  CG 

Normal  acceleration 
component 


ZBAR 

ZBODY,  ZB 
or 
AZ 


Angle  of  attack 
Angle  of  side  slip 

Flight  path  angle 


ALPHA 

BETA 

(not  used) 


Aileron  total  deflection 


DAIL  or  DA 


Body-flap  total  deflection 


FLAP  or  DBF 


Elevon  total  deflection 


ELEV  or  DE 


Rudder  total  deflection 


DRUDor  DR 


Speed  brake  total  deflection 
Time  increment  (time  frame) 


BRAK  or  DSB 


DELTA 


Pitch  angle 

Angle  between  surface  of  the 
fuel  and  Xg  axis 


THETA 


(not  used) 


Unit 


. L, 


Table  8-1  : Continued 


F 

Symbol 

Explanation 

Computer 

Representaion 

Unit 

X 

Longitude 

LNG  or  LN 

rad 

P 

Colatitude  ( ) 

MU 

rad 

P 

Atmospheric  density 

RHO 

slug/ft^ 

cp 

Roll  angle 

PHI 

rad 

r 

Yaw  angle 

PSI 

rad 

We 

Angular  velocity  of  earth 

OMEGA 



rad/sec 

Some  variables  appear  both  in  the  Main  program  and  the  SHTFCS. 
If  they  have  been  previously  listed  on  Table  8-1,  they  will  not  be 
repeated  here  on  Table  8-2. 


Table  8-2  : FCS  Variables 

Computer 

Variable 

Description 

Unit 

DACMD 

Aileron  command 

deg 

DATMPAN 

Aileron  panel  trim  command 

0,-1 

DATMRHC 

Aileron  stick  trim  command 

O.  + l 

DATSUMI 

Aileron  trim  integrator  command 

deg 

DBFMAN 

Manual  Body  Flap  command 

O.  + l 

DBFRC 

Body  Flap  rate  command 

deg /sec 

Table  8 -2  : Continued 


Computer 
V ariable 


DECMD 


DETRIM 


DETMPAN 


DETMRHC 


DPJET 


DRCPF 


DRCMD 


DRJET 


DRLP 


DRRP 


DRTRIM 


DSBPC 


ELFBK 


ETRIM 


Description 


Aileron  command  due  to  stick  and  rate  feedback 


Elevator  command 


Left  elevon  deflection 


Right  elevon  deflection 
Elevator  position  trim  command 
Elevator  panel  trim  command 


Elevator  stick  trim  command 


Pitch  Jet  Command 


Elevator  Command  due  to  stick  and  rate  feedback 


Rudder  command  due  to  pedal 


Rudder  command 


Yaw  Jet  Command 


Left-half  rudder  deflection 


Right-half  rudder  deflection 
Yaw  Trim  Integrator  output 
Speedbrake  command 
Elevator  Feedback  signal  ( Mach^l2  ) 
Elevator  Position  Trim  ( Mach>12  ) 


Lateral  Acceleration 


Normal  Acceleration 


! 


Table  8 -2  : Continued 


Computer 

Variable 

Description 

Unit 

PE 

Roll  & Yaw  stability  feedback 

DL 

PRHCSOP 

Pitch  Stick  Command 

deg 

QDOT 

Pitch  Acceleration 

rad/ sec^ 

RDOT 

Yaw  Acceleration 

rad / sec^ 

RPTASOP 

Rudder  Pedal  Command 

deg 

RRHCSOP 

Roll  Stick  Command 

deg 

SBHP 

Speedbrake  Handle  Command 

deg 

SUM23 

Roll  & Yaw  Command  due  to  stick 

DL 

1 

_ . i 
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PKOGWAM  SHILOFT 


SPACE  SHUTTLE  OKHiUL  FLlGHf  TEST  SlMULATlUN. 

IMPLICIT  H t A l ***(A  - l) 

REAL*8  DATAQBL 

logicalm  0 A I A H Y I 

COMMON  /ArtKAVl/i)ATA*HDM(400) 

COMMON  /AKWAY2/DAI AOBLL17) 

COMMUN  /AWhAVi/OATABYT (17) 

COMMON  / AHWA  Y**/ANGLE  ( 3 ) 

BIT  iNBIT(lb) 

LUG  1 CAL*  1 I C,  HULL),  POISE  I , A HOW  I , I HEOEG,F  I VOEG 
LOGICALM  FCSOA,FCSDE,FCSOK,FCSDF ,FCSSB,*1ND 
LUUICALM  KEAOY,L)UMPFLAG,  AUTOBF 
LOGICALM  TOHNCOKO 

LOGICALM  UtPULSE,DAPULSE,OWPULSE 
LOGICAL*1*  tASEH 
logical***  oiscwriN 

I N 1 EGEW*2  IEMPI , TEMP2,PWEVTEMP 
IN1EGEK*4  SAVFL 
IMEGEK*«  LINK, II, JJ 
INTEGEW*4  handle* 

WEAL *8  DP  IT  1, OP  I l 2, DPI  1 3 , CUSL A T , COSMU, COSPS I H, L A , LN, MU , * AO, 

* WCUSLAT,SHAHE,SINLA| , SINMU, S1NPSIH,UA, UH,UE, VA, VB, VE , Kb, WE , 
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* *HP,XB,  Y B, Zb, PSlM,OELH,PELR PUT, UHAV IT  Y , L A I Du T , LNGOU 1 , PSIHOUT, 

* HUH,  HUHDU  l,UH,*H,XEAH|H,XH,YEARIH,YM,ZEAHTH,ZH 
HEALM  UUTDAT A(S243) 

INTfcGtRM  AEHOS1ZE 

THE  FULLOwJNC.  TYPING  STATEMENTS  UImENSIUN  T Mt  DATA  TAbLt  ARRAYS. 

INTEGERS  DHL  S I ZE , HS I ZE , (SIZE »hAFSIZE*BYT  SIZE, SIZE, PLACE 
INTEGEH*2  ELEVPI (2) ,PELPT (2) ,UERPl ( 2 ) , ALPhAP T ( 2 ) , maCHP T ( 2 ) , 

* BETAPT (2 ) ,USdPT (2) 

INTEGEHM  LIST (434) 

THE  F OLLPw 1 NG  AHKAYS  FURM  The  1/U  TABLES, 

UATAUUT  ANP  DAIAIix  form  the  44  Channel  BATA  UUTPUT  AMP  32  CHANNEL 
PATA  INPUl  TABLtS,  RESPEL l 1 VtL Y , 

SCALUUT  ANP  5 C A L i N AxE  THE  C URRE SPUND I NG  OUIPUT  ANP  INPUT  SCALE 
FACTORS,  KESPEC I IVELT, 

IN  I EGtK*4  CUUNT 

I N I hl»ER*2  UAFAUUT  ( 44 ) , DA T A I N ( 32 ) 

REALM  SCALUUT (44) ,SCALIN(32) 


REALM 

ccur (63) 

REALM 

PCCE  T (37B) 

REALM 

CCPHF I T (63) 

HE AL*4 

CLPBF2T (63) 

REALM 

CCPSBl T (63) 

HE AL*4 

CCDSB2T (63) 

REALM 

CLb  T ABL ( 252 ) 

REALM 

CLHPE  1 T (63) 

REALM 

CLBUE2T (63) 

REALM 

CLBPSB1 I (63) 

REAL*4 

CLB0Sb2T (63) 

RE  AL  * 4 

CLBOSH3T (63) 

RE ALM 

CLP A T ( 3 16) 

HE  AL  * 4 

CLPRT (252) 

RE AL*4 

LLPT (63) 

HE AL*4 

CLRT (63) 

HE ALM 

L Hp  T (63) 

H E A L * 4 

DCPE  T (376) 

REALM 

CHDBF IT (63) 

HEAL *4 

CMUBF2I (63) 

REALM 

CMPSH1 f ( e 3 ) 

REALM 

CMOSB2  T (63 ) 

REALM 

CMuT(2b2) 

RE  AL  * 4 

CNb  T ABL ( 252 ) 

HE  AL  * 4 

CNBPE l l ( 63) 

REALM 

CMBPE2T (63) 

HE AL*4 

CNBUSb 1 1 (63) 

HE  AL  * 4 

CNBDSB2  f ( 63 ) 

HE  AL  * 4 

CNBUSB3  T ( 63 ) 

REALM 

CnDA T (378) 

HE  AL  * 4 

CNORT (252) 

RE AL*4 

CNPRBT (63) 

realm 

CnP  T (63 ) 

RE  ALM 

CNR  T ( 63 ) 

HE  AL  * 4 

CN0T(63) 

RE  ALM 

CNPET (63) 

HE AL*4 

CNPBF I T (63) 

REALA4  C YU  A T (378) 

REAL*4  C YOR  T (63) 


REALA4  C YURSo 1 T ( 7 ) 

WEALA4  CYURSB2I  (7) 

REAL  *4  SND  I AdL (31 ) . SNUAWi,(3D , ZE  T ABL (41) , ZE  T ARG(41 ) 

HE  AC *4  MNUARG(20)  , Nw  r ABL (20)  ,EMA«L(2v) 

REAL *4  ELEv  1 (0) ,OELT (2) ,0ERT(2) p ALPHA! ( 1 1 ) , M ACH  T ( 9 ) , BE T A T ( b ) , 
* DSBI(fe) 


EQUIVALENCE  (INoll pOISCKI In) 

EQUIVALENCE  (SnDARGU  ) , ZETARG(  1)  ) 

EQUIVALENCE  (ELEV  I (2)  pELfcvPT ( 1 ) ) , (UtLT  (2)  pDELPT  ( 1 ) ) p (L>ERT(2)  , 
aUEKPT ( 1 ) ) , (ALPHA  T (2) , ALPHAPT ( 1 ) ) , ( MACH T ( 2 ) p M ACHP T ( 1 ) ) , ( BE T A T ( 2 ) 
a p BE  T APT 1 1 ) ) p (USBK2)pOSBPT(l)) 

EQUIVALENCE 

* (UUTDAT A ( 1 ) ,CCU r ( m 
a, (0UTUATAC64) ,0CCE  H1D 

a,  (OUIUA  I A (442)  pCCOBUT  (1  ) ) 

*,(UUTUAlA(b05),CCUBF2I(l)) 
a, (OUTOAT  A ( 560 ) pCCUSBl T ( 1 ) ) 
a,  (UU  rDAT  A (e>31  ) .CCUSB2  f ( 1 ) ) 
a . (UUTUATA(694) p C L B I ABL ( 1 ) ) 
a,  (UU I DATA (946) ,CLBOEl 1 ( 1 ) ) 

*,(uuiuArA(iu09),cLoi)E2r(i)) 

* p ( OU I U A I A ( 1 0 72) pCLBuSBl  I ( 1 ) ) 

*,  (UUTDAT A ( U 3b)  , CLBU5B21  ( 1 ) ) 
a,  (OUTUA f A ( 1 198) . CLBUSB  3T ( 1 ) ) 

*,  (OUIi)A!A(  126  1 ) » CLUA  T (1)  ) 

A, (OUTUA T A( 16  39) , CLURT ( 1 ) ) 
a,  (UUI DAT  A ( 189  1 ) , CLP  ( ( 1 ) ) 

a,  (0uTUATA(1954),CLHT  (D) 
a, (UUTUATA(2017) ,CMUT  (1)) 

*, (OUTUA T A (2080  ) , DC MET ( 1 ) ) 
a, (OU I DA  I A (245B) pCMOBE 1 T ( 1 ) ) 
a,  (UUTUAT  A (2521  )pCmi)bE2T(1  )) 

A, (uul DA f A (2504) pCMUSdl T ( 1 ) ) 

*, (UUTOATA(2H4  7),CMUSB2T(l  ) ) 
a,  (OUIUAT A (271 0) , CMQT ( 1 ) ) 
a, (UUT0ATA(2962) ,CN0TABL(  1)  ) 
a,  (UUTUA  ( A ( 3214)  pCNBOEl  HD) 
a,  (OUTUA  T A ( 32  7 7 ) « CNB0E2  1(D) 
a,  (UU1  uA  1 A ( 3 340)  pCnbUSBI  HD) 

*p (UUTUAT A (3403) ,CNB0SB2T (1)  ) 

Ap (UUfUA I A( 34  66 ) pCNBUSb3I ( 1 ) ) 
a,  (UUTUAT A<3529) ,CnDA  HD) 

*p (UUTUA  T A(3907) pCNUR  H 1)  ) 
a , (UUTUA1 A (4159) p CNOkh f ( 1)  ) 

A,  (UU1  UA  r A (4222)  pCnP  HD) 
a,  ( UUTUA  [ A (4205)  .CNR  HI  ) ) 
a, (UUI DATA (434B) pCNUTl  1)  ) 
a,  (UUTi.A  I A (44  l 1 ) pCNUET  ( 1)  ) 
a,  (UUTUA  I A ( 4 4 7 4 ) , CNUHE  11(D) 
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* $ l UU  1 0 A | A l«b37)  , CN0BF21  (U) 
*,  (CIUTOAT A(afeOO) » CnUSB 1 1(1)) 
*,  (UU1 OAT  A (<466  j)  , CNOSH2  I ( 1 ) ) 
*, (0UTDATAU726) ,CY«T ( 1 ) ) 

*,  (UUlOATA(**7  89),CY0AT ( 1 ) ) 

*,  (OUTDAT A(5167) ,CYORT ( 1 ) ) 

*,  (OUlOATA(S230),CYOWSB11  (1)) 
*,  {OUTUA?A(5237),CYORSB2T(l)i 
EUUl V ALt^CE 

* (Da  I AOBL  ( 1 ) > DELR ) 

* , (OATAOBL(2),OELROOI ) 

* , (DaTAObL  (3)  ,liWAVl  TY) 

* , (UA  TAOBL(R)  ,LA1U0D 

* , (DATAOBL(B) , LOGOUT) 

THERE  13  A HOLE  AT  U A T AOBL  ( 6 ) 

* , (DATAOBL(7),PSlhD0T) 

* , (OATAOBL C » ) , HUH) 

* , (OA  I AOBL  (9)  t RUHOO  T ) 

* , (DA f AOBL ( 1 0) ,UH) 

* , (OAT AOBLl 11 ) , AM) 

* , (oatadbl(1£),xearth) 

* t (OATAOBL! 13) t xh) 

* , (OA 1 ADBL ( 1«) , YE ARTH) 

* t (OATAOBL(lb) < YH) 

* / (OATADBL( lb) >ZEAkTH) 

* , (OA  T AOBL ( 1 7 ) , ZH) 

*,  (OATAwRuRU  ) ,C  ALPHA 0) 

* , (UA  I AwROR (2)  , Ax  ) 

*, (DAlAwROH(3)» AT) 

*, (OA1 AwROR (A) , AZ) 
*,(OATAwHOR(S),CALRHA) 

*, (OA (ARROW (to) ,ALl 1 IUOE) 

* , (OAT  A wROH ( 7 ) , H ) 

* t (0ATAwR0R(.6)*ChETAQ) 

*,  (OAJAhkOR(9),BE  IA) 


* , ( 0 A 1 AwHOR  (10) 

.DRAM 

* , (OAT AhROR (11) 

,C) 

• , ( OA  T AAKUR ( 12) 

»CC  ) 

*, (OAT AwROR (13) 

,CCBF) 

* , ( OA  T AWHOK  ( 1 <4 ) 

,CCSB) 

* # (OAT  AwHOR ( 1 b) 

,CO) 

*,  ( U A 1 ArtROR ( 16) 

*LL) 

*, (OATAwHON(  1 7) 

, CLB  T ) 

* , ( OA  r AWHOK ( 1 «) 

,CM) 

* , (OAT  A wRDK ( 19) 

,CMBF) 

* , (OA  1 AwRUR  ( 20  ) 

,CMSB) 

*, ( 0 A I ARROW (21) 

,CN) 

*, (OAT  AWROR ( 22  ) 

» CObF ) 

*, (OATAwHD R (23) 

» CNSB ) 

* , (OAT  AWROH  ( 2<i  ) 

,CNB  ( ) 

*, (0ATAhRDR(Z5) 

,CY) 

*,  (DATAwWOH(2to) 

t CLBOt  ) 

*, ( 0 A 1 A wroR (27) 

* CLoSR ) 

*,  (OA  I AWROR  (2«) 

» CNnOE ) 

*, (OATAwROR(29) 

, RUNTIME) 

* , (OATAWRLIK(iO) 

* TIME) 

* , (OA  r AWRUR ( 31  ) 

« CNBSt)  ) 

*,  (OA)  AwR0R(32) 

»CYORbtj) 

*, ( DA T AwkDk ( 33) /CKDSSkNU) 
*,  (DAT  A WkDk  ( 34  ) ,DAIL) 

*, (DAT AwkDk  C 35 ) »DfcU 
*, (DAT  AARDk 136) * AL  T 0 ) 

*, (DA  TAwkDR(37) /PkHCSUP) 
*,  (DArAkRQK(38), Dfck) 

*, (DAT AkROK ( iV) /DHIC) 

*,  (DA  ( AwkDK(40)  ,DkjO) 

*, (DA  I AWWDK ( 41 ) ,DBFIC) 

*,  (OATAWKDK(4«>)  ,DEHC) 

*, (DAT AkKDk (43) * D A C M ) 

*, (DAI AkkUK(4«) /DfcLHIC) 

*, (DATArtkDk(45)/DfcLIA) 

*, (DA  I a wkDk (46) z Ofck  1 C ) 

* z (DAfAkkDk(47)  , DUWNRNG) 

* t (DA  I AkK0R(48) » USB  1 C ) 

*, (DAf AkRDk(49),tLEV) 

*,  (DA  I AkWDk(30) zfcklkD) 

*, (DAI AwkDk  (31  ) , FLAP ) 

*,  (DATawRD«(32) /GO) 

*, (DA TAkHDN (33) /GEAR) 

*, (UAI  AWKDM34)  zOPJET) 

* z (dAT  a k w Dk  ( 33 ) z DK  J E T ) 

*, (DAI AkkDk (3b) /MDOT ) 

*, (yA)  AkKDK(37),lXI“1Y) 

* , (DA  I A«RDK(38) ,GtUAL  f) 

*, (DA  I AkkQH (39) , IXX) 

*,  (DATA«HDH(60)  , 1XY) 

*,  (DA  I AkkD*  ( t>\  ) / 1 XZ) 

* , (DAT  AkkDk  ( 6<?  ) / 1 YMZ  ) 

*, (DA  I AwkOk(ft3) / I YUY) 

*,  (DAI  AWkDK(t>4)  / 1 Y Y ) 

*, (DA  I AkkDk (63) / 1 ZMX) 

* z (DATAkkUk  (66)  / I ZUZ  ) 

* / (DA  I AkWOK (67 ) / I ZZ ) 

*/ (DAI AkkDK(b8) / I P UNC  ) 

*,  (DA  I ArtkDk  l 69)  , J<?) 

* , (DATAkkDk(70)/KCAS) 

* , (DA  1 AkkDk (71) / KfcAS) 

(DA  T AkROk  (7^) /kkhCSuP  ) 
*, (DATAWkDH(73)/*N0fCUN) 
*, (DAT AkkDk ( 74)  , kPTASOP) 
*,  (DA  I AwkDK  ( 73)  /CLAD 
* , (DA T AkRDK ( 76) zCLAJO) 

* , (DA  I AkkDk (7  7) / C L k G ) 

* , (DAT  AkROk  ( 78 ) ,CLNGu) 

*, (DAfAkkDk(7V),UZCMD) 

*, (DA  I AkkOk(ttO) zUXCMO) 

*z (DAT AkkDk (61 ) /LBDDY) 

* t (DA  ( AkkDk(8<?)  / LtNGfH) 

*, (DAI AkkDk (65) zUYCMU) 
*,(DAIArtkDK(64j,MACH) 

* / (DA  T AwkDK  (83)  / ^ASJi) 
*/(DAIA«vKDk(66)/CMD) 

* z (DAT  AkkD«(67) > MbODK) 

* , (OAT AwRDK (88) , NBUDY ) 

*, (DA  I AkkDk (69) , NklND) 

*, (OAT AkkDk (90) / D6L  A f t ) 

* / (DAI AkkDk (91)/ (JMfcGA  ) 
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*,(0AlA*H0K(9£),Rt)lC) 

*,  (DATAWN0K(93)  ,ROUT) 

*,  (OAT  ArtRUR(9<4) 

*, (OAI AwR0H(95) ,CPH10) 

*,  (0ATAWK0K(9t>),CPSI  J 
*,  (DATA*R|)R(97)  # CHS  I 0 ) 
*,(0ATA*HDR(98)#CPSlH) 

*, 10AT ArtK0R(99) #CHSlK) 

*,  (OATAAkORdOOdRHOOY) 

*, (OA I AwHDH( 1 01 ) , PH 100  I ) 

*,  (OAT  AmROR  ( 1 Q«J  ) , Pfc  ) 

*, (0AIAwRUK(103)#UBAHJ 

* , (OA  T A AH OH ( 1 04 ) , WBIC  ) 

*, ( OA  f A wKOR (105)#  000  T ) 

*, (OAT A*KOK ( 1 Ob) # UHUO  Y ) 

*, (OAT  ArtkDH ( 10/)#K0) 

*, (oatawhoh( io«)#kbiu 

* # (OAT AaKOH ( lu9 ) # ROO  r ) 

* # (OA 1 AnkOR ( 1 1 0 ) , KHO) 

*,  (OAT  ArtROR ( 1 1 1 ),RhUO) 

* , (OA  I aakOR  ( 1 t<2 ) , SWHP) 

*, (OATAwHOK ( 1 1 3) , RAO l US) 

* # ( 0 A T AwWDK ( 1 1 a) , kdOO Y) 

THtRE  IS  A MULE  AT  OAT  ArtROR (115) 

* # (OA I AwRDR ( lib) , WUHIC) 

*, (OATAwKO«( 1 1 7),S) 

*» (OATArtHuHt 1 18) #(,P5lHU) 

*, (DA  TAaHDK  U 19) ,SOUND) 

*,  (OATAwHi>K(i20)»CTMfclAO) 

*,  (OATArtR()R(  l <21  ) ,CTHE1  A) 

* # (0A|  A*KDR  ( U<2)  # i HE.  I AOUT  ) 

*,  (DATA*ROK(  1<23)  »UAI«) 

*,  (OATAWHOK( 1241 ,UYGJfcT) 

*,  (OAT AwRDR ( 125) ,UbOOY) 

* , (OA  I AwRDR ( 1 <26  ) /DEARTH) 

*,  (OA  TAwWOR(  127)  . VAIR) 

*,  (Da  r aahok ( , vbody ) 

*,(OATA*KOK(l<29),VfcANTH) 

*, (OATawROK(130), VELIO 

* , (OA  T AwRDR ( 1 31 ) , VEL  TRUE) 

* # (OAJAwndR(  1 3 <2 ) , h I TCMRAivl 

* # ( 0 A I A wk  DR  (133) ,wHlC) 

* # (OAl  A«Hl)h  ( 130)  , wbODY) 

*,  (DAT  AwRDR ( 135), DEARTH) 

*,  (OA  f AwROHd  36)  »«EH3hI ) 

*, (0A1 AwRDH(137),wHD0T) 

*,(0ATAWkDKd3«)#Xt>AR) 

*,(UATAwR0R(t39),X«KX) 

*,  (OAT  ArtRDH  ( 1<40),R1TCHSTK) 

*,  (DAT AWRDR ( Id ) , XBARNOM) 

*, (OA I AWROK ( ia^) , X800Y ) 

* , (OAT  AwRDR ( 1 U3) /ROLLPAN) 

*, (DAT  AwRDR (1««) , YbAR) 

*, (OATAWRDR  C l aS ) , YBRX) 

*,  (OAT AWRORt 146) , ROLLS  IK) 

*,  (DAT  AwRDR  ( 1 <1  7 ) # YBOQY) 

* # (OA  I A wROR ( 148)  # £ H Ak  ) 

*, (DAT AwRDR l 1 A9  ) # ZBR  X ) 

*,  (DAlAwhOKl  1 50 ) , ZE  I A ) 


*,  t UA  1 ArtRDK  US1) 

* , (DA  I AWHDM  ( 1S2  ) 
* , (DA  I awWOW  (153) 
* , (DA  r AwRUR  C 1 sa ) 
*, (DAI A«KDP(155) 
*, (DAI ArtR0H(156) 
*, (DA  I ArtRDH ( 157  ) 
*, (DAIAWKDR(15«) 
*, ( DA  I AwRDR ( 1 59 ) 
*, (DA  I AwRDw ( 1 toO  ) 
*,  (DAI  AwRDR(  161  ) 
*,  (DATAWHDHU62) 
*, ( DA  f AWKDR ( 1 63 ) 

* t ( D A I AaRQN (1 6A) 
*, (DA  I AwRDK( lbb) 
*, (DA r ArtMDK ( 166 ) 
*, (DAI A*KOW( 167) 
*,  ( DA  I AtvKDK  (160) 
* , (DAI AWKDK (16V) 
*,  (DAI  A aHDW  ( 170) 
*, (UAIA*RUH(171  ) 

* » (DAI A*RDK( 172) 
+, (DAI AWKDK( 1 73) 

* f (DAI ArtKDK (17a) 
*,  (UA  I A«HUH( 1 75) 
*, (DAI  A AROR ( 1 76) 
*, ( DA  T AaRDR ( 1 7 7 ) 
* , (DAT  A aKDR ( 1 78) 
a# (DA | AARDH (17V) 
*, (DATA*WDk( 160) 
*, (DA  I AaHuD ( 181  ) 
*,  (DA  r AaKDK (182) 
*, (DAT  AWHU* (163) 
*, (DATAWRDk( 18A) 
*,  (DAI  AkvRDR  (185) 
*, (DAT  A wRDR ( 186) 
*,  ( DA  I A*vKuR  (18  7) 

* r (DAT  AwRUK ( 188) 
*, (DATAAHDR (18V) 
*, (DATAaHDW ( 1 VO) 
*, (DAT AWKUK( lsl  ) 
* , (DA  I ArtKUR ( 1 V2  ) 
*, (UAlAAKDK(lVS) 
* , (DA  T aaRDK ( 1 VA  ) 
*, (DATArtROK( 1V5) 
* , (DATAWRUK ( 1 V6) 
* , ( DA | AaRDK ( 1 V 7 ) 

(DA  I AaRDK (1 R8) 
* , (UA | ArtHDR ( 1 VV) 


YA *1  RIM) 
ZBAWNOM) 

zauDY) 

FLAPCMU) 

CCD) 

CCE ) 
ccoaFi) 

CCDBF2) 
CL0SB1 ) 
CC0S82 ) 
CLB) 
CLdDtl ) 
CLBDE2) 
CLaosai ) 
CL  BDSH2 ) 
CLdDSb3) 
CLDA) 

CLDR  ) 

CLP) 

LL») 

CMO) 

CmE) 

CMOBF 1 ) 

CMDBF2) 
CMDSH1  ) 
C*0SB2) 
LMD) 

CwB) 
LNdDEl  ) 
LNB0E2) 
LNBDSB1  ) 
CNBDSB2) 

C NBDSb  3 ) 
CNDA) 
C'YDR) 
CNDR8) 

CMP) 

CNR) 

LNO) 

CNDE) 
CNDBF  1 ) 
CndBF  2 ) 
CNUSB1  ) 
CNDSB2  ) 
CTB) 

C YD  A ) 
LYDK) 
CYDRDSbl  ) 
CYDKDSB2) 


* , (DATABY  I ( 1 ),FC3DA) 

* , (DAIABYI  (2UFCSDE) 

* r (DAIA6YI  (3)fFCSDF) 

* , (DA f AbYT («) ,FCSDR) 

* , (DA  TAB Y I (8) ,FCSSB) 

* , (DAIABYI (6), F1VDEG) 

* , (DA  I ABY I ( 7 ) , I KtDtU) 

* , (DA  I AHYT (6) , *IND) 

* , (DA  I ABY I (V) ,Ht)TSET) 


'T  U-U 


* , ( 0 A 1 AbY  Ido;  , HULl)  ) 

* , (DA  T AB  Y I (11), IC) 

* , (DATABY  I (12)  , AdUWT  ) 

* , (DA  I AbY T ( 1 3) , T OKNCUKD) 

* , (OAT AbY K 1«) , AUTOBE ) 

* , (DAT AbY T ( lb) ,DEPULSE J 

* , (DAT  AbY  T ( lb)  ,OAPULi>t  J 

* , (DA  I AbY  T ( 1 7 ) ,DKPULSt) 
*, (ANGLE ( 1 J , ALPHA ) 

*, ( ANGLE ( 2 ) , PH  I ) 

*,  (AnGLE(3), THETA) 

DATA  Pl/3.  141592654/ 

DATA  AEHOSIZE/S243/ 


DATA  SCALuUT/ .6727, 3.441 , .6727,  . 87  2 7 , 1 . 7 45 , . 4383 , 1 2b  . 7 , 6 4 . 35  , 

* 64, 35* 1,745,25, 0,2040.0,32/68. 0,16,0,5. 142,3. 142,3. 142, 

* 102400.  0,1024.0,1024.0,512.0,  10.24, 0.0. 0,0, 3. 142, 3.  142* 

* 35. 0,38.0,38.0, 400000.0,35. 0,30. 0,90. 0,25. 0,4. 096, 9*0.0/ 

DATA  SC AL I N/25. 0,20, 0,30.0,00.0,1. 0,1. 0,1. 0,1. 0,1. 0,1.0, 

* 1.0, <71*0.0/ 

DATA  ALPHA)/ 0.0, 0.0, 0.0, 0.0 893, 0.1 7 45, 0,26 18, 0.3 49 1,0. 4363, 
*0.5236,0.6481,0.6727/ 

DATA  MACHT/O.O, 0,0, 1.5, 2. 0,3. 0,4, 0,5. 0,8. 0,10,0/ 

DAfA  ELEVT /U.O, 0.0, -35. 0,-20. 0,-1 0.0, 0.0, 10.0,20.0/ 

DATA  DSbT/O.O, 0.0,0.0,25. 0,55. 0,87. 2/ 

DATA  BE TAT/o. 0,0. 0,0. 0,2. 0,5. 0,10.0/ 

DATA  ALPHAH T / 1 , 4/ 

DATA  MALHPl/1,7/ 

DA  I A ELt  VP  1/1,6/ 

DATA  DfcLP  T / 1 , e / 

DA  1 A DEPP  T / 1 , 6/ 

DATA  DS8PT/1,4/ 

DATA  bt  T AP I / 1 , 4/ 

DATA  ZtrAPG/O., 10000, ,20000. ,30000. ,40000. ,50000. ,60000. ,70000., 

1 80000. ,40000. ,100000. ,110000. ,120000. ,130000. ,140000, ,150000., 

1 160000., 170000. ,180000. ,190000. ,200000., 210000. ,22 0000. ,230000., 

1 240000., 250000. ,2b0000., 270000. ,280000. ,29 0000. ,300000. ,310000., 

* 320000. ,33000  0 , ,340000. , 350000. , 36 0000. ,370000. ,360000. , 39 0000., 

* 400000./ 


DATA  SND I AbL/1  1 16. , 1 07  7. , 1 037. ,995. 
1 991. ,1004. ,1022. ,1040. ,1058. ,1075. 
1 1020. ,494. ,96/. ,934. ,911. ,884. ,884 

DATA  ZfclAbL/. 292330,. 303149,. 314793 

* ,.34  38  32, -.407019, .417173, .424332,. 

* , .439288, .436588,  .436120, .433096, . 

* ,.420577, .420275, .420996, .422748, . 

* ,.44694 9,. 451434,. 4 54 810,. 457721,. 

* , .46U  1 34, .458847, .457300/ 


,968. ,968, ,968, ,9/1. ,978 
, 1082. , 1 082. , 1071  . , 1058. 

. , 884 . , 884. , 904  . / 

,.32 7/19,. 350427,. 376466 

429406. . 433938. .437363. . 

429438. . 426540. .424377.. 

425734..  431314.. 4 36443.. 

459620..  460466.. 461 186.. 


. ,985.  , 
, 1046.  , 


4 38942 
421935 
441316 
461258 


DATA  AlNDAWG/0., 5000. ,10000. ,15000. ,20000. ,25000. ,30004. ,35000., 
1 UOOoO., 45000. ,50000. ,55000. ,60000. ,65000. ,70000. ,75000. ,80000., 
1 85o00., 40000. ,95000./ 

DA  I A DbLSi/E,K31ZE,  1 S I Zt , H AKS 1 Zt , b Y I S I Zt , S 1 Zt / 1 7 , 400 , 0 , 0 , 1 7,434/ 
LIST  CONTAINS  iHt  ASCII  D0U6LE*UHDS  WHICH  DtPI'Vt  V Aw  1 A8Lt  NAMES  Tu 


EASE  AND  DATASTORE 


DA  I A CIST/ 

REAL  OOUBLErORD  VARIABLES  ARE  DEFINED  HERE. 

* ' DELR ' 

* , 'DELROUf  ' 

* .'GRAVITY' 

* , 'LATDUT ' 

* t ' LNliDOT  ' 

THIS  HOLt  CORRESPONDS  TO  DAT ADBL (6) 

* , • ' 

* t ' PSIHDOT ' 

* , • RUH ' 

* r ' RUHOO I ' 

* t 'UH' 

* , ' HH  ' 

* t ' XE  AR  T H ' 

* , ' XH' 

* t ' YE AR  TH  1 

* , ' YH' 

* / 'DEARTH' 

* , 'ZH'  * 

THE  REAL  WORD  VARIABLES  ARE  DEFINED  HERE. 

* ,'alphao' 

* .'AX' 

* , ' AY  ' 

* , • AZ  * 

* , ' ALPHA  ' 

* , 'ALTITUDE* 

* , 'b' 

* ,'BETAO' 

* , ' Bt  T A ' 

* ,'BRAK' 

* # 'C  ' 

* ,'CC 

* ,'CCBF' 

* , 1 CCSB ' 

* , 'CD' 

* » ' CL  ' 

* . 'CLBT ' 

* ,'CM' 

* ,'CMBF' 

* ,'CMSB' 

* » ' LN ' 

* ,'CNBF' 

* ,'CNSB' 

* , • CN8  T ' 

* , 'CY  ' 

* t ' CLBDE ' 

* ,'CLBSB' 

* ,'CNBOE' 

» , ' RUN  TIME' 

* # • I 1 ME ' 

* ,'CNBSB' 

* , ' C YURSB ' 

* t 'CROSSRNO' 

* , ' DA IL ' 

* .'DEL' 
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* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

Ik 

* 

* 

* 

ik 

* 

* 

* 

* 

* 

* 

* 

* 

* 

Ik 

Ik 

* 

Ik 

* 

* 

Ik 

* 

* 

* 

* 

* 

Ik 

Ik 

* 

Ik 

Ik 

Ik 


. . 


, ' ALTO* 

, 'PHMCSOP  ' 
k ' DEW  ' 
k 'OHIC  ' 

, 'OHUD' 

, ' UttF  I C ' 
k 'DEUC  ' 

, ' 0 ACM  * 

# * DELH1C  ’ 

, 'DELTA' 

, 'DEKiC ' 

, 'DOwNHNG ' 

, 'DSblC ’ 

» 'ELEV 
, ' EwINl) ' 

, 'Flap' 

, 'go • 

, 'GEAH' 

, 'UPJET  ' 

# 'UPJET ' 

, 'HUOT  ' 

, ' IXMY ' 

, 'GEUAL  T ' 
k • lxx ' 

, ' ixr  • 

, ' ixz* 

» 'IYMZ' 

. * I YG  y • 

» ' i y y ' 

, '1ZMX' 

k ' izoz* 

, ' IZZ' 

, • IFUNC ' 
k ' J2' 

k ' KCAS ' 
k 'PEAS' 

, 'KWHCSUP  ' 
k ' KNOTCON ' 

, ' HP  1 ASUP  ' 
k 'LAT ’ 
k 'LATO' 

» ' LNG ' 
k ' LNi»0  ' 
k 'UZCMD' 
k ' UXLMU ' 

, ' LBOU Y ' 
k 'LENGTH ' 
k ' UYCMu ' 
k ' MACH ' 
k 'MASS' 

, 'MU' 

, 'MttOOY ' 

, ' NtiUU  Y ' 
k ' Nrt I NU ' 

, 'ObLA {E  ' 
k 'OMtGA ' 

, 'Pbic  ’ 

, ' PDOT  ' 

, 'PHI' 
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* , • PM  I 0 * 

* , 'PSI ' 

* ,'PS10' 

* ,'PSIH* 

* , 1 PS I P * 

* , 'PbQUY ' 

* ,'PHIUUT* 

* , 'Pfc' 

* ,'UbAP' 

* ,'UBIC' 

* r ' OUUT ' 

* ,'gbooy' 

* , 'PO1 

* , 'PB1C  ' 

* ,'PUOT' 

* , 'PHO' 

* , • PMUO ' 

* , ' S8HP  ' 

* / 1 RADIUS ' 

* t • PBODY  1 

THIS  HULt  CORRESPONDS  fU  DA  l AwPDR  (115) 

* , • PUH 1 C ' 

* , 'S' 

* ,'PSlHO' 

* , ' SOUND ' 

* ,'ThETAO' 

* , ' IhETA  • 
nr  , ' THE1  AOOT  ' 

Hr  , ' U A IP  • 

* , 'urujE  i ' 

* ,'UBODY' 

* ,'UEAPTh' 

* i ' V A I R * 
nr  ,'VBODY' 

Hr  , ' Vt  API  H' 

* ,'VEUC' 

-»  , ' VtlTPUfc  ' 

Hr  ,'P1ICHPAN' 

* , ' ftHlC ' 

* , ' *BODY  ' 

* ,’WEAPTH' 

Hr  ,'rtElOHT' 

* , ' hhHuOI  ' 
nr  , ' XBAP  * 
r*  , ' XBPX  ' 

Hr  ,'PIICHSTK' 
nr  , • XBAPNOM  ' 

* ,'XBUDY' 
nr  t ' POLLP  AN  1 

* , ' YoAP ' 

* , ' YBPX ' 
rt  , * PULLS  1 K 1 

* , ’ Y8UUY  ' 
h-  ,'ZbAH' 

* , ' Zbp  x ' 

* , ' Zfc  T A • 
nr  ,'YAhhTPIM' 

* ,'ZBAPNOM' 
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* , ' IbODr' 

* ,'HLAPC MU' 

THIS  PUlNr,PEFtPfcNCtO  by  THE  0 A T AS  T 0W£  SutJPGuTINfc,  IS  LISMl7<?> 

* ,'CCO' 

* ,'CCfc' 

* ,'CCOdf-l' 

* , 'CCOMF«>' 

* , ' CCOSb 1 ' 

* t * CCOSB<£  * 

* ,'CLb' 

* t ' CLbOfc 1 ' 

* f’CLtiDt*1 

* i ' CLbDSB l 1 

* / ' C LttDSbtJ ' 

* t ' CL0OSB3 ' 


* 

•CLDA ' 

* 

'CLO«' 

* 

•CLP* 

* 

'CLP' 

if 

'CMU' 

it 

'CMt  ' 

if 

'CMOtt£ 1 ' 

* 

' CMOBF  2. ' 

* 

'CMOSbl  ' 

* 

'CMUS82' 

if 

'CMu ' 

it 

' C N 8 1 

if 

' CNbDt  1 ' 

it 

' CNB0E2 ' 

* 

' CNBDSb 1 • 

★ 

' CNbOSb^ ' 

* 

' CNbDSbi ' 

* 

' CNDA  ' 

if 

'C*UR  ' 

if 

'CNDP8' 

* 

' t>P ' 

* 

*CNK  ' 

♦ 

'CNU' 

* 

' CNOE  1 

* 

' C NOBP  1 ' 

* 

'CND0F2' 

* 

' C'NOSb  1 ' 

♦ 

'CNOSBc?' 

if 

' c Yb ' 

if 

' C YDA ' 

it 

' c yok  ' 

* 

' C YUPOSts ! 

* 

'CY0K0S82 

* 

' OE.MAN  ' 

* 

' ESHAPE ' 

* 

'GPS' 

* 

'OEMS' 

* 

' Q ' 

if 

' PHUt)<?b  ' 

it 

1 1 

if 

•GJET ' 

if 

i i 

it 

' t T P 1 M I N ' 
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r ' fc  TP  I M ' 
p 'kLFbKIN' 

- 'ELFBK  ' 

r 'UklMPAN' 
r ' DSBPC ' 
r 'DSBXTP' 

. 'OSBXTPS* 
r 'OfcTUSB' 

. 'UtTKIM' 

- 'GUSB' 

> 'bTRfc ' 

'OfcTMKHC ' 

'Ufc.TP' 

'OEMSP' 

' P T PHI  * 
'PTANPHI • 

' L)  J P TP  * 
'BCSL' 
•UCSL ' 
'OuCT  ' 

' DtCMD ' 
'UFFIL  ' 

•rjcsu' 

'OCSLLO' 
•OCSLHi ' 

•midpick  * 

• K P l T ' 

' GDU  • 

’ UuLU ' 

*P  ' 

' DUHI  ' 

' DUM  1 A ' 

• Oy  h i n 1 

•nz  a ' 

' N Z P * 

' DULON  1 
' OULU  A • 

' ALPHMIN ' 

' ALPHMAX  » 

' 6 T PEDDC  T ' 

• P • 


, • i 

, ' U A f MHrtC  ' 

, • P*UD2« ' 

, ' PPU025  * 

, £* ' • 

# 'PPODZ1  • 

/ 'PPl^^* 

, 'PPUDZi' 

, *PkuDZ/' 

* 'SUM^l ' 

i i 

# ' 3 ' 

, 'DLSP* 

, • DACMD' 

* 't»A  fSUMl  ' . 
, ' KyAML  IN ' 

, ' KOAMPAP  * 

UPS' 


9 


* 

* 


* 

* 

* 

* 

* 

* 

* 

A 

* 

A 

A 

A 

A 

A 

A 

A 


• Gl)  AC  • 

•GUA  * 

1 D A ( SUM  1 
,PKUU42t 
•□AIR1 
•OKPR*1 

• PS  I Ab' 
•SUM^a  • 
•OKHHl  • 
•RHODES1 
•SUM22* 

G T X • 

YAwXFttD  1 
PN0D41  • 

GR  1 
GUM  1 
• 

i 

S JM<?5 ' 
SUM^bS  ' 

• 

SUM£6 * 
PKO0«3 ' 
KSIAB' 
BfcTAU' 
BtTAf-  1L  r • 

OA 1 MP AN ' 

0 A T P ' 

i 

phuoup  * 

OKI  ' 
OKTmSF ' 
OHMS  ' 

OP  IMS ' 
PKUU  1 ' 
auMi  • 

N YP  ' 

NY  A ' 

GK  AY  ' 
PKUD2  ' 
PH003' 

SUM2' 

PkuOb ' 
OKCPF  ' 
HOPC ' 

GO«C  ' 

GNYDKm • 
OKTKW' 
PHUO  1 ' 

G IKK  • 
UKTKIm' 
SUM3* 
OKC.MO ' 
DACMC  ' 
GKXFl>  ' 
l»Okt  ' 

i 

nk  * 

3*'  ' 


m 

* lHE  LOGICAL  BY  Tfc  V AK 1 ABLES  ARE  DEF 1 NED  hERE. 

* ,'FCSDA' 

* i * F LSDE 1 

* / * FCSDF ' 

* ,'FCSOK' 

* t ' F C S S B ' 

* t ' F 1 VDEG 1 

* , 'TxEDEG' 

* , ' W 1 ND ' 

* ,'PorsET' 

* , ' HOLD ' 

* , ' 1C  ' 

* , ' A BUR  l ' 

* , ' 1 URnCURD  ' 

* ,'AUFOBF' 

* » ' OEPULSE ' 

* , ’OAPULSE  * 

* f * DKPUL St ' / 

•EXECUTIVE  PORTION  FULLUWS 

•THIS  IS  THE  EnIRY  PUINI  FUR  PRUGrAM  ACTIVATION 

** 

•INITIALIZATION 

* 

• THIS  SECTION  INITIALIZES  THE  INTERVAL  TIME*  HANDLER. 


) 1020 


I nL l Nt 


01 

43 

disable  interval 

T IMEK 

L* 

7.4ZB 

FETCH 

ADDRESS  OF 

ULO  HANDLER 

S T m 

7 t HANDLER 

STORE 

FOR  LATER 

USE 

LM 

7 t ) 1000 

FETCH 

ADDRESS  UF 

nEw  handler 

STr 

7,4*8 

STORE 

ADDRESS  IN 

handler  LOC 

E NO  I 


THE  AERODYNAMIC  OAIA  IS  PLACED  In  CORE  HERE 

II*l 

JJ=1 1S2 

3 CALL  BUFFERING, 1,0UTDATA(II),JJ) 

4 CALL  STATuSCS, TEMPI) 

IF ( TtMPl ,nE.2)G0T0« 

11*11*11 S4 

IFdl.GT.AEROSIZEJGO  TO  5 
IF  ( Il  + USZ.GT  .AERUSIZE)JJ*AEKOSIZE*l-ll 
GO  TO  5 
b CONTINUE 

1C*. T RUE . 

POTSET*. FALSE. 
huLDs.F ALSt. 

THE  F OLLOR 1 NO  VARIABLES  ARE  SET  10  INDICATE  A DESIRED  CONDITION 
AS  INDICATED.  THESE  VARIABLES  ARE  EASE  ACCESSASLE. 

AUTO  BODY  FLAP  (NUJ  AVAILABLE) 

AU100F*. false. 

TURN  COORDlNAl ION 
TURNCORDs.TRUE. 
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ORDERL  Y PROGRAM  f E H M X N A I I U N 

abor is. false* 

FIVE  DEGREES  OF  FREEDOM 
F IVDEG*. TRUE. 

THREE  DEGREES  OF  FkEEOOM 
TREUEG*. false. 

ATMOSPHERIC  «INO 
WIND*. FALSE, 

CONTROL  OF  VEHICLE  SOKFACE  OEF LE L T I ONS . . . . I RUE . INDICATES  THAT  T Ht 

VALUES  ARE  DETERMINED  BY  DEL  1 C , UER I C # OR  I C , USB  1 C , DBF  I C . 

F C SUE sFCSuAsFCSURs. false. 

F C SOF  *F  C SSB* , TkuE. 

surface  pulse  for  oynamic  response  chec*....true.  indicates  that 
THE  SURFACE  is  fo  be  pulsed. 

DEPULSE=UAPULSEsORPULSE*.F alse. 

PRESET  PROGRAM  RUN  IIME  IN  SECONDS, 

RUNT IME  *5 . E 75 

C ALRH  A0*20  . U 

CBE  T A0*0 , 0 

CThETA0*20.U 

CPhIOsO.O 

CPSIusO.O 

CL  A T 0*0 . 0 

CLNGOaO , 0 

CPSIK*«!6.0 

S*2b9u,0 

B*/8.o3 

C*39.5b 

LENGTH* 10  7,5 

wE1GhT*155UOO,0 

iXXsbOOGOU.O 

1 Y Y *550  0 0 U 0 , 0 

122=5500000.0 

1 X Y * 0 . 0 

1X2*100000.0 

XBAKiO ,6b 

YB AR*0 , 0 

ZBAR*0.290  7 

X6AHN0M*0.6S 

ZbAHNUMsO.2907 

GO*  32 . 14b54b 

J2*. 0010823 

KNUTCONaO, 59210526 

UHL A TE*  7 0 1 SO , 0 

OMEGA*.  000072«2U59 

RO  *20  925  7 38,0 

RHUOs. 0023769 

DBF IC*-5.0 

DEL IC*« 1 6 . 33 

DERIC*-lb.33 

DRICsO.O 

DSB 1 C *25 , 0 

Al  T0  = 49a80.u 

PBICSO.O 

UBICsO.O 

RB I C*0 , 0 

v El  I L *29  7 0 , o 

DELTA*. 04 
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HALFP1  s Pi/2 
TWOPI  s 2 * P 1 

a 

* T Mi  S IS  THE  EASE  HtlUHN  PUiM 

10  I F ( A dUP  I ) GCJ  TO  1100 

* 

* A OI&ifAL  INPUT  IS  PEPFUPMEo  HEPE  TU  INITIALIZE  THE  MODE  CunTpOL 
a PAPAMEIEPS, 
lb  CONTINUE 
INLINE 


Lw 

4,BZ00o3FF'FF 

LEA 

1 .7920 

SI  MW 

1,256 

CD 

O,4ZSOQ0 

E NO  1 

* 

MASS*w£I&HT/32.174 
pe*po-ohla I t 
HALFdELsDELTA/2 

* 

•PI TSl ,PI TS2,PI f 3 awe  TEMPOPAPY  LOCATIONS 
PI  I 1=MASS*LENG1M 

XbPX=(XbAKNUM-XdAP)APlT  1 

YdPXsYBAKApl  T 1 


Z6*X  S ( ZdAKNOM-ZdAP) API  T 1 


I 


1 XMY=1XX-I YY 
1 YMZ=I YY-iZZ 
IZPKSIZZ-IXX 
IYUYsIXY/  I Y Y 
1Z0Z=IXZ/IZZ 

IFUNC«1/ ( IXX-IZUZAIXZ-IYQYAIXY) 

SC»SaC 

SCZbY2s.b*SACA*2 

Sfj*S  Ad 

SB2BY2s.5»S*B**2 

* THE  FOLLOWING  SECTION  CALCULATES  T mE  INITIAL  CONDITION  VALUES 
a FOP  IME  INTEGPAL  EUUATlUNS.  EASE  DEFINES  f FiE  T A , PH  I , PS  I , ALPM  A , 

a BETA,VELTPUE  AND  ALTITUOE.  THESE  VALUES  APE  PPOCESSED  TO  OBTAIN 

a PUH,PS1H,WH, AND  DELP. 

A ALL  ANGLES  APE  1NTEPNALLY  DEFINED  IN  PAD  I ANS . ALL  COMMUNICATIONS 
a WITH  LAKDS  AND  EASE  APE  WITH  ANGLES  DEFINED  IN  OE&PEES . A "C" 

A PPECEEDS  ALL  UEGPEE-DEFINEO  ANGLES. 

THE  T AOeC  THE  I AO/57. 295779b! 

PHI 0=CPH 1 0/57. 2957 7951 
PS 10=CPS 10/57. 2957 7951 
ALPHA o*C AlPHA 0/57. 295 7 79b 1 
BET AOsCHE I AO/57 ,295779bl 
PS  I P=CPS IP/ 5 7, 2957 7951 
L AT 0*CLATO/57. 2957 7951 
LNG0=CLNG0/57.29577951 
PI T 1=S1N( The  TAO) 

PIT2sC0S( THETao) 

PITiaSlN(PHlO) 


PI  UsCuS(PHlO) 


PI15»SIN(PSiO) 
PI Td=COS(PSIO) 
L1*PI T2*PI Tb 
L2*P 1T2aP1T5 
L3»-P1  I 1 


A 
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11 11  m*1'.  ” 


Ml a-PI T tt*PI IStPI T3*PITl*PITb 
M2sPl la*PI l6tPI13*PlTl*PlT5 
M J*P 1 T 3*P  1 12 

N»«P1  T3*PI  T‘>  + P1T«*P1T  1*P1  lb 
N2s-PlT3*PlTbtPlT«*Pll  1*P1 lb 
N3»P 1 T2*P  1 1 R 

DELKIC*ALTO-RO*( , 00lb7b28-.001b761b*COS(2*LA TO)-.  OOOOV2U  * 
* CUS(R*LA10)) 

RADsROfOELRK 

U0UDYaVEHl*COS(  ALPHAt))*COS(bE I AO) 

VbOOY*VELlC*SlN(bE  TAO) 
wBUUYsvELlC*SiN( ALPH AO ) *COS ( BE T A 0 ) 

UAIH«UBODt  *Ll  + VbUDY *Ml  +rtbUDY*Nl 
V AlRsUbODY*L2* VbODY*M2+wBUUY*N2 

wEARl  H*ubuDY*L3*  VBQDY*M3+wBUOY*N3 
UEAKIHsUAIH 

veawth*vaik+rau*cuS(latu )*omeua 

PSIH0*ATAN2( VEAWTH.UEAR  TH) 

RUHIC*RAD*UEAR I H/COS ( PS I HQ ) 
wHICswEARTh 

CPSlHO*PSIHO*S7.2VS77Rbl 

* 

* I H 1 S SECTION  STARTS  THE  INTERVAL  TIMER  OPE.RAI1UN. 

* 

C0UNT*DELTA/.00u0012 

INLINE 


Lw 

O.CUUNT 

LOAD  FRAME  TIME 

CD 

127. bO 

START  TIMER 

El 

«3 

tNAbLE  INTERVAL  TIMER 

ENDI 

a 

•PATCHES  WILL 

00  HERE 

*E ND  OF  INITIALISATION  ROUTINE 

* I Ml S IS  THE  LOOP  RETURN  ENTRY  POINT 

u 7 PUTSETsINblT  u ) 

ICslNblT (2 ) 

HOLO= 1 Nb I T 13) 

* 

IFtPUtSEI)  00  T u 1 1 OS 
I F ( HOLD ) 00  TO  170 
IFCIC)  TIMEsO 

* 

* ANALOO  INPUT  IS  STAHTED  HERE. 

INLINE 

L*  *4,  B20003FFFF 

lea  i,mo 

STM«  1 t 2Sb 

Cu  0 » U26000 

END  1 

• 

* THIS  SEC  I ION  TRANSFERS  THE  INPUI  VALUES  CONTAINED  IN  THE  DATA 

* HALFWORD  TABLES  TO  THE  appropriate  weal  WORD  variables. 

JJ*t 

INLINE 

BL  ) RbO 

ENDI 

PRHCSOPsDA I A 

J J»2 

INLINE 


BL  I960 

ENDI 

KRMC3UP*UA1 A 

JJ«3 

INLINE 

BL  ) 96  0 

ENDI 

RP  T A SOPsD A I A 

JJ*0 

INLINE 

BL  I960 

ENDI 

SBHPsOA  T A 

JJ*3 

INLINE 

ttL  ) 9b  0 

E Nl)  I 

PI ICMPANsOA! A 

J Jsb 

INLINE 

BL  ) 9b0 

ENDI 

PI  I CMS  TK»D A 1 A 

JJ*7 

INLINE 

BL  )9t>0 

ENDI 

HULLPANsDA  T A 
j Jsb 
INLINE 

BL  ) 960 

END! 

ROLLS  f**DAIA 

JJ39 

INLINE 

BL  I960 

ENDI 

YAhTRImsDA 1 a 

JJ«10 

INLINE 

bL  I960 

ENDI 

FLAPCMDsDA  i A 
JJ*1  1 
INLINE 

BL  ) 9b0 

ENDI 

PuLSEsO A I A 

IF(IC)  GO  U)  «9 
GO  ID  170 

«9  PBUuYaPBIC IUBODYsUBIC; RBODYsRBIC 

T ME  1 A = TmEIAO  J PM  I ■ PM 1 0 > PS1  S PSIO 

*hs*mic;uelrsdelmic 

PSIMsHSIMO; MOMaRUMIC 
LASLA rOJLNsLN&O 
170  C0N1INUE 


C IMIS  S1AWTS  THE.  DA  f A LOOK-UP. 

C FIND  AHGUMEN  f BKt AKPO I N I S 
ABSbtl AaAbSIbt  I A) 

CALL  P0lNfF(MACM,MACMPT(2)#MACHPT (1), MACH  I U),MACHT(3)) 

CALL  POINTFl ALPHA, ALPMAPI (^) , ALPHAHT (1 ) # ALPHA! 1 1 ) , ALPMA1  ( 3)  ) 

CALL  POIMF  (fcLtV/tLtVPI  (2)  #ELEVPT(I)#EL£VF  ( 1 ),ELEVT (3) ) 

CALL  PUINTF (DEL,0ELPT(2) ,DELPT ( 1 ) • DEL T ( 1 ) ,ELEVI  (3) ) 

CALL  PUlNfF(OEK#DEKPT(2)#DEHPT(l)»OERT(l)»tLEVT(3)) 

CALL  MU  IN  IF IBkAK.DSBPT 12) ,OSbPT ( 1 ) , DSB T C 1 ) , USHT  ( 3 ) ) 

CALL  POlNTf  (Ab8BETA,BETAPT(2)fBETAPTU)»BETATll),BETAT(3)) 

C FIND  UEHlVAiivt  VALUES 

CMO»DERIVt2lCMO!,ALPhAPT (2) ,MAChPr (2) , ALPhAPTU ) ,maChPT  C 1 ) > ALPHA  T ( 
*1)#MACHT(D) 

CNUsDEPl V£2(CNOT, ALPHAPT (2) ) 

CNOt*OEHlV£2(CND£T, ALPHAPT (2)) 

CC0*DEKIVE2CCCUT, ALPHAPT (2) ) 

CLP»DfcWI VE2ICLPI /ALPHAPI (2) ) 

CLH*DEKIVE2(CLKI , ALPHAPT (2) ) 

CNDKbsDEKlV£2(CNUKbT .ALPHAPT (2) ) 

CNP*DEP1VE2(CNPT, ALPHAPT (2) ) 

CNH»DEPlvt2(CNK!, ALPHAPT (2) ) 

C YbsDEMl VE2(C YBI , ALPHAPT 12) ) 

C YDNsdEPI VE2(C YDK r , ALPHAPT (2) ) 

IF  (FLAP. Gt .0.0)  GOTU  220 

CCbF*tCUBF 1 *UtP I VE2 CCCDbF 1 1 , ALPHAP H 2 ) ) 

CMBF  3CMDBF 1 autK I VE2(CMOBF 1I,alPhAPT(2)) 

CNBF st NDBF 1 *OEH 1 VE2CCNUBF1 T, ALPHAPT (2) ) 

GOTu  230 
220  CONTINUE 

CCbF=CCUBF2»DEPIvE2(CCDdF2T, ALPHAPI  12)  ) 

CMbF*CMDbF2»DERI VE2(CMDBF2T, ALPHAP I (2) ) 

CNbFsCNDdF  2*DEH I VE2 ( CN0BF2 T , ALPHAPT  (2) ) 

230  CONTINUE 

IF  (bPAK.GT .25.0)  GUT0240 

CLdSbsCLBUSblsDtHl VE2CCLBDSblT,ALPHAPTC2)) 

CNbSbsCNBOSB isOEPl VE2(CNBDSB1 T , ALPHAPT  C2) ) 

GOTU  245 
2 A 0 CONTINUE 

CL«Sd*CLbDSB2*DEKI V£2(CLBUStt2T,  ALPHAPT (2) ) 

CNBSb*CNbDSB2=DEHI VE2(CNBDSB2T , ALPHAPT (2) ) 

245  CONTINUE 

CLdDSo3*CNbDSdi=0.0 
IF  (dWAK.Ll ,60.0)  GOTU  240 
CLHDSB3=0EK1 VE2(CLdDS83T, ALPHAPT (2) ) 

CNBDSB3=DEH I VE2CCNBDSB3T# ALPHAPT (2) ) 

249  CONTINUE 

IF  IELEV.GT.0.0)  GOTO  250 

CLB0EsCLdUtlaUERIVE2(CL6DElT, ALPHAPT (2) ) 
CNBDEsCNbDElsDEHIVE2(CNBDEl T , ALPHAPT (2) ) 

GOTU  255 

250  CLBUEsCLbUt 2sDEK I VE2(CLdUE2T» ALPHAPT (2) ) 

CNBDE=CnBDE2=DERI VE2ICN6DE2T/ ALPHAPTI2)) 

255  CONTINUE 

IF  (BHAK.GT .55.0)  GOTO  258 
CCSBsCCUSolsDEWi VE2ICCUSB1 T , ALPHAPT! 2) ) 

CMjjbaCMUSHlsDEPI  Vt2(CMUSbl  T , ALPHAPT  (2)  ) 

CNSbaCNDSdI»UEWlVt2(CNDSBlT , ALPHAPT (2)) 

CYDHSbsCYDPDSblsDEKI VE1 C C YuRSb 1 T , M ACHP T (2) .MACHPT ( 1 ) ,paCh1 (1 ) ) 

GUIU  259 
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258  CONTINUE 

CCSbaCCDSBBaUtRI  VE21CCOSB2T,  ALPHAPU2) ) 

CMS«stMDSb2sutKl Vt2(CMDSb2T , ALPHAPT (2) ) 
CNSHaCND3b2at)EKIvE2(CNUSB2T,  ALPHAPT (2)) 

Cn>K3b«Cn>WUSb2*L»fcHI  VfcJ  (CYDRSB2  T .MACHPT  (2)  /MACHPT  U ) /MACHT  (1)) 

259  CONTINUE 

CMELaOER I VE3(0CM£T, ALPHAPT (2),MAChPT(2)/ DELPT (2), ALPHAPT (l ), 

* MACHP1 ( 1 ) .DELPT 1 1 ), ALPHA! U ), MACH T U ), DEL T ( A) ) 
CCELaDEKIVE3(DCCET, ALPHAPI (2) ) 

CMEh*OE«I VtilOCMtT » ALPHAPT (2), MACHPT (2), DERPT (2) . ALPHAPT (1), 

* macHPT(1),OEHPI Cl) » ALPHA! (l)#MACHl (l).DERT(l)) 

CCEkauERI VE3(DCCET, ALPHAPT (2)) 

CLDAaOERI VE 5 (CLDAT .AlPhAPT (2) .MACHPT (2) .ELEVPT (2) , ALPHAPI ( 1) , 
*MACHPT(1),ELEVPT(1),ALPHAI ( 1 ) , MACH T ( 1 ) , ELE V T { 1 ) ) 

CNOAsDERIVE3(CNDA I , ALPHAPT  (2)  ) 

CTOAaOER I VEitCYOAT, ALPHAPT (2)) 

CLbaOtHlVti(CLBTAbL» ALPHAPT (2) #MACHP1 ( 2 ) , bt I AP T ( 2 ) , ALPHAPT  ( 1 ) 
•.MACHPT ( 1) ,tttl APT  U ) . ALPHA  T ( 1 ) .MACHT  C 1 ) , BE T A 1 ( 1 ) ) 

CNbsOEPI Vt 3 (CNHTABLf ALPHAPT (2) ) 

CMuaDtKl Vt3(CMUT, ALPHAPT (2) , MACHPT  C 2 ) , OSbP ! ( 2) , ALPHAPT ( 1 ) , 

•MACHP  f ( 1 ) , USBPT ( l ) , ALPHA  I ( l ) .MACHT ( 1 > , OSBT ( 1 ) ) 

CLUKaoEKI YE  3 ( CLOP T . ALPhAP 1 (2) ) 

CNOKaUtKlVt3(CNu»HT  .ALPHAPT  (2)) 

CMta(CMfcL+CM£K)/2.0 
CCts(CCELtCCEH)/2.0 
C THIS  CONCLUDES  THE  DATA  LOOK-OP. 

* 

•CALCULATION  of  DERIVATIVE  coefficients 
•CL  and  LbODY 
•CALCULATION  of  llbt 

CLbT  aCLb*CLBDE*ELE VfCLbSb*  ( BR  AK-2S ) +CLBDSb3*  ( BK  AK-t>0  ) 

•CALCULA1 ION  OF  LbODY 

LbUDraubAK* ( (CLH I *B£ 1 AtCLDA*DAIL ♦CLDR*DKUU) *Sb 

* ♦Sb2BY2*(CLP*PbU0Y+CLR*RB0DY)/VELTRuE)-YB«X*ZBU0Y*ZBRX*YB0DY 

* -uZCHD*6123 

*CM  and  MBUOY 

•CALCULATION  uF  cm 

CMaCMO  + CME^CMBF  *F L AP tCMSB* ( BK AK-55 ) 

•CALCULATION  OF  MttUDY 

MB0dYsuBAk*(CMasC+CM0*SC2BY2*uBUDY/VELTRUE)-ZBRX*XB0DY+XBRX*ZB0DY 
*C N AND  NHUDY 
•CALCULATION  of  cnbt 

CNb  f aLNb*CNbOE*tLEV*CNBSB*  CBRAK-25) ♦CNBDSB3* (BRAK-bO) 

•CALCULATION  OF  NBODr 

NbODY*UBAR* C (CNb! *BE 1 A + CND A*D A I L* ( CNORtCNUKb* AbSBE 1 A ) *DROU ) * SB t 

* (CNR*RBOUY  tCNP«P60DY)*8B2bY2/ VEL T RUE ) -XBRX*YBODY +YBRX*X BODY 

* •fUZLH0»339b9 
•calculation  OF  CC 

CC=CCU*CCE*CCBF*FLAPfCCSB*(bRAK-35) 

•CALCULATION  OF  CY 

CYaCYb*bETA+CYUA*UAILtCCYDR+CYDRSB*(bRAA-55) )*DRUD 

* CALCULATION  OF  CN  (NORMAL  FORCE  COEFFICIENT) 

CNSCNO  + CNUt*tLt v+CNBF *F  LAP tCNSB* (BR Ak-55) 

• CALCULATION  OF  CL  (lift  COEFFICIENT) 

CL  a CN*C08ALF  -CC*SINALF 

• Calculation  of  cu 

Cu  S CC*CU3ALF  ♦ CN*SlNALF 

* 


* 1 Hfc  POLLUTING  lb  A f-UUKlM  ORDtR  RUNGt  KUITA  INTEGRA!  1(JN  RUUIINE 
*FUR  FORMATION  UF  HbuD  Y , UBUU  Y , AND  RtiODY 

*Th1S  StCUON  UbtS  FUNCTION  TYPE  SObRKUGHAMS  FuR  CAlCULAUUn  OF 
*DEW  I V AT  1 Vd  TERMS 
*DELTA=FKAMt  TIME  IN  SECONDS 

* 

ROOT *POOf  T ( P B 0 D Y , U B Q D Y , R B 0 D Y ) 

KOaOEUT A*RDUI 

UDOf  auDOT  T (PBODY,UBODY  ,RBOUY  ,PU0! ) 

COaUEU I A «UDU T 

RDUTskUUTT  (PbuDt»UtJUOY*HdGDY»PDUI  ) 

MOaOEL  T A *RyUT 


P I T 1 *POO  T T (PbOUY+,5*KO,UdUDYt.S*LO,RHODYt.5*MO) 

K 1 *OtL  T A *P  1 T1 

C 1 *DEC  T A *QDU  T T C PdOO Y * . b*K0 , UdOD Y ♦ . S*L0 , RBUO Y ♦ . 5*M0 , P IT  1 ) 

Ml aOELT  A*HOUT  T ( PBOD Y t . b* K 0 , QdOO Y ♦ . S*L 0 , RduD Y ♦ . S*MO , P I T 1 ) 

* 

PITl=PDOT  T (PbOOY  + .B*M , U«OD Y ♦ . 5*L 1 , RbOOY+.S*Ml ) 

K2sDEL  T A *P  I I 1 

L2»0EL  | A* QUO  T T (PHU0Y  + .5*K 1 , UBOD Y ♦ . 5*L 1 , RoUD Y t . 5*M 1 , P I T 1 ) » 

M2aD£U  A*hOOT f ( PBOO Y ♦ , b** 1 , ObOO Y ♦ . 5 *L 1 , RBOO Y ♦ . b*M 1 , P I T 1 ) 

* 

PI TlsPOOT  T (PdUUY+K2,GdODY+L2,KdODY+M2) 

*3*0ELT A*PI 1 1 

L 3*0t  L T A * uDO  T T (PBODY«’K2,OdUDY  + L2,RBOOY  + M2,PITl) 
M3sO£LTA*KDOrT(PeOOY+K2,OdOOYtL2,K«ODYtM2,Pin ) 

* 

* S T OP  INTEGRATION  IF  IC  OR  HULO  MODES 

I F ( IC , UR  « HULu ) GO  TO  210 

* 

PBOUYaPbODY tKO/b+K 1/3+K2/3+K3/6 
wHODYsUHODY+LO/b+t l/3+L2/3+L3/b 
RBODY*RBODY*MO/b+Mi/3+M2/3tM3/fe 

* 

210  CONTINUE 

* 

•THIS  SECTION  uStS  THE  FOURTH  ORDER  RUNGE  KUTTA  TECHNIOUfc  TO  GENERATE 
*ThE  EULER  ANGUES  THETA, PHI,  AiyD  PSI 

»fht  RESULTANT  VALUES  FOR  T HE T A , PH  I , A NDPS I ARE  IN  RADIANS 

* T HE  SIN  AND  COS  VALUES  FOR  THE  EULER  ANGLES  AhE  ALSU  CALCULATED 

SINTHETAsSlNl THt TA) 

COSThET AsCUSI ThE  TA) 

S I NPM I =S I N ( PH  I ) 

COSPhIsCUSCPHI) 

SlNPSIrSIN(PSI) 

COSPSlsCOSIPSI  ) 

* 

THE  TADUT=UbUOY*COSPHI-RBODY*Sl NPHI 

koo«delta*the  r adot 

PSIDUT* IKBOU Y*COSPHI*GBODY*SINPH I ) /COS THE  I A 
LOO*DELTA*PSIDOT 
PHlDOI=PdODY*PSiDOT*SINTHEl A 
MOOsDElT  A*PHIDOT 

• 

K01=DELT  A* the  1 DOT  T (PHI t.b*MOO ) 

PIT 1=PSID0TT (PHIt,5*M00, THETA*. S*KOO) 
loi*delta*pi n 

M01aDELTA*PHIDOTI(PHIt.b*MOO,THETA-v,5*KOO,PITl) 
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K02=DELT  A*  I ME  I DUTT (PH1+.S*M01 ) 

PI  I 1*PSI00T  1 (PHl+.S*M01  , THfcl  A+.S*K01  ) 

L02sDELT  A*P  I T 1 

M02*DELTA*PHIOUI l (PH It ,B*MOl , THE  T A + . $*K 0 l , PI T1 ) 

* 

K03aDELTA*TrtETDUTT  (PHI+M02) 

PI T isPSIDO r T (PH1+MQ2# l ME 1 A+K02) 

LO  3*DEL  T A*P  I T 1 

M0  3si)fcLTA*PHIOUT  I lPHltM02»  THETA*K02»PIT  1 ) 

* 

• INI  EGHA  I ION  IS  NOT  PERFORMED  IN  THE  1C  OH  MOLD  MODES 

IF  ( IC  .UH.HOLOJ  GO  JO  270 

* T HHEt  DEGREES  UF  FREEDOM  MODE  HEQUIHES  PhIDOT  ANO  PSI  DOT  sO 
IF ( 1REDEG)  L00=L01*LO2*LQ3sM00sMO1sMO2xMO5sO 

* 

THETA=IHETA+KOO/6+KOl/3+KO2/itK03/6 
PHlsPHl tMOO/6+MO I /3+M02/3+M0 J/6 
PSI*PSI+L00/6tL01/3»L02/3+L03/6 

* 

270  CONTINUE 

* 

• CALCULATION  OF  GHAV I T Y» INCLUDING  E AH  IMS  OBLATENESS 

MU=HALFPI-LA 
COSMUsdCOS (mu) 

SINMUSUSIN(MU) 

•IMIS  TERM, SHARE,  IS  USED  BOTH  HERE  ANO  IN  THE  CALCULATION  OF  XEARTH 
*J2,GC0Nl , ANDGCUN2  AHE  CONSTANTS 
ShARE=-3*J2*GO*(HO/RAD)**9 

GHAV I T YxGO* ( HO / R AD ) • *2 ♦ , 5* SH ARE  * ( 3*C0SMU**2- 1 ) -UH**2/H AD 
•GRAVITY  HAS  UNITS  FI/SEC**2 

* 

* BODY  AXIS  ACCELERATION  CALCULATIONS 

DPI T 1 sUBAH/MASS 

•TEST  FOR  DtSIRED  FIVE  DEGREES  UF  FREEDOM 
IF  (FIVDEG)  GO  TO  290 
XBs-DPl I 1 ‘CC • S 
GO  TO  320 

290  DPIT2s(RB*ZB+VB*YB)/UB 

DPI T3sSlNThETA-CUSTMET A* ( w9*C0SPH I ♦ V8*S I NPHI )/UB 
xrtsGRAv I T Y *OP 1 T 3-DP IT  2 
320  AXsXB 

XBODYsXB 
YBsDPI T 1 *CY«S 
A Y a YH 
V HUD  Y s Y B 
ZBa-ijPI T I «CN«S 
*1*10 
ZBODYsZB 

* EARTH  ANO  H-FHAHt  AuCELERA  T IUnS 

DPIFi=ZB*CuSPHl>YB*SlNPHI 

DPIT2=-ZB*SlNPHi>YB*C0SPHI 

DPI T 3*XB*CUSTHE  T A ♦DP I T 1*S IN THETA 

XE AHThs-DP I (2*S 1 NPS ItDP I T3*COSPS I* SHARE *COSMu*S I NMU 

YEARTHXDPIT  2*COSPSI*DPI T3*SINPS1 

ZE AHTH2-Xh«SlNTHE  T A top  I T ! *COSTHE T AfGH A V I T Y 

XH  s xEAR I h*COSPSIH  ♦ YEARTH*SINPSIH 

YH  e -XEAHTH*SINPS1H  ♦ Y E AR T H*COSPS I H 

ZH  = ZEAHTh 


* THE  FULLOwlNO  SECTION  CALCULATES  IHE  H-FkAME  DERIVATIVE  VALUES. 

RUHDUT = N AD*XH 

PSlHUUIavE*SINLAT/RCOSLATtYH/UH 

«vHl)s^H 

WHDU  F :»H|) 

LA  I DU  1 *Ut/K AD 
LNGDU  r»vE/RCOSLA I -OMEGA 
•OMEGA  IS  THE  EAHImS  ANGULAR  VELOCITY 
•FOR  FIVE  DEGREES  UF  FREEDOM  DELRDOT  MUST  BE  2ERU 
DELRUO T s-rtH 
lF(FIVDEG)  UELRDOT  *0 

* 

* THE  H-FRAME  DERIVATIVE  VALUES  ARE  INTEGRATED  BY  A FIRST-ORDER 

* EULER  InIEGRAFION  ROUTINE. 

•INTEGRATIuN  IS  BYPASSED  In  IC  AND  HOLD  MOOES 

* 

IF ( I C . UR . HOLD ) GO  FO  570 
RUH3Ruh+-RUHD0T*0£L  FA 
PSIH=PSIHtPSIHDOT*OEL F A 
WHsAH-fWHD^OELT  A 

la»la  + la  F DO T *OEL T A 

LN«LN+LNGUUF •DELTA 
OELHsDELR+DELROOT* DELTA 

* 

570  LAFsLA 
LNGslN 

rad*ro+delr 

RAOiUSsKAO 
UHsHUH/R AD 

* 

•THE  FOLLOWING  SECTION  CALCULATES  TRANSLATION  OF  THE  BODY  PROJECTED 
•ONTO  THt  EARTHS  SURFACE, IN  FEET. 

•LATO  AND  LN&O  ARE  REFERENCE  LATITUDE  AND  LONGITUDE  UF  THE  STARTING 
•POINT  OF  FLIGHi 

•PSIR  IS  THE  AnGlE  IN  RAOIANS  BETWEEN  TRUE  NURTH  ANO  TOP  OF  MAP 

• POSITIVE  R S I R IS  WITH  MAP  HOTAFED  EAST , (CLOCKWISE)  FROM  NORTH 
•NOTE  THAI  [nt  CUSLaI  TERM  ACCOUNTS  FUR  CHANGE  IN  DISTANCE  BETWEEN 
•LONGITUDE  LINES  AS  A FUNCTIUN  OF  LATITUDE 

SINLAT=OSIN(LA ) 

CUSLA F sUL  US  (LA) 

RCUSLAT  xRAD*COSL AT 
SINPSIR=SIN(PSIR ) 

CUSPSJRsCUSTPSIR J 

DOwnr  NbaH  0 • ((LAT-lA!0T*CUSPSIR+ILNG-LNG0)*SINPS1R*CDSLAT) 
CROSSRNGsRu*  ( ILAT-LATO)*SINPSIR-»(LNG-LNGO)*COSPSIR*COSLAT  ) 

ALT  IIUDf.  *L'ELRtR0*(, 001678^8-.  0 01o761G*C0SCB*LAT)-.00000^I1* 

• COS 1 4*L A T ) J 
COSPSIH=DCUS(PSIH) 

SINPSIH*0S1N(PSIH) 

Ut  = UH*COSPS  IH 
UE AR  T nsuE 
Vt»UH*SINPSiH 
Vt AR  T m* VE  . 

WE*wH 
w£  AR  T hswE 

* 

• WINDS 

* NORTH  wlNDS(NwlNU)  ANS  EAST  wlNO  ( E w I NO ) Ah E TABLE  LOOKUPS 


* AS  A FUNCTION  OF  ALTITUDE  AND  AKE  added  TU  CORRESPONDING  EARTH 

* AXIS  VELOCITIES  TO  GET  AIK  SPEED 

IF (.NOT .WIND)  GU  10  41 0 
WINDALI  a ALTITUDE 

IF  IwlNDALT.GT.  9BOOO)  wINUALT  » 95000 
CALL  POINTF(wINDALT,20,PLACE,FRACTION,wINDARG) 
NwIND=NwTABL(PLACt )‘(NwTABL (PLACE* 1 )-NwTAbL (PL ACE) ) ‘FRACTION 
EwINUsEwTABL (PLACE )*(£WTABL( PLACE* 1 )-EwTAbl (PLACE )) ‘FRACTION 

GO  TO  420 

410  NwlNu  s EwIND  s 0 

* EARTH  AXIS  AND  AIR  VELOCITIES  CONSIDER  H-FKAME  VELOCITIES, 

* WINDS,  AND  ROTATING  EARTH 
420  UAsUE-NRINU 

OA IRsua 

VAsVE-EwINO-RCOSLAT*OMEGA 

VAIRsVA 

HDOI  = -RH 

*HDOT  REFERENCES  POSITIVE  UPWARDS 

* 

* calculation  OF  bOUY-AXIS  velocities,  alpha,  HE T A 

* SIGN  CONVENTIONS  FUR  VELOCITIES  ARE,  U IS  POSITIVE  FOR  FORWARD 

* V IS  POSITIVE  RIGHT,  AnD  w IS  POSITIVE  DOWNWARD 

DP  IT 1=VA*C0SPSI-UA*SINPSI 

DPI T2=UA*C0SPSI + VA*SINPSI 

DPI  I 3 = w£*CUS  THE  T A*uPI T2*SINTH£  TA 

UB3DPI 12‘COSThE I A-wE*SINTh£TA 

U60D  Y = UB 

VB3DPiTl*COSPHI*UPIT3*SINPHl 

vbuoy  = vb 

wbsDPi J 3*CUSPHl-uPI T 1*SINPHI 

WBUOYswB 

* ALPHA  IS  ANGlE  OF  ATTACK  IN  RADIANS 

ALPhAsATAN2(wBUDY, UBUD Y ) 

SINALFsSlN( ALPHA) 

COSAlFsCOS(ALPhA) 

PIT!  * SORT  (UBUDY**2  * wdODY**2) 

* BETA  is  ANGLE  OF  SIDESLIP  In  DEGREES 

BET As AT AN2( VBODY,PI Tl) 

BE  TAsBETA*b7. 29577951 

AVELTRUE  is  true  airspeed  in  ft/sec 
aTEST  for  FIVE  deghee-of-freedom 
IF'(FlvOEG)  GO  Tu  a30 
VELTRUE3SuRT(VBUDY**2*PI11**2) 

GO  TO  440 
430  VEL  TRiJEsVELlC 
440  CONTINUE 

* 

ARHO  AND  uBAR  CALCULATIONS  reuuire  table  lookup  uf  zeta 
GEOAlTsREaalT I TuDE/(REaALT ITUDE) 

RHOAL  T=GEUALT 

call  P0INTF(RH0ALT,41,PLACt,FRACTI0N,ZETARG) 

ZETASZE TABL (PL  ACE) t(ZETAbL( PLACE  a 1 )-ZETABL( PLACE )) ‘FRACTION 
RHOsHHO0*txP(-ZETA*GEUALT/10000) 

0 B A K = . S*RHU» vEL TRUE* *2 

•EQUIVALENT  AND  CALIBRATED  AIRSPEEDS  ARE  CALCULATED  IN  KNOTS 
*KEAS  IS  E GUI  VALENT  AIRSPEED,  KCAS  IS  CALIBRATED  AIRSPEED 
KEAS  = VEL I rUE*SORT (RhU/RHUO)aknOTCDN 
* THE  CALCULATION  OF  RACh  REQUIRES  A TABLE  LOOKUP  OF  SPEED  OF 
‘SOUND  BASED  UN  ALTITUDE 
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CALL  PUINTF  ( RHOALT , 5 1 , PL AC t , FR AC  II  ON , SND ARl. 1 

SOUND aSNDT AbL ( PLACE )♦ (SND TAbU PL ACEtl  ) -SNU  T AbL  ( PL  AC  t ) ) *F  P AC  T I UN 
MACH  * VEL 1 RUt /SOUND 

• REAL  FUNCTION  SPEED  IS  USED  TO  DERIVE  kCAS  FROM  KEAS,  AND 

• IS  PPESEN TED  IN  THE  FUNCTION  DEFINITION  SECTION 

P1T1  a SPEEO(MACH) 

P I T <?  a SPEED! .00151 1*KCASJ 
KC AS  * PI  11 *KEAS/PI W 


* analuu  uutput  and  digital  input  are  STARTED  here. 

* 

INLINE 


Lrt 

4,6Z0003FFFF 

Lt  A 

1 , J615 

S T MW 

1 ,<?56 

CD 

0»4Z8000 

END  I 


* IF  DYNAMIC  CHtCKS  ARE  DESIRED  THE  APPROPRIATE  SURFACES  ARfc  DEFLECTED 

* AND  FCS  PROCESSING  IS  BYPASSED. 

* 

IF ( .NOT .DtPULSE ) GO  TO  443 
DELaDtLIC  + PULSEMb.O 
OERaDERIC+PULSE* 16.0 
GO  TO  446 

443  IF ( .NOT .DAPUUSE ) GO  TO  444 
DEL  * DEL IC+PULSE* 16. 0 
0ER=DER1C-PULSE* 16.0 

GO  TO  446 

444  IF (.NOT .ORPULSE)  GO  TO  445 
URUDsDRlC*RuLSfc*20.0 

DERaOERlC 
DElsDELIC 
GO  TO  4 4 6 

445  CONTINUE 

* 

•EXECUTE  SHUTTLE  FlIGhI  CONTROL  SYSTEM 
CALL  SMTLFCS 

* 

IF C .NOT. (FCSDE.OR.FCSDA) ) GO  TO  446 

DELSDELIC 

DERsUERIC 

448  I F ( FC  SDR ) DRUD  = DRIC 

446  IF(FCSDF)  FLAPSD6FIC 
IF (FCSSb)  bRAwsDSbIC 
ELEVs(OEL*DER)/<> 

DAILs(DEL"DERI/2 

* 

* THIS  SECTION  BUILDS  THE  DATA  HALFWORD  TABLES  TO  BE  OUTPUT  THRU 

• THE  U-A  CONVERTER. 

JJ*1 

DA  T ASUDOT 
INLINE 

bL  1650 

E NO  I 
JJai? 


oa  r ashoot 

INLINE 

6L  )950 

ENOI 

i J J*3 

DA  1 AswDO  T 
INLINE 

BL  )950 

ENDI 
JJ  = « 

DA  T A suBUD i 
INLINE 

BL  mu 

ENDI 
JJsb 

DAT AsPBOOr 
INLINE 

bl  mu 

END  I 

J J36 

DAT ASKBUOY 
INLINE 

bl  mo 

ENDI 
JJ  = 7 
DA ! A s A Z 
INLINE 

BL  mo 

ENDI 
J J *B 
DATAsAV 
INLINE 

BL  mo 

END! 

JJs9 
DATAsAX 
INLINE 

bl  mo 

ENDI 
JJslO 

DAT  A* ALPHA 
INLINE 

bl  mo 

ENDI 
JJa  l l 
PATmsBEIA 
INLINE 

bl  mo 

ENDI 
JJslB 
DA  f ASUBAK 
INLINE 

BL  mo 

t nd  I 

J J s 1 3 

DA  r AxVELTWUE 
INLINE 

hi  mo 

ENDI 
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JJ314 

OATAsMACH 

INLINE 

BL  )950 

END I 
JJ  = 15 
DAI As^HI 
INLINfc 

bL  )9b0 

fc  NO  I 
JJ316 
OAIasHSI 
INLlNt 

BL  ) 950 

fc  NO  I 
JJ=1  7 

DATAsTMfcT A 
INLINfc 

BL  )950 

E NO  1 
JJsift 

OA  f ArALl I TUDfc 
INLINfc 

BL  )950 

fcNOI 
JJslV 
OATAaKt AS 
INLINE 

BL  )950 

ENOI 
JJs20 

OA  T AsVfcLT&Ufc 
INLINfc 

BL  ) 950 

ENOI 
JJ*2l 

DA  r A3HOUT 
INLINfc 

BL  )950 

ENOI 
JJ=22 

OATAswHOOT 

INLINfc 

BL  )950 

fcNOI 
JJ  = 2J 

OA  T AsCHUSSBNO 
INLlNt 

BL  )950 

ENOI 
JJ*2« 

OA  r A = OOrtNKN(i 
INLINfc 

BL  ) 950 

ENOI 
JJ325 
OATAsLAT 
INLINE. 


BL 


) 950 
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ENOI 

OAT  A*LNG 
INLINE 

BL  ) 950 

ENOI 
JJ*£7 
DATa*0ER 
INLINE 

BL  ) 950 

ENOI 
JJa*S 
DAT A*ELEV 
INLINE 

BL  1950 

ENOI 
JJ  = <?9 

uata*dail 

INLINE 

bL  ) 950 

ENOI 
JJs3u 

DATAsALT I TUOE 
INLINE 

bL  ) 950 

END! 

JJ*31 
OA  T AsDEL 
INLINt 

BL  ) 95 0 

t NO  I 
JJs3<! 

DAT  AsORUD 
INLINE 

bL  ) 950 

ENOI 
JJa3i 
OA  f AsbRAK 
INLINE 

BL  ) 950 

ENOI 
JJ*3« 

OAT AsKLAP 
INLINE 

HL  ) 95  0 

ENOI 
J J * 3 5 

OAT A=OZCM0 
INLINE 

BL  ) 950 

ENOI 


* 'READY  ' IS  USED  TU  IN  OILATE  THAT  A NEW  LOOP  IS  TO  START. 

* IT  IS  SET  TO  .TRUE,  BY  THE  INTERVAL  TIMER  INTERRUPT  HANDLER. 

READY*. FALSE . 

* 

•TEST  FOR  DESIRED  EASE 

•EASE  IS  REQUESTED  BY  SETTING  SENSE  SWITCH  0 * .TRUE. 
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500  CALL  SSwTChC  12,EASER) 

IF(.NOT.EASEH)  GU  10  550 
505  CONTINUE 
INLINE 

0 1 43  D I SAULE  INTERVAL  UMt« 

E NO  I 

* C ANGLES  ARE  FOR  COMMUNICATIONS  TU  AND  FROM  THE  OPERATOR  VIA  EASE 

* THEY  ARE  EXPRESSED  IN  DEGREES. 

CTHET  AsThE I A *57. 295 77 951 
CPHI=PHI*57 .29577951 
CPSI*PSI*57. 29577951 
CALPHA*ALPHA*57. 295 77951 
CPSlHsPSIH*57. 29577951 
CPSIR*PSIR*57. 29577951 
CMU=MU*57. 29577951 
CLAT*LAT*57. 29577951 
CLNG«LNG*57. 29577951 
CALL  mjluad  ('EASEFCS  Ml) 

CALL  EASE(OATADBL,DBLSIZE,DATAnRDR, RSI ZE, DAT AwRuI , I SIZE, 

* DATAHAF,HAFSIZE,DATABY1 , BYT  SIZE, LI  ST, SIZE) 

GO  TO  10 

* 

*#AIT  FOR  INTERVAL  TIMER  INTERRUPT 

550  TIMEaT iMEtOELTA 
IFCTIME.GT. RUNTIME)  GO  TU  505 

551  CONIINUE 

IF (READY ) GO  TO  47 
GO  TO  551 

* 

* THIS  IS  THE  ANALOG  INPUT  IOCB 

INLINE 

) 9 1 0 OAT  Aw  bZC 0010020 

ACH  ) 9 1 1 

ACM  0 A T A 1 N 

RES  IN 

) 9 1 1 RES  IN 

* 

* THIS  IS  THE  LINKED  IOCB  FOR  BUTH  ANALOG  INPUT  AND  DIGITAL  UUTPUf 

) 9 1 5 DAT  Aw  BZ4200002C 

ACh  ) 9 1 6 

ACH  DATAOUT 

RES  1« 

) 920  DATA*  BZE0000001 

ACH  ) 92 l 

ACH  OISCRTIN 

RES  lw 

) 9 1 6 RES  lw 

) 92  1 D A I A *v  BZ40000000 

END! 

* 

INLINE 

) 950  ST*  0 , L 1 NK 

END! 

IF CSCALOUT IJ J ) ,L T .0)  GU  ID  954 
IF(DATA.LE.SCALOUT(JJ))  GO  TO  951 
DATAOUT  (JJ)s527t>7 
GO  10  953 

951  IF COATA.GE.-SCALOOT(JJ))  GO  TO  952 


DA  f AOOT  (JJ  )*"327b7 

GO  TO  953 

952 

DAT AOOT (JJ)*DAT A/SC ALUOT(JJ) *32767 

INLINE 

) 953 

L»v 

3,  JJ 

SLA 

3,1 

LH 

7,DATA0UT-2,3 

SHL 

7, a 

STH 

7 , DA  T AuU  T-2 , 3 

L« 

0 , L I NK 

r ns* 

0 

END  1 

95« 

IF (DAT A.GE 

.SCALOUT ( JJ) ) GO  TO  955 

0 A T AOOT  (JJ  ) *32  767 

GO  10  953 

955 

IF (DATA.LE 

,-SCALOOT (JJ) ) GO  TO  952 

DATAOUT ( JJ)s-32767 

GO  TO  953 

INLINE 

) 960 

S T ft 

0 , L I NK 

Lft 

3,  JJ 

SLA 

3,1 

LH 

7 , D A T A I N*2 , 3 

sll 

7,1 

STH 

7 , D A T A I N«2 , 3 

ENDI 

DATAaDATAINT JJ)*SCALIN( JJ1/32767 

INLlNt 

Lift 

0 » L I NK 

THSft 

0 

* 

ENDI 

* 

INLINE 

) 1000 

AC  ft 

) 1 0 1 0 

) 1010 

HE  S 

1 ft 

S T ft 

7,SAVFL 

Lb 

7 , 2ZFF 

STo 

7, READY 

Lft 

7 , S A VF  L 

SHI 

* ) 1010 

* 

ENDI 

* 

INLINE 

)1  100 

D 1 

93  DISABLE  INTERVAL  TIMEH 

Lft 

7 , handler 

S T ft 

7,428 

ENDI 

STOP 

INLINE 

11105 

01 

43  DISABLE  INTERVAL 

LEA 

6,  ) 1 105 

LW 

7,8ZFFFFFFEC 

calm 

84 

ENDI 

GO  TO  10 

END 
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APPENDIX  B 


EASE  PROGRAM  LISTING 


1 

i 


# 

* t A S t 


* 

* 


* 


********************************************************************** 

SUBROUT  I Nt  EASE  (OAT ADSL,  OBLSIZE,  DA T A wROR , R S I ZE » D A T A wRO I , I S I Zfc  * 

1 DATAHAF ,HAFSIZE, DAT ABY  T, BY TSIZE, LIST, SIZE  ) 

* THE  PURPOSE  OF  THE  EASE  PROGRAM  IS  TO  ALLOW  INPUT  AND  OUTPUT  UF 

* OATA  VALUES,  OATA  MAY  BE  REFERENCED  BY  EITHER  VARIABLE  NAME  OR 

* numbered  memory  LOCATION. 

* 

* DATA  IS  COMMUNICATED  IN  THE  FOLLOWING  FASHION* 

* 

• 

* NAMED  VARIABLES  ARE  ADDRESSED  BY  NAME,  AnD  MEMORY  LOCATIONS  ARE 

* ADDRESSED  BY  A FULLWORO  ADORESS. 

* 

* TO  EXAMINE  a DATA  VALUE  a QUESTION  MARK  SHOULD  FOLLOW  THE  LUCATIUN 

* DESCRIPTION. 

* 

* TO  CHANGE  THE  VALUE  OF  THE  SPECIFIED  LOCATiuN  AN  EQUAL  SION  SHOULD 

* FOLLOW  THE  LUCATIUN  DESCRIPTION. 


* INQUIRIES  OF  UR  INPUTS  TO  MEMORY  LOCATIONS  THAT  ARE  INTEGER  VALUES 

* SHOULD  BE  PRECEDED  BY  AN  ASTERISK. 

* INQUIRIES  OF  OR  INPUTS  TO  MEMORY  LUCATIOnS  THAT  ARE  REAL  VALUES  may 

* BE  LEFT  UNPREFIXtU. 

* 

* IMBEOUtU  BLANKS  ARt  IGNORED  IN  ALL  CASES. 

* 

* TO  RETURN  TO  ThE  MAIN  PROGRAM  An  EXCLAMATION  MARK  SHOULD  Bt  TYPED. 

* 

* HEX  ADDRESSES  BEGINNING  wITh  A LETTER  (C.F.  A/2C)  MUS  BE  PHtCEEOfcU 

* BY  A ZERO  (C.F.  0A72C) 


* ThE  FOLLOWING  Input  WILL  CAUSE  ALFA  TU  bE  SET  TU  30  UEGREES- 

* ALFA»30 


* 

THE  FOLLOWING 

WILL 

CAUSE 

THE 

PROGRAM 

TO 

uu 

IPUl 

ZETA- 

* 

Zfc  T A ? 

* 

ThE  FOLLOWING 

WILL 

CAUSE 

The 

PkOGRAM 

TU 

uu 

IPUT 

ARE  REAL  V 

ALUE 

F ROM 

* 

LOCATIUN  12345 

* 

12345? 

* 

THE  FOLLOWING 

WILL 

INPUT 

THE 

INTEGER 

VALUE 

UF 

546  TO  LOCA 

TION 

12345 

* M2345»456 

* The  DATA  STORAGE  ARRAYS  ARE  SET  UP  AS  FOLLOWS 

* * 

* * 

* OATADBLOMSKEAL  doubleword  array  * 
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* 

* 


kskmaxDHL 


UATArt«l)H(K)=KhAL  WORD  AkRAY 


KsKMAXKt AL 


DATAwRUi (K)slNTEGER  WORD  ARRAY 


KsKMAXwRD 


DAI  AHAF (K )=lNTEGtR  HALFWORD  ARRAY 


KsKMAXHAF 


OATABY1  (K)=InTEGER  BYTE  ARRAY 


KsSIZE 


THE  LIST  SYMbUL  ARRAY  CONTAINS  A DOUBLEWORD  CHARACTER  NAME  FOR  EACH 
VARIABLE  IN  the  DATA(K)  ARRAY.  The  VALUE  OF  K in  LIST ( K ) CORRESPOND 
TO  THE  VALUE  OF  K IN  THE  DATA(K)  ARRAYS 

THE  FOLLOWING  VARIABLES  MUST  BE  DIMENSIONED  AS  INDICATED  TO  INSUkt 
PROPER  OPERATION  UF  THE  PROGRAM  IN  EACH  CASE 
OATABYT (VALUE  UF  K M A XHAF , S I ZE ) 

DATAHAF (VALUE  OF  KMAXWRD,  VALUE  OF  KMAXHAF ) 

DATAwRDI (VALUE  OF  KMAXKEAL, VALUE  OF  KMAXWRD) 

DATAWRDR( value  of  kmaxdbl, value  OF  KMAXREAL) 

DATADBLdr  VALUE  OF  KMAXDBL) 

LISHTOTAL  NUMBER  OF  VARIABLES) 


ALSO,  THE  FuLLOwING  PARAMETERS  MUST  BE  INITIALIZED  IN  A DATA 
STATEMENT 

size=iotal  NUMBER  OF  variables 

KMAXDBLaVALUE  OF  K AT  END  OF  DOUBLEwORD  ARRAY 
KMAXREAL*V aloe  OF  K at  END  OF  REAL  WORD  ARRAY 
KMAXwRD=VALUE  OF  K AT  END  OF  INTEGER  WORD  ARRAY 
KMAXHAFsVALUE  OF  K AT  END  UF  HALFWORD  ARRAY 


fall  intgr, letter, decimal 

Char  WURUIBI ) , TEMP(B) , NUMB (5) , 0UT8UF(2b) ,NAMEIN(2S) 

ChAk  MSSG1  ( IB) ,MSSG£(18) ,MSSG3(3«) ,MSSG4 (33) ,MSSG5(3B) ,MSSG6( 31  ) 
CHAR  MSSG7  (37)  ,MSSG8(i6)  ,mSSGR(<?9)  *MSSG10(«?^) 

INTEGERS  dblsize,rsize, isize,hafsize,bytsize,size 

I N T EGER* 1 DATABYT (BYTSIZE) 
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INTEGER*2  DATAHAF IHAFSIZE) 

IMEgEH*4  OATAwkDl ( lSIZfc) 
lNTEGtR*4  AUDH  » ttOX  I 
INTEGER* 8 LIST (SIZE) # INbUF I , SYMBOL 
neal*4  da r awrdrirsizei # boxk 

KEAL*tt  DAI  AUbL(DBLSIZE)  > I NBUF  R 

EQUIVALENCE  C SYMBOL t TEMP ( II 1 , C TEMPI 1 ) , NUMB! 1 J ) , I BOXR , BOX l ) 
DATA  MSSG1/ ' WELCOME  ID  EASE'/ 

DATA  MSSG2/ ' WRONG  INPUT  FORMAT'/ 

DATA  MSSG3/ ' VARIABLE  NAME  EXCEtDS  6 CHAkACTERS'/ 

DATA  MSSGR/'NO  SUCH  SYMBOL  IS  DEFINED  TO  EASE'/ 

DATA  MSSGB/ ' MEMORY  ADDRESS  MUST  BE  A HEX  NUMBER'/ 

DATA  MSSG6/ 'MEMURY  ADDRESS  EXCEEDS  5 DIGITS'/ 

DATA  MSSG7/'MEMORY  ADDRESS  IS  OUTSIDE  CURE  LIMITS'/ 

DATA  MSSG8/ 'MEMORY  DATA  MUST  BE  A OECIMAL  NUMBER'/ 

DATA  MSSG9/ 'MfcMURY  DATA  EXCEEDS  25  DIGITS'/ 

DATA  MSSGIO/'ENO  OF  EASE  PkUCESSlNG'/ 

KMAXDBL*DBLSIZE 

KMAXREAL=DBLSIZE*RSIZe 

KMAXwROeDBLSIZE+RSIZEtlSIZE 

KMAXHAF=D«LSIZE*RSIZE+ISIZE*nAFSIZE 

CALL  CARRIAGE 

CALL  ms IELEaIMSSGI, 151 

CALL  carriage 

10  CONTINUE 

LETTER=INIGR=.FALSE. 
call  CARRIAGE 
call  ms TELER(RORU,80J 
L*2 

TEST  FOR  BLANKS  IN  INPUT  LINE 
15  IF  IwORDTL J ,NE. ' *1  GO  TO  20 
L*L*  1 
GO  TO  15 

20  IF (wORO(L) .fcu. ' 1 ' ) GO  TO  100 

IF ( (wuRD(L) .LT . ' A ' ) .UR. (wORU(L) .GT . ' Z' 1 ) GO  Tu  22 
GO  TO  200 

22  IF  (WORU(L).NE.  '*'  ) GO  TU  24 
INTGR*. TRUE. 

23  L*L*  1 

IF (wORD(L) .NE. ' ')  GO  TO  24 
GO  TO  23 

24  IF(  (wuRDCL) ,LT . ' 0 ' ) .OR. (hORD(L) .GT. '9' 1 ) GO  TO  30 
GO  TO  300 

40  CONTINUE 

IF ( WORD ( L ) ,EU. '? ' ) GO  TO  400 

IF(wOkO(L) .EQ. 1 GO  TO  500 
GO  TO  30 


THE  FOLLOWING  SECTION  CONTAINS  ERROR  MESSAGES 


30  call  carriage 

CALL  ms  I ELEw (MSSG2» IB) 
GO  TO  10 

50  CALL  CARRIAGt 

CALL  msTELEw(msSG3,34) 
GO  TO  10 


60  CALL  CARRIAGE 

CALL  Mj TELfc*(MSSG«, 33) 

GO  IU  10 

* 

70  CALL  CARRIAGE 

CALL  MjTELEw(MSSG5, 35) 

GO  TO  10 

* 

80  CALL  CARRIAGE 

CALL  M: TELEwIMSSGb. 31 ) 

GO  TO  10 

* 

90  CALL  CARRIAGE 

CALL  *:TELE*lMSSG7,i7) 

GO  TO  10 

* 

55  CALL  CARRIAGE 

CALL  M:TELE«(MSSG8»36) 

GO  TO  10 

* 

65  CALL  CARRIAGE 

CALL  Ml  TELEw(MSSG9,29) 

GO  TO  10 

* 

* 

* 

* NAMED  VARIABLE  PROCESSING 
200  CONTINOE 

* LETTER  SIGNIFIES  TO  THE  I/O  SECTION  THAT  THE  EFFECTIVE  VARIABLE 

* ADDRESS  IS  DEFINED  BY  LOCATION  IN  Th£  VARIABLE  TABLE. 

LE I TERa, TRUE. 

* initialize  TEmP  TO  blanks 

TEMRs'  ' 

J*1 

210  TEMPI J)aWORU(L) 

220  L*Lt 1 

I F ( WORD ( L ) . NE . ' ’)  GO  TO  250 

GO  TO  220 

230  IF l (WOHD(L) .EU. ’?•  ) . OR . I NORD ( L) . EO . ' * • ) ) GO  TO  250 
J*J-Fl 

1F(J.EU.9)  GO  Tu  50 
GO  TO  210 

* TEST  FOR  SYMBOL  IN  SYMBOL  LIST  AND  FIND  RELATIVE  LOCATION  IN  DATA 

* ARRAYS 
250  K s 1 

255  IF (SYMBOL. EU.LIST(K) J GO  TO  260 

K3RM 

IF(K.GF.SIZE)  GO  TO  60 
GO  TO  255 

260  IF (K .GI .KMAXREALJ  INTGK*. TRUE, 

* RETORN  TO  MAIN 

GO  TO  UO 

* 

* 

* 

* numbered  address  processing 

300  CONTINUE 

* FILL  NUMo  rtlTrt  BLANKS 

NUMB*'  • 
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J*1 

310  NUMB(J)*wOKD(L) 

320  L*L  + 1 

IF  (wORD(L)  .NE.  ' ')  GO  To  330 

GO  TO  320 

330  IF ( («ORD(L) ,EU. ' ) . UR , ( WORD ( L ) .EG . ' = ' ) ) GO  TU  350 
1F( (WONU(L) .GE. ' O' ) . ANU. (rtURU(L) ,LE. ’9* ) } GO  TU  340 
IF ( (rtORD(L) , GE . ' A* ) .AND. (nOROCU) . Lt . * F * ) ) GO  TO  340 
GO  TO  70 
3«0  J*J*1 

IFU.EU.6)  GO  TO  BO 
GO  TU  310 

* NUMB  MUST  Bt  WIGHT  JUSTIFIED  FOR  CONVERSION 
350  IF(J.EU.S)  GO  TO  357 

M0VEs5-J 

DO  355  N»  1 » MOV E 

NUMB(5)»NUMB(4) 

NUMB(4)rNUMB(3) 

NUMB(3)»NUMBC2) 

NUMB(2)sNUMB(l  ) 

355  NUMB ( 1 ) s ' ' 

* CONVERT  ASCII  INPUT  TO  HEX  AUURESS 
357  DECOOE(5,3bO,NUMB) AUDK 

360  FORMA  T ( Z5 ) 

* TEST  FUR  EXISTENT  MEMORY  ADDRESS 

IF(ADOH.GI .5Z27FFC)  GO  TO  90 

* RETURN  TO  MAIN 

GO  TO  40 


* CRT  OUTPUT  RROCESSING 

* TEST  TU  DETERMINE  IF  NUMBERED  LOCATION  UW  NAMED  VARIABLE  UUTPUT  IS 

* DESIRED 
400  CONTINUE 

IF (LETTER)  GO  TU  450 

* RETIEVE  VALUE  FROM  EFFECTIVE  ADDRESS 

INLINE 

Lw  7 , * ADDW 

STw  7 » BOXR 

END  I 

* TEST  FOR  INTEGER  OR  WEAL  CONVERSION 

ifuntgrj  go  to  4io 
ENCODE (25, 405, OoTBUF) BOXR 
405  FORMA  I (F25. 10) 

GO  TU  490 

410  ENC00E(25,415,UUTBUF)B0XI 
415  FORMAT (125) 

GO  TU  490 

450  IFUNTGR)  GO  TO  470 

* TEST  TU  DETERMINE  WHICH  ARRAY  TYPE  IS  TO  BE  ACCESSED.  AND  PROCESS 

* ACCORDINGLY 

IF(K.bT.KMAXDBL)  GO  10  4 b 0 
ENCODE (25.405. OUT BoF )D AT  ADBL(M 
GO  To  49u 

46  0 ENCODE (25, 405. UUTBUF) DAT AwRDR(K-KMAXDBL  ) 

GO  TU  490 

470  IF (K.GT .KMAXHAF ) GO  TO  4P5 
IF (K.GT .*max«Rd)  GO  TU  460 
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ENCUDt(25,4l5,UUlBUP)0ATArtRUl(K-KMAxKEAL) 
GO  TO  490 

4 8 0 ENC00t(25,4l5,U0t8UF)DAjAHAF  (K-kmaxwRD  ) 

GO  10  440 

4 85  ENC0DE(25,415*0UIBUF ) DAT  Art  1 T (K-KMAXHAF ) 
OUTPUT  ANSWER 
490  CALL  LArRIAgE 

CALL  M:TtLEw(0UTBUF,25) 

WE  T URN  TO  MAIN 
GO  TO  10 


CWT  INPUT  PROCESSING 

FINISH  INPOI  OF  DESIWtU  VALUE  ANO  STOWE 
500  CONTINUE 

DECIMALS. PALSE. 

J*  1 

NAME  I NS  • ' 

510  L = L+1 

TEST  FUW  BLANKS 

IF (wORDlL ) .EU. ' ')  GO  TO  510 

TEST  FUR  CAnWIAGE  WETOWN 

IF  (wOwO(L  ) .tU.2Z0D)  GO  TU  520 
TEST  FOR  vALIU  NUMBER 

IF ( (kQWOCL ) .GE , ' 0 ’ ) . ANO. (wORD(L ) ,LE . ' 9 ' ) ) GO  TO  515 
IP (wURD(L).EO. I GO  TO  515 
IF  (wORD(L) ,N£. • . ' J GO  TO  55 
DECIMALS. TRUE. 

515  NAMtlN(J)SrtUHU(L) 

■J  s J ♦ 1 

IF  U.tU.^o)  GO  TO  85 
GO  TO  510 

NAME  I N NOw  CONTAINS  THE  VALUE  TO  BE  STORED  IN  ASCII 
TEST  FOR  INTEGER  OR  REAL  CONVERSION 
520  IF (INTGR)  GO  TO  525 

IF C ,NU1  .DEC IMAL  ) NAMEINI JI*'  . • 

OE CODE  (25, 4 0 5,  NAME  INHNBUFR 
GO  TO  530 

525  IF (J, EG. 251  GO  TO  529 

NAMEIN  MUST  Bt  RIGHT  JUSTIPItO  FOR  CONVERSION 
MOVE326-J 
00  527  1 = 1 , MO V t 
DO  52b  Nil, 24 

526  NAMEIN(26-n)sNAMEIn(25-n) 

527  NAME JN( 1 )='  • 

529  DECuut(25,4l5,NAMElN)IN«UFl 

INHUPR  CONTAINS  THE  REAL  DATA  VALUE,  ANO  InhUFI  CONTAINS  THE  INTtGER 
DATA  VALUt 

TEST  WHETHER  DATA  IS  FOR  NAMED  VARIABLE  OR  NUMBERED  MEMUHY  LOCATIUN 

530  IF(LtTTER)  GU  TO  550 
IF ( I N I GR I GO  TO  540 

STORE  RtAL  VALUE  IN  EFFECTIVE  ADDRESS 
INLINE 

L*  7,1 NbOF  R 

STw  7 , * AuOR 

t NO  1 

GO  Tu  599 

STORE  INTtGER  VALUE  IN  EPPtCUVt  ADDRESS 


7,  INbUF  If  lw 
7 » * AQOR 


540  CONTINUE 
INLINE. 

LW 
ST  w 
ENl)  I 

GO  TO  599 

* CHOOSE  DATA  I ABLE  IN  WHICH  VARIABLE  IS  TO  BE  STORED 
550  IF (K.GT .KMAXHAF ) GO  TO  560 

IF  (K.GT.KMAXwRO)  GO  TO  562 
IF (K.GT .UMAXREAL)  GO  TO  569 
IF (K.GT.KMAXDBL)  GO  TO  566 
DATADBL(K)alNBUFR 
GO  10  570 

560  OATAbTT(K-KMAXHAF)*lNBUF I 
GO  TO  570 

562  OATAHAF (K-KMAXwRD) »INBUF I 
GO  TO  570 

564  DATAwRDI (K-KMAXHEALJ*INbUFl 
GO  TU  570 

566  DATAwKORCK-KMAXDBLjsINBuFR 

* THE  FOLLOWING  SECTION  TEST  FOR  VALUES  OF  K THAT  HEUUlKt  ADDITIONAL 

* CALCULATIONS  FOR  EACH  CHANGE  IN  VALUE 
570  CONTINUE 

* 

* 

* insert  ADDED  CALCULATIONS  here  ****************** 
************** 

*************** 

* RETURN  TO  MAIN 
599  GO  TU  10 

* 

* 

100  CALL  CARRIAGE 

CALL  m:TELEw(mSSG10.22) 

RETURN 

* 

* END  OF  EASE  PROGRAM 
ENU 
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FUNCTION  GENERATION  PACKAGE 


* 

* 

* 

* 

* 


PROGRAM  GATASTUR 
REAL  DATAUOOO) 

INTEGERS  NAME  , NAMELIST  (**5) 

NOTE  THAT  OOTDATA  MUST  HE  OIMENS10NE0  TO  'SIZE'/  OR  THE  TOTAL 

NUMBER  OF  ftORDS  IN  ALL  DATA  ARRAYS. 

REAL***  OUTuATA  (S2**3) 

INTEGER***  LENGTH,  NUMBER,  1)1  FF,  ADDRESS,  SIZE,  PIT,  ST  AT 

REAL***  CC0K63) 

REAL***  OCCET  (378) 

REAL***  CCOBFlT(bi) 

REAL***  CCDbFZT  (63) 

REAL***  CCOSB  IT  ( 63 ) 

REAL***  CCOSB2T  (631 

REAL***  CLB  T AbL  1252  ) 

REAL***  ClBOE  1 ( ( 63 ) 

REAL***  CLBDE2  T (63) 

HEAL***  CLbOSB  1 T ( 63  ) 

REAL***  CLBOSB2T  (63) 

REAL***  CL60SB3T(63) 

REAL***  CLOAT  (378) 

REAL***  CLOR  T (252) 

REAL***  CLP T ( 63 ) 

REAL***  CLRT(b3) 

WEAL***  C MO  T ( 6 3 ) 

HEAl***  OCMET(37B) 

REAL***  CM08F  1 T ( 63 ) 

REAL***  CMDHK«?T(&3) 

heal***  CMDSB1T  (63) 

REAL***  CM0SB2TC63) 

REAL***  CMUM252) 

REAL***  CNB  I ABL  (252 ) 

REAL***  CNBOE  1 T ( 63 ) 

REAL***  CNBDE2T  (63) 

REAL***  CNBDSB1  T (63) 

REAL***  CNB0SB2T  (63) 

REAL***  CNBDS83T  (63) 

REAL***  CNDAT  (37B) 

REAL***  CNOw  T ( 252  ) 

REAL***  CN0RBTC63) 

REAL***  CNPT  (63) 

REAL***  CNR T C 63 ) 

REAL***  CnOT  (63) 

REAL***  CN0EH63) 

REAL***  CNDBMTC63) 

REAL***  CnOBF2T(63) 

REAL***  CN0S81  T ( 6i  ) 

REAL***  CN0Sd2T(63) 

REAL***  C YB T (63) 

REAL***  CYOAT  (378) 

HEAL***  CYOnTIbi) 

REAL***  C YORSB 1 T ( 7 ) 

REAL***  CY0HSb2T(7) 
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CHAR  BUFFEH(81  ) , COMMtN T 3 ( 65 ) ,MSSG(80) , TEST 
CHAR  MSSGl (45) ,MSSG2(44) ,MSSG3(64) 
EQUIVALENCE  (BUFFER (2 ) , MSSC(l),  TEST) 
EQUIVALENCE 

* <OUTDATA(l),CCOT(l)) 

*, (0UTDATA(64),DCCET(1)  ) 

*, (OUT DAT  A (442) ,CCD8F 1 T ( 1 ) ) 

*, (UUTDATA(50S) ,CCO0F2T (1 ) ) 

*,  (OUTOATA (568)  ,CCOSBl  T (1)  ) 

*, (OUTOATA (631 ) ,CCDSB2T (1 ) ) 

*, ( OUTOAT  A ( 694  ),  CL8  TABLU  )) 

*,  (OUTOATA (946),CLB0E1T(1)) 

*, (0UTDATA(1009) ,CL80E2T(1)  ) 

*, (OUTOATA( 1072) ,CL60SB1 T (1 ) ) 

*, (OUTOATA( 1 135) ,CLB0SB2T (1 ) ) 

*, (OUTOAT A ( 1 198),CLBDSB3T( 1 ) ) 

*, (OUTQATA( 1261 ) ,CLDAT (1)  ) 

*,  (OUTl)ATA(  lt>39)  ,CLl)KT  ( 1)  ) 

* , (UUTOAT  A ( 1891  ) , CLP  T (D) 

* t (OUTOATA  ( 1954) , CLRT  (D) 
*,(OUTl)ATA(2017),CHUT  ( 1 ) ) 

*, (OUT  DATA (2080), OC«ET(l)) 

*,  (0UT0ATA(2458) ,CMDBF1  T ( 1 ) ) 

*, (OUTOAT A (2521 ) ,CMOBF2T ( 1 ) ) 

*, ( OUTOAT A ( 2584 ) ,CM0SB1 T ( 1 ) ) 

*, (UUTDATA(264  7),CMDSB2T(1 ) ) 

* » (OUTOATA(2710) , CMU  T (1 ) ) 

*, (OUTOATA l 2962) ,CNBTABL(  1 ) ) 

*, (0U10ATA(3214).CNB0E1T(1  ) ) 

*,  (OUTOAT  A(  3277  ) , CN8DE2  T (D) 

*, (OUTOAT A (3340 ),CNB0SB11 (1)) 

*, (0UTUATA(3403) ,CNBOSB2T ( 1 ) ) 

*, (OUTOAT A (3466) ,CnBDS83T ( 1 ) ) 

*, ( OUTOAT A ( 3529) ,CNOAT ( 1 ) ) 

*, (OUTOATA (3907 )»CNQRT(1)) 

*, (OUTOAT A (4159), CNORBT ( 1 ) ) 

*, (OUTOAT A ( 4222), CNPT(l)) 

*, (OUTOAT A ( 4285) , CNR T ( 1 ) ) 

*, (OUTOAT A( 4348 ) ,CNOT( 1 ) ) 

*, (OUTOAT A (44  1 1 ) ,CN0ET ( 1 ) ) 

*, (0UT0ATA(44  74),CN0BF1T(1  )) 

*, (OUT  OAT  A (453 7) ,CN08F2T ( 1) ) 

*, (OUTOAT A14600) ,CNDSB1T (1)) 

*, (OUT DAT A( 466  3), CN0SB2T(1 ) ) 

*, (0UTDATA(4726),CYBT(1)) 

*, (OUTOAT A (4789) ,CY0AT(  1 ) ) 

*, (OUT  DATA (5167), CYORT(l)) 

*, (OUTOAT A (5230), CVORSBlT(l)) 

*, (OUTOAT A( 5237 ) , C Y 0RSB2 T ( 1 ) ) 
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. .... 


* 


DATA  M$SG1/'IF  HEADER  CAND  CONNECT  TYNE  C,  IF  NUT  TYPE  A?*/ 

DATA  M$SG2  /'PLEASE  INPUT  CONRECT  HEADER  CARC  INFORMATION'/ 

DATA  MSSG3/ ' NAME  HAS  NOT  FOUND  IN  NAMELIST  PLEASE  ENTER  CONNECT  IN 
1FONMATION. '/ 

DATA  SIZE/5243/ 

DATA  NUMBEN/45/ 


DATA  NAMELIST/ 

» 'CCO' 

* ,'CCE' 

* » ' CC  DBF  1 ' 

* , ' CCDBF  2 ' 

* ,'ccosei' 

* » ' CCDSB2 ' 

* , ' CLB ' 

* , ' CLBDE 1 ' 

* » ' CLB0E2 ' 

* , ' CL bDSB  1 ' 

* , ' CL6DSB2 ' 

* , 1 CLBDS63 ' 

* , 'CLDA' 

* » ' CLDN ' 

* * ' CLP ' 

* #'CLH' 

* ,'CMO' 

* ,'CME' 

* , 'CMDBF 1 ' 

* , ' CMDBF 2 ' 

* , ' CMDSB 1 ' 

* » ' CMDSB2 ' 

* , ' CMG ' 

* , ' CNB ' 

* / ' CNBDE 1 ' 

* » ' CNBDE2 ' 

* , ' CNfeDSB 1 ' 

* ,'CNBDSB2' 

* , ' CNBDSB3 ' 

* , 'CNDA' 

* , 'CNDR ' 

* , ' CNDNB ' 

* ,'CNP' 

* , ' CNN ' 

* ,'CNO' 

* , ' CNDE ' 

* , ' CNDBF  1 ' 

* , 'CNDBF  2 ' 

* ,'CNOSBl' 

* , ' CNDSB2 ' 

* ,'CYB' 

* ,'CYDA' 

* ,'CYOR* 

* ,'CYDRDSB1' 

* , ' C YDNDSB2 ' / 

10  DATA  s 0.0 

READ  (4,20. END* 150)  NAME, 
20  FORMAT  (A8,1X,I6,65A1) 

30  ENCODE  (60,20, MSSG ) NAME, 
IF  (TEST  «EU. ' C ' ) GO  TO  50 


LENGTH,  COMMENTS 
LENGTH, COMMENTS 
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* THIS  SECTION  OUTPUTS  THE  HEADER  CARD  TO  THE  CRT 


CALL  CARRIAGE 

CALL  M:TELE*  (MSSG,80) 

CALL  C ARR I AGE 

CALL  MiTELER  (MSSG1.4S) 

CALL  MJ TEL£R(BUFF£R,80) 

IF  (TEST.EU. ' l • ) GO  TO  150 
IF  (TEST.tu. *C»)  GO  TO  50 
CALL  CARRIAGE 
CALL  MjTELE*  (MSSG2,44) 

40  call  carriage 

CALL  M : TELE*  (MSSG,  80) 

CALL  MjTELE*  (MSSG,  80) 

CALL  CARRIAGE 

CALL  MjTELER  (BUFFER, 81) 

DEC0DE(80,20»MSSG)  NAME , LENGTH,  COMMENTS 
GO  TO  30 

• 

* SEARCH  EASE  TABLE  OF  DERIVATIVE  NAMES  FOR  NAME  AND  RETURN  RUS1IION 

* 

50  J*0 

51  J=Jf 1 

IF  (NAME. EG. NAMELISI (J)  ) GO  TO  55 
IF  (J.LT .NUMBER)  GO  TO  51 
CALL  CARRIAGE 
CALL  MsTELEw  (MSS63,64) 


GO  TO  40 

CONTINUE 

INLINE 

L* 

3 , J 

SLA 

3,2 

L« 

7, )200-4,3 

ST* 

7 , ADDRESS 

E NO  I 

♦ This  SECTION  READS  The  OATA  CARDS  UNTIL  THE  NUMBER  OF  POINTS  RE AO*LENG T 

* 

59  READ  (4,60)  ( DATA(I), 1=1, LENGTH) 

60  FORMAT  (8F10.6) 

* 

DO  160  1=1, LENGTH 
DATA*  x OATA(I) 

INLINE 


L* 

3,1 

SLA 

3,2 

ADM* 

3,  ADDRESS 

L* 

7, DATA* 

ST* 

END  I 

7,-4, 3 

160 

CONTINUE 
GO  TOIO 

150 

call  carriage 

131 


1*1 

PIT*1 15* 

180  CALL  BUFFEROUT(5,l,OUTDATA(I),PIT) 

185  CALL  ST  ATU3 ( 5*  ST  AT ) 

IF(STAT.NE.2)  GO  TO  185 
1*1+1152 

IF(I.GT.SIZE)  GO  TO  186 

IF  (1  + 1152. GT. SIZE)  PIT*SIZE+1-I 

GO  TO  180 

186  CONTINUE 
STOP 
INLINE 

!00  ACm  CCOT 


ACM 

ACM 

ACM 

ACM 

ACM 

ACM 

ACM 

ACM 

ACm 

ACM 

ACM 

ACM 

ACm 

ACM 

ACM 

ACM 

ACM 

ACM 

ACM 

ACM 

ACm 

ACM 

ACM 

ACM 

ACM 

ACM 

ACm 

ACm 

ACM 

ACm 

ACm 

ACM 

ACM 

ACM 

ACM 

ACM 

ACm 

ACm 

ACm 

ACm 

ACm 

ACm 

ACm 

ACm 

ENDI 

END 


CCOT 

OCCET 

CCDBF1 f 

CCDBF2T 

CCOSB1T 

CCDSU2T 

CL8TABL 

CLBDE1T 

CLBDE2T 

CLB0SB1 T 

CLBOSB2T 

CLB0SB3T 

CLDAT 

CLDRT 

CLP  T 

ClRT 

CMOT 

DCMET 

CMDBF1T 

CM0HF2T 

CM0S81T 

CM0SB2T 

CMOT 

CNBTABL 

CNB0E1T 

CNB0E2T 

CNBOSB1T 

CNBDSB2T 

CNB0S63T 

CNOAT 

CNORT 

CNORB I 

CNPT 

CNR  T 

CNOT 

CNDET 

CNOBF1 T 

CNDBF2T 

CNOSB1T 

CN0SB2I 

C VBT 

C YOAT 

C YORT 

C YDRSB1 I 

CY0RSB2T 
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FUNGI  ION  DtNlvtl  (ARRAY,  I,  ARGPN T1  , ARFRAC1  ) 

INTEGEk*2  AkGPNU,! 

DIMENSION  ARH  A Y ( I I 

OERi  VE1*(AWRAY  (ARGPNT  1 ♦ J I » ARR A Y ( ARGPNT  1 ) )*ArFRAC  1 * ARR  A Y ( ARGPNT  l ) 

RETURN 

END 

* 

******** 

* 

FUNC1ION  DERI VE2( ARRAY2, I , J, ARGPNI 1 , AKGPNT2, AHFRAC 1 , ARFRAC2 ) 

D I MENS  I UN  ARRA  Y2 ( i » J ) 

I NTEGER*2  ARGPN T 1 , ARGPNT2 

TEMP3*(ARRAY2( ARGPN  Tl*l , ARgPnT 2+  1 ) - A RR  A Y 2 ( ARGPN  T 1 , ARGPnT 2*  1 ) ) *4RFk 
*aC1+ARRAY2(AKGPnT1,ARGPnT2*1) 

T EMP2* ( ARK A Y2 ( ARGPN  T 1 ♦ 1 , ARGPN I 2 ) • ARK A Y2 ( ARGPNT 1 » ARGPN T 2 ) )*ARFRAC1* 
*AWrAY2(ANGPNT 1 , ARGPN 12} 

OEk I VE2* ( IEMP3-TEMP2)*AkFHAC2+TEMP2 

RE  TURN 

END 

* 

******** 

* 

FUNGI  ION  DERI  VE3(  ARRAY  3,  I , J , A , ARGPN  M , AR&PN T 2,  ArGPN T3,  ARFkAC  1 , 
*ARFRAC2»ARFRAC3) 

OIMENSION  ARRAY3(I,J,k) 

INTEGER* 2 ARGPN FI , ARGPN I 2, ARGPNT  3 

TEMP7=( ARRA Y3( ARGPN T 1 ♦ 1 , ARGPNT  2+1 , ARGPNT 3* 1 ) - ARR A Y 3 ( AR&PM 1 , ARGPN 
*12+1 , ARGPN I 3 + 1J ) *ARFRAG1 +AKRAY3( ARGPNT 1 , ArGPN T 2*1 , ARGPN T 3*1 ) 
TEMPbat ARRA Y3( ARGPNT It  1 , ARGPN T 2 , ARGPNT  3 ♦ 1 ) - ARR A Y 3 ( ARGPN T 1 , ARGPN T2, 
* ARGPNT  3+ 1 ) )*ARFRAC 1 ♦ ArR A Y 3 ( AR&PN Tl , A RGPM 2 , ARGPN T 3* 1 ) 

TEMPS* ( ARRA Y3( ARGPNT 1 + I » ARGPN T2+1, ARGPNT 3) -ARRA Y 3 ( ARGPN I 1 , ARGPNT  2 + 
* 1 , ARGPN  T3))*ARFRAC 1 ♦ ARRA Yi ( ARgPnT 1 , ARGPNT 2 1 1 , ARGPNT  3) 

TEMPAsIA RR AY3( AkGPNllt l,ARGPNT2,ARGPNT3)-ARRAY3l ARGPNT 1,ARGPNT2,AR 
*GPNT3) J*ARFRAC1+ARRAY3( ARGPNT 1 , ARGPNT 2 f AR&PN T3) 

TEMP3*( TEMP7-TEMP6) *ARFRAC2*TEMP6 
TEMP2* ( TEMPS- TEMPa J*ARFRAC2*TEMPa 
DERI vti*( TEMP3-TEMP2J*ARFRAC3*TEMP2 
RETURN 
END 

* 

******** 

* 

FUNCTION  DERIVE«(ARRAY4,I,J,K,L,ARGPNT1, ARGPN T2 , ARGPN T 3 » ARGPNT a , 
*ARFRAGUARFRAG2>ARFRAg3«ARFRAG4J 
DIMENSION  ARRAY4(I, J,K,L) 

INTEGEK*2  ARGPNT 1 , ARGPN T 2, ARGPN T3, ARGPNT 4 

* DERIVE^  IS  THE  VALUE  RERURNEO  TO  THE  MAIN  PROGRAM. 

* ARRAY4  IS  A FQuh  D I MENS  I UAL  DUMMY  AkRAY  CONTAINING  FUNCTION  DATA  U 

* I IS  AN  INTEGER  VALUE  REPRESENTING  THE  SIZE  OF  THE  FIRST  DIMENSION  OF 

* J IS  AN  INTEGER  value  REPRESENTING  1 HE  SUE  UF  ThE  SECON  DIMENSION  OF 

* k IS  AN  INTEGER  value  REPRESENTING  ThE  SIZE  OF  THE  THIRD  DIMENSION  UF 

* L IS  An  INTEGER  VALUE  REPRESENTING  FIh  SIZE  OF  THE  FUORTH  DIMENSION  0 
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1EMP1S=(APPAY4(AHGPN1  1 + 1 , APGPNT  2 + l , APGPN  T 3 + l , APGPN14+1 ) 

* -APPAY4(APGPNI1,  APGPN  T2+  1 , APgPM3+1,APGPN  I4+1))aAPFPAC1 

* tAKHAYUlAKGPNl l , APGPN T 2 ♦ 1 , A PGPNT 3 + 1 , A PGPN 1 4 ♦ 1 ) 
TEMP14s(APPAY4( APGPNf 1 ♦ 1 , ARGPN I 2 , ARGPN T 3 + 1 , APGPNT 4+ 1 ) 

* -APKAY4 (APGPN  f 1 , AKGPN I 2 , APGPN I 3*1 , APGPNT  4+1 ) ) aaPFPAC 1 

a ♦APPAY4(APGPNTi, APGPNT2*  APGPNT 3+1 , APGPNT4+1 ) 

TEMP13*(APPAY4(  APGPNT  1 + 1 , APGPN T2+1  , AHGPNT3,  APGPM  4+ 1 ) 
a -APR  A Y4(  APGPNf  1 , APGPN 12+1 , APGPNT 3,  APGPNT4+1  ) )*APFPACi 

a +APPAY4(AkGPN?1  , APGPN  I 2 ♦ 1 , APGPN  T 3 , AkGPN f 4 ♦ 1 ) 

TEMP12=(  AKPAY4C APGPNT  l + l , APGPNT 2,  APGPNT 3, APGPNT4+1) 

* -APPA Y4(  APGPN  ( l , APGPN T 2, APGPN 13, APGPN T 4+1 ) ) aANFPACI 

a +APPA Y4( APGPN U,APGPNT2» APGPNT 3# APGPNf 4+1  ) 

ItMPl 1=(APPAYP( APGPN T l+l , APGPN T2, AKGPN T 3+1 , APGPN TP) 
a -APPA Y4 ( APGPN  1 1 , APGPN  T2  + 1 , AkGPN  T 3+  1 , ARGPN T 4 1 ) +APFPAC  1 

a *APPAY4( APGPN T 1, APGPNT2+1 , APGPN T 3+ 1 * APGPN T 4 ) 

TEMP10s(APPAY4( APGPN T 1+1, APGPN T2 , APGPN T 3+ 1 *APGPNf4) 

* -AHPAY4(AHGPNT 1 , ARGPN T 2 , A P GPN 1 3 + 1 ,APGPNT4))*APFPACl 

a +AHPAY4  (APGPN  f 1 , A WGPN  I 2 , AKGPN  T 3 + 1 , APUPNT4) 

TEMPRs(APPAY4(APGPNTl+l,APGPN12+l , APGPN T3,APGPn14) 

A -APPAY4CAPGPNI1, APGPN 12+1, APGPN T 3 , APGPN T 4 ) ) a APF  WAC 1 

a ♦ APPAY4(APGPMl,  APGPNT  2+1  ,APGPNT3,ApGPNT4) 

TEMP8*(APPAY4(APGPnT1+1, APGPN T 2 , APGPN I 3 , APGPN  T 4 ) 

* -APPAY4(APGPNfl, APGPNT 2, APGPN (3»APGPnT4))*APFPAC1 

a +APHAY4(AkGPnT 1, APGPNT 2 , AKGPNT3,APGPnT4) 

T£MP7s(T£MP1WEMP14)*APFPAC2  + TEMP14 
TEmPos (l tMPl 3-1 EMP 12) aakFpAC2+TEMP  12 
TEMP^s( TEMPI  1- TEMPI 0)AAPFhAt2  + T£MP 10 
lEMp4=(TEMpP-TtMPbJAAP(-PAC2+TEMP8 
T E M P 3 s ( | EMP7-TtMPo J AAPFPAG3+  TEMP6 
TEMP2*(TtMPb-ltMP4)AAPFPAC3+IEMP4 
OtKi VE4s ( TEMP3- TEMP£) AAPF  PAC4+ IEMP2 
RETUPN 
END 


a 1 ME  FOLLOWING  ape  the  FUNCTION  TYPE  SUBPROGRAM  DEFINITIONS 
PEAL  FUNCTION  PDOTT(A,B,C) 

anote  Thai  a, a, and  c represent  p,w,and  k PtSPfc C f i vELY 
implicit  p£AL»4  (A-Z) 

COMMON  /APPAYl/UAl AwPDHlaOO) 

EUUIOALENCE  (UA TAwPOP(57) , IXMY) 

*, (OATAWPDP (60 ) , 1XY) 
a, (OATAWPDP (61 ) , I X/ ) 
a, ( DA  T AwPDP ( 62 ) , IYMZ) 
a, (DAT AwPDP (63) , 1YOY) 
a, (DAT AWPDP(65) , IZMX) 
a,  (DATAWPi)P(60)  , IZOZ) 

*, (UATAwPDP(68) , IF UNC) 
a,  (DAT  AWPDP  (81  ) , LBUD  Y ) 

A, (DA T AWPDP (8 7 ) ,M8UDY) 
a , (uA I AWPDP (86) / N600Y ) 

PDO 1 T = 1FUNC*( (LB00Y-IXYAAAC  + IXZAAA8+I YmZ*B*C  ) 

1 + IYOY*(MHUDY+iXY*6*C-IXZ*(A«r*2-C**2)  + IZMX*A*C) 

1 +IZOZa(NBOOY-1X  Y*  (B**2-A**2)-IXZ*B*C+IXMY»A*tt)  ) 
RETUPN 
END 


!. 


ri 

i; 
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REAL  EUNCTIUN  UOU  T f ( A , H , C » 0 ) 

•A,B,C*D  REPRESENT  P,U,H,POUT  RESPECTIVELY 
1MPLIC1I  REAL*4  (A-Z) 

COMMON  /ARRAY l/OATA*ROR(4O0) 

EQUIVALENCE  (OATAwROK(bO), IXY) 

«, (OA  I ArtROR (bl ) , I XZ  ) 

* # l U A I AwRl)H(b4)  , 1 Y Y) 

* » (0ATAWR0R(65) , 1ZMX) 

* » (DAI ArtRDP (87) ,M»OUY) 

QDU  | T*(l/I  YY)*(MB00Y  + IXY*(d*C  + U)-IXZ*(A**<?-C**2)*lZMX*A*C) 

RETURN 

ENU 

* 

REAL  FUNCTION  RO0TT(A,B,C,0) 

*A,B,C»0  REPRESENT  P,Q,R,POUT  RESPECTIVELY 
IMPLICIT  HE AL*4  (A-Z) 

COMMON  /AHWAY1/DATAwROk(400) 

EUU I VALENCE  (OATArtRDR(b7) , 1XMY) 

*, (OATANRDR(bO) , IXY) 

*, (0AIA*RDR(61), IXZ) 

*,  (OA  t AwROR (b7) , III ) 

*, (DATA*WDR(«8) ,NBODY) 

RDUTT=(1/IZZ)*(nb0DY-Ixy*(B**2-A**2)-IXZ*(6*C-D)+IXMY*a*H) 

return 

ENU 

* 

REAL  FUNCTION  THETUUT  T (A) 

* A REPRESENTS  phi 

IMPLICIT  WEAL*4  (A-Z) 

COMMON  /ARRAY1/UATAnHUR(400) 

EQUIVALENCE  (UATAWRU«(106) , UBODY) , ( DA T A wkDR ( 1 1 4 ) , RBOD Y ) 

THE  (UUTTsQBOOY*CUS(  A ) -RBUD Y * S I n ( A ) 

RETURN 

ENU 

* 

REAL  FUNCTION  PH  1 DO  I T ( A , a , C ) 

*A,B,C  REPRESENT  PH  I , T HE T A , PS  I DO T RESPECTIVELY 
IMPLICIT  Kfc  AL*  4 (A-Z) 

COMMON  /ARHAY1/DATANRDR(400) 

EQUIVALENCE  (OA  T AwHQK( 100) .PHUOY) 

Phi  Out  T*PBOOY*C*SlN(b) 

RE  I URN 
END 

* 

real  Function  psiootua,b) 

* A , B RE  PKE  SENT  phi, THETA  RESPECTIVELY 

IMPLICIT  RE Al * 4 (A-Z) 

COMMON  /ARRAY1/OATA«HDk(400) 

EQUIVALENCE  (0A1 A*RDW( 10b) , UBODY) , ( 0 A T AWRDR ( 1 1 4 ) , KBOD Y ) 
PSIOUT  Ts(HBODY*COS( A)>UBU0Y*SIN(A) )/COS(B) 

RETURN 

ENU 

* 

REAL  EUNCTIuN  SPEtO(A) 

* A REPRESENTS  AN  UNNAMED  F'JNCT  I t),x  / V A R I AbLR 
IMPLICIT  RE  AL*  4 (A-Z) 

IF(A.QE.l)  OU  TO  10 
STUREsl+A»*2/4+A**4/40*A*«b/lbOO 


135 


ooooooooo 


GO  TU  2 o 

10  STOWE3l.83P-,772/A**2t.lb4/A**4+.035/A**b 

20  SPEEDsSUWT (STOWE ) 

WE TOWN 

END 


POINT  FUNCTION 


SUkJWOuT  INt  POINT  F ( AWGV  ALOE  , AWGLNTH , A WGPN  T , AWFW  AC  , AWGL I ST  ) 

* THIS  IS  THt  FLOATING  POINT  VERSION  OF  POIM 

SOBKUUT1NE  PONT  COMPAHES  THE  FLOATING  POINT  VALUE  IN  AWGV  ALUE  I T H 
A list  OF  FLOATING  POINT  VALUES  startimg  in  location  AWGLIST 
AND  WE1UWNS  AWGPNT  AND  ARFWAC  TO  THE  M 1 1 N PWU7WAM 

rtHEwE  AWGLIST  IS  A LIST  OF  (AWGLNTH)  hEAL  NUMHEWS  IN  ASCENDING  UWDE 
ANO  whEFE  AWGPNT  IS  THE  POSIIIQN  OF  1 HE  LARGEST  VALUE  IN  AWGLIST  TH 
IS  LESS  THAN  AWGvALUE  AND  ARFRAC  IS  A SCALED  FRACTION  VALUE  REPWESE 
THE  DISTANCE  ttETwtEN  THE  TrtO  POINTS  AW7PNT  ANO  ARGPN1PLUS  ONE 

IN  I EGEW*2  AhGLNTh, AWGPNT 
WEAL  AWGVALUE, AWGLIST (AWGLNTH) 

2 IF  l AWGPNT ,LT . 1 ) GOTU  ISO 
4 IF  (AWGPNT, GT.AWGLNIh)  AWGPNTsAWGLN T H 
10  IF  (AHGVALUE-AWGLISTIAWGPNT))  70,130,20 
20  I TEHPsAWgPN I 1 1 

30  IF  (I  TEMP. G! .AWGLNTH)  GO  TO  130 
AO  IF ( AwgVALUE-AWgL 1ST ( l TtMp) ) 100,120,50 

SC  ARgPnTsITEMW 
bO  GO  TO  20 

70  IF  (AWGPnT-1)  ISO, 130,80 
80  AHGPN!=AWGPNT-1 
PO  GO  TO  10 

100  AWFWACSI AWGV ALOt-AWGL I ST (AWGPNT ))/ ( AWGLIST ( I T £MP ) - AWGL I S T ( A WGPN I ) ) 
110  GO  TO  140 

120  awgpnt*itenp 

130  AHFWACsO 

132  IF (AWGPNT ,L T .ARGLN fH)  GUTO  140 
135  AWGPNl sAKGLNTh-i 
137  AWFWAC*1.0 
140  RETURN 
ISO  ARGPNTsl 
160  GO  TO  10 
END 


# ################## 


NOTE:  Pages  136  to  148  inclusive  contain  all  the  aerodynamic  coefficients 
used  within  this  program. 
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-.006 

-.019 

-.021 

-.024 

.005 

-.007 

- , 020 

-.022 

CNH 

63 

-0.460 

-0.450 

-0.385 

-0.405 

-0.620 

-0.425 

-0.450 

-0.530 

-1.730 

-1.730 

-0.340 

-0.370 

-0 .645 

-0.520 

-0,520 

-0,290 

-0.985 

-0.815 

-0.580 

-0.580 

-0,620 

-0.655 

-0.490 

-0,400 

-0.240 

-0.220 

-0.250 

-0.250 

-0.230 

-0.240 

-0.220 

-0.250 

CYB 

63 

-.0182 

-.0167 

-.0171 

-.0168 

-.0078 

-.0195 

-.0180 

-.0165 

-.0070 

-.  0070 

-.0182 

-.0150 

-.0095 

-.0105 

-.0090 

-.0155 

-.0095 

-.0100 

-.0100 

-.0075 

-.0085 

-.0080 

-.0080 

-.0090 

-.  0080 

-.00  70 

• . 0068 

-.0068 

-.0100 

-.0060 

-.0070 

-.0068 

.00064 

-.00053 

-.00042 

.0002  1 

.00046 

-.0004 2 

-.00034 

. 00027 

.00040 

-.00035 

-.00030 

.00023 

.00037 

-.00033 

-, 00026 

. 000  1 9 

.00036 

-.00030 

-.00024 

. 000  15 

.00000 

-.00036 

-.00030 

.00024 

.00000 

. 00000 

-.00108 

.00102 

.00084 

-.00072 

-.00060 

.00092 

.00070 

-.00064 

-.00052 

.00040 

,00049 

-.00040 

-.00031 

. 00007 

.00042 

-.  00034 

-.00025 

, 000  1 b 

.00039 

-.00031 

-.00023 

. 000  15 

.00031 

-.00025 

-.00017 

. 00008 

.00035 

-.  00029 

-.00022 

. 000  1 4 

.00000 

-.00040 

-.00039 

- 

.00038 

.000  18 

-.00004 

-.00038 

. 00036 

.00022 

-.00008 

.00000 

. 00037 

.00018 

-.00014 

-.00004 

. 00000 

.0001  7 

-.000  13 

-.00009 

,00001 

.00020 

-.0001  3 

-.00010 

. 00006 

.00024 

-.  000  1 6 

-.0000/ 

.00004 

.00029 

-.  00023 

-.00015 

,0000b 

. 000U36 

.000035 

.000034 

,000009 

.000024 

.000028 

.000027 

.000026 

.000012 

.0000  1 4 

,0000  16 

. 0000  1 6 

.000005 

.000007 

.000010 

,00 

.000000 

.000002 

.000004 

, 000006 

,000000 

,000000 

,000000 

, 0 0 o 0 0 0 

.000000 

,000000 

,000000 

, 0O0000 

.000000 

, 000000 

.000000 

.053 

.026 

.020 

-.005 

.062 

.033 

.020 

,006 

.054 

.032 

.010 

.001 

.049 

.034 

.015 

-.005 

.051 

,0  36 

.021 

,004 

.024 

.034 

. 020 

, 006 

.025 

-.026 

.030 

.017 

.024 

-.024 

-.028 

0,620 

•0,620 

•0,620 

-0,620 

0.700 

-1.010 

-1 . 390 

- 1 .665 

0.420 

-0,660 

-0.910 

-0.900 

0.310 

-0.380 

-0.550 

-0.770 

0.220 

-0 .240 

-0.300 

- 0 . 4 4 0 

0.400 

-0.335 

-0.240 

-0.230 

0.250 

-0.250 

-0.335 

-0,240 

0.250 

-0.250 

-0.250 

.0  135 

-.0104 

-.0088 

-.  0078 

.0145 

-.0113 

-.0100 

-.  0084 

.0135 

-.0125 

-.0115 

-.0105 

.0140 

-.0130 

-.0116 

-.0  105 

.0140 

-.0128 

-.0120 

-.0100 

. 00  75 

-.0113 

-.0110 

-.0100 

.0068 

-.0068 

-.0113 

-.0110 

.0068 

-.0068 

-.0  068 
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CYOA 

-.00 JUQ 

378 

-.00326 

-.00318 

-.00284 

-.00250 

-.00025 

-.00025 

-.00025 

00025 

-.00222 

-.00171 

-.00121 

-.00086 

-.00051 

-.00025 

.00000 

.00000 

.00000 

-.00031 

-.00023 

-.00015 

-.00007 

.00000 

.00006 

.00016 

.00027 

.00036 

-.00028 

-.00021 

-.00014 

-.00010 

-.00006 

-.00003 

,00000 

.00023 

.00023 

-.00010 

-.00008 

-.00006 

- , 0 o 0 05 

-.00003 

-.00002 

.00000 

,00006 

,00016 

-.0001  1 

-.00007 

-.00004 

-.00004 

-.00004 

-. O00O2 

.00000 

.00007 

.00014 

-.000  15 

-.00010 

-.00006 

-.00004 

-.00002 

,00000 

.00002 

.00006 

.00014 

-.00270 

-.00222 

-.00174 

-.00118 

-.00062 

-.00012 

-.00012 

-.00012 

-.00012 

-.00081 

-.00056 

-.00031 

-.00016 

-.00002 

.00008 

.00016 

.00016 

,00016 

-.00006 

-.00008 

-.00007 

-.00003 

.00001 

.00004 

.00006 

.00017 

.00017 

-.00007 

-.00005 

-.00002 

.00000 

.00003 

,00003 

,00003 

.00006 

,00006 

-.00008 

-.00004 

.00000 

.00002 

.00004 

.00006 

.00006 

.00014 

.00020 

-.00006 

-.00005 

.00000 

.00000 

,00001 

.00004 

.00007 

.00012 

.00018 

-.000  1 1 

-.00006 

-.00001 

.00000 

.00001 

.00004 

.00007 

.00013 

.00016 

-.00180 

-.00116 

-.00058 

-.00040 

-.00022 

.00012 

.00012 

.00012 

.00012 

-.00042 

-.00024 

-.00006 

.00006 

.00018 

.00020 

.00022 

.00022 

.00022 

.00000 

.00000 

.00000 

.00002 

,00005 

.00007 

.000  1 0 

.00021 

.00021 

-.00004 

-.00001 

.00002 

.00005 

.00006 

. 000  1 1 

,00013 

.00018 

.00018 

-.00005 

-.00001 

.00003 

,00006 

.00006 

.0001 1 

.00013 

,00017 

. 00022 

-.00007 

-.00003 

,00002 

. 00O04 

. uOOOb 

.00008 

.00012 

.00017 

.00022 

-.00009 

-.00005 

.00000 

.00002 

,00005 

.00008 

,00011 

.00016 

,00022 

-.00054 

-.  00036 

-.00016 

-.  000  1 0 

-.00002 

.00078 

.00078 

.00076 

.00078 

.00000 

,00030 

.00030 

.00036 

.00046 

.00064 

.00060 

.00080 

.00080 

.00015 

,0002b 

.00038 

. 00046 

,00060 

.00072 

.00085 

.00110 

.00110 

-.00001 

.00006 

.00020 

.00025 

.00030 

.00036 

.00044 

.00060 

.00060 

-.00008 

.000O1 

.00010 

.00017 

.00025 

.00030 

.00035 

.00042 

.00050 

-.00010 

-.00001 

.00008 

.00015 

. 00022 

.00027 

,00032 

.00042 

.00052 

-.00012 

-.00004 

,00004 

.000  1 1 

.00016 

.00025 

.00032 

.00042 

.00052 

.00000 

.00000 

.00000 

.00000 

.00000 

-.00001 

-.00001 

-.00001 

-.00001 

.00001 

.00001 

.00001 

,00001 

.00002 

.00002 

.00003 

.00003 

.00003 

.0000  1 

,00001 

0.00001 

.00000 

.00000 

.00005 

.00011 

.00066 

.00066 

-.00012 

-.00003 

.00006 

.00011 

.00016 

.00021 

. 00026 

.00040 

,00040 

-.00005 

,00006 

.00017 

.00024 

.00031 

.00036 

.00042 

,00050 

.00058 

-.00004 

. 00005 

.00016 

.00025 

.00035 

.00036 

.00041 

.00050 

.00060 

-.00006 

.00001 

.00008 

.00015 

.00022 

.00031 

.00040 

.00050 

.00060 

.00034 

.00023 

.00012 

. 000  12 

,00012 

-.00168 

-.00168 

-.00168 

-.00168 

-.00020 

-.00045 

-.00071 

-.00066 

-.00100 

-.00125 

-.00151 

-.00151 

-.00151 

-.00038 

-.00052 

-,0006b 

-.00107 

-.00148 

-.00124 

-.00100 

-.00005 

-.00005 

-.00060 

-.00048 

-.00036 

-.00025 

-.00014 

-.00012 

-.00010 

-.00001 

.00000 

,00000 

.00010 

. 00020 

. 00030 

.00040 

,00045 

.00050 

.00055 

. 000b 1 

-.00001 

.00006 

. 00016 

.00026 

.00036 

.00045 

.00052 

.00061 

.00070 

.00000 

.00005 

.00011 

.00014 

.00028 

.00038 

.00048 

.00058 

C YDR 

.00205 

.00066 

63 

,00160 

.00175 

.00175 

.00174 

.00170 

.00165 

.00155 

.00145 

.00160 

.00150 

.00136 

.00134 

.00126 

.00122 

.00  1 14, 

. 00 094 

.00084 

. 001  12 

.00106 

.00066 

.00060 

.00081 

, 0 0 0 6 4 

,0005b 

.00017 

.00000 

.00083 

.00077 

.00070 

.00058 

• ,0004>» 

.00034 

.00023 

.00006 

.00000 

.00068 

.00062 

.00055 

.00041 

.00026 

.00016 

.00010 

.00002 

.00000 

.00060 

• 0 0 u5 1 

.00043 

.00026 

.00015 

. 00007 

.00000 

.00000 

.00000 

. OOOoU 

.00051 

.00043 

.00026 

.00015 

.00007 

. noooo 

.00000 

.00000 

CYDROSttl  7 
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APPENDIX  D - FLIGHT  CONTROL  SYSTEM 
PROGRAM  LISTING 


SUHPUUTINt  SHTLFCS 
IMPLICIT  PEal  *«(A-Z) 

I N r t Gfc K *«>  1 A1  , IPI  , lot  » 1GB 

L0GICAL*1  oatahyt 

COMMON  / APPA Y 1 /OA T ArtPOP ( «UU ) 
COMMON  /APPA Y3/OATABYT ( 17) 
COMMON  /APpAY«/AnGLE(3) 
LOGICAL* 1 AUTOttF 
L OG I C AL  * 1 IL  » HOLU 
LOGICAL*!  TUPNCOPD 
tool  VALENCE 

* (DAUwPDK(3),AY) 

* , (OATawpOP ( a ) , AZ  ) 

*, (0ATAwP0B(9) ,Bt T A) 

*, COAT  A wPOP ( 1 0 J r OSH  J 
*, (OATA»POP(3«),DA) 

* r (OATAwPOP ( 35) « GEL ) 

*, (DATAWPOP(37),PRHCSOP) 

*, (OATAwPDP  (3t>) ,0EK) 

*, (OA) ArtHOH(aO),OH) 

*, (OAT  AwKOP (43) #0ACM) 

*, COAT  ArtPOP (49) # Ot  ) 

*,  (OAT ArtP0«(5!  ) ,OBF  ) 

*, (uATAwB0P(54),DPJET) 

*, (OAT ArtPOP (55) > OP  JE  f ) 

* * (OAT  AwkOP ( 72 ) » PphCSOP ) 

*, (OAT ArtPOP (74)  , PPTASUP) 

*, (0ATA*P0P(79),UZCMD) 

*# (OAf ArtPOP (HO)# UXCMU) 

*, (OAT awP0P(63) »UYCMO) 

*,(0ATArtB0P(H4),MACH) 


( DA  I AwkDk ( 1 00  ) 
(DAT  AwkDk (102) 
(DAT AwkDk( 104) 
(DA (AwkDk( 10b) 
(DATAwkDK(106) 
(DATAwkDk( 109) 
(DAI  A wRDk  (112) 
(DATAwRDK (119) 
(DAlAwRDH(124) 
(DAT  AwkDk ( 131) 
(DAT  AwkDk (142) 
(DA  I AwkDk (148) 
(DAI ArtKDK(laO) 
(DAT AwkDk ( 1 44) 
( D A T A wkUk ( 1 46 ) 
(DATAwkDk(151) 
(DATAwkDk( 154) 
(DAT  AwkDk (200) 
(DAT AwHDH(201  ) 
(DATAwRDk(202) 
(DAI AwkDk(203) 
(DAT AwRDk(204) 
(DAT  A wkDR ( 205 ) 
(DAI AWRDk (207) 
(DAT  AwkDk (209  ) 
(DAI  A wkDk ( 2 1 0 ) 
(DAIAwHDR(21 1 ) 
(OATAwkDk (212) 
( DA  I AwkDk ( 2 1 3 ) 

( DA  I AWRDk  (214) 
(DAI  A WRDk (215) 
( DA  T A wkDR ( 2 1 6 ) 
( DA  T A wkuk ( 2 1 7 ) 
(DAT AWkDk (218) 
(DATAwkDk (219) 
(DAT  AwkDk ( 220  ) 
(DAT AwRDk(22i ) 
( dA I AwkDk ( 222 ) 
(DAT  A WRDk (223) 
(DAT  AwkDk ( 224 ) 
( DA  1 AwkDk ( 225 ) 
(DATAwRDk (22b) 
(DAI AwkDk (227) 
( OA  T A wkDR (228 ) 
(DAT  AwkDk ( 229 ) 
(DAT  AwkDk (230 ) 
( DA  T A wkDR ( 23 1 ) 
(OA I AwkDk (232  ) 
(DATA WRDR (234) 

(DATa*RDR(234) 
(DA  I AwkOk (245) 
( DA  T A wRDk ( 236 ) 
(DAT  A wkDR (237) 
(DAI AwHDK(238) 
(DAT AwRDR(239) 
( DA  I AwkDk (240) 
(DAI AwkDk (241 ) 
(DAT  AwkDk ( 242 ) 
(DAI AwkDk (243) 


> PBUD  Y ) 
,PE) 

, UUAk) 

, DDO  I ) 

, UbUDY) 
,KDOT  ) 
t SBHP ) 
,RdODY) 

» UYliJfc.  I ) 

# IAS) 

,PI TCHkAN) 
, xbAk ) 

,PI ICHSIK) 
/ kOCLPAN) 
.KULLSTtV) 

, YAwIklM) 
,PLAPCMD) 

, DEMAN) 

» e Shape ) 

f GPS) 
,DEMS) 
t u) 

i PM0D26 ) 

,GJtn 

,E  rkIMlN) 
rETkIM) 

, ELFBk IN) 

» ELPriK ) 

,DE  Tmpan) 
,dS8PC) 

» DSb  X T k ) 

, DSbXTkS) 
t DE I DSb ) 
,0t ikiM) 

, GDSd ) 

» GTRE ) 

, DE IMkHC) 
,DETR) 
t DE  MSP ) 

» R TPH  I ) 

, k 1 ANPHI ) 

,DJkIP) 

,tttSL) 

» DC  SL  ) 

# DQC  I ) 

t DEC"D ) 
.uM-  ID 
,OCSU) 

» DCSLLU ) 
,DCSLH1 ) 

# Ml  DPI Ck  ) 

# KP  I T ) 

, GDQ) 

t DULO  ) 

» k ) 

» DQh  1 ) 
t DUM I A ) 
.OUrtlN) 
,N2A) 
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* , IDA  I AaROH  (24  <4  ) , N2P) 

* , IDA I AnWOK (tub) i OULUN) 

*,  (OAT  A AH  OH  (24b)  , D (JLO  A ) 

* , (DA  I AwNDk (247) , ALPHMIN) 
*,  (OA I AWHDK (248) , ALPHMAX) 
*,  (OATAwHok (249) , GTKfcOOCT ) 

* , (OAT AwH0H(2bQ),P) 

* , (DATAbHUH(2b2) ,UA  MWHC) 
*, (OATAaHoK(2S3),PKQ024) 

*,  (DAIAAkOK(2b4),PHU025) 

*, (0ATAWMUH(2S7) ,PH0021 ) 

*, (OAI A*RDH(2b8),PH0D22) 

*, (OAI a ahoK ( 2b9 ) »PR0023) 

*, (OAT  A aRDH ( 2b0 ) , PHUD27 ) 

*,  (OAT  A*KOR(2bl  ) ,SUM21  ) 

*,  (0ATAinP0N(263)»SUm23) 

*, (OAI  ArtKOW ( 2b4 ) ,DCSP) 

*, (OATAwRok (265) ,OAC*D) 
*,(DATA*ROk(2bb),OATSUMl) 

*,  (DAT  AwWL>P(267),KDAMLIN) 

* , (DA  T A aHDR  (268)  »*DAMPAR) 
*,  (DA  F AaHDH  ( 2b9 ) » GRS) 

*, (OAT AwRDk (2/0) ,GDAC) 

*, (DA) A aROH (2/1 ) , GO  A ) 

*, (DAIAaHDH(272),DAT5UM) 

*,  (0ATA»vk0H(2/i),PK0042) 

*, (DAI AaKDH(274),DATR) 

*, (OA  T Artkuk (275) , OkPRM ) 

*, (DATAwHDH(276),PSTA8) 

*, (OAI A*ROR(27  7 ) ,SUM24) 

*, (DATA WHDH(27»),DH PHI) 

*, (OA I A a Hurt (279) ,PRG029) 

* , (DAlAwRON(2ttO),  SU*^) 

*,  (DAT  AaHI)h(281  ) » G T X ) 

*, (DATAwKDK (282) , YAwXFfcEO) 
*, (DATA *RDR( 283) ,PK0041 ) 

*, (DATA  aHOH ( 284 ) ,GP) 

* , ( OA  T A aHDR ( 28b  ) ,GRH) 

* , (OAT  A WHOM (28b) , SUM25) 

*, (0ATArtHDk(289),SUM2b5) 

*, (0 A T Aahdh (291 ) , SUM26) 

*, (OATAaHOH (292) ,PH0043) 

*, (DA  I AaHOH (293) , HST AB) 
*,(0ATA*WuH(294),t)tTAo) 

*, (OAT  AaHOH (29b) ,ttt  TAFILT  ) 

* , (OATAaHDH (296) ,DATMPAN) 

* , (OAT  ArtHOH ( 29  7 ) ,OATP) 

*, (DA  T A a R D N (299) , P R 0 0 4 4 ) 

*, (OA I AWHDR ( 300) ,DK f ) 

*, (OAT  A aPDH (301 ) , OR  T MSF ) 
*,(OATAaHoW(302),0HmS) 

•, (OATAaHOH(303),ORTM3) 

*, (OAI AWKDR(  304) , PROD  1 ) 

* , (OAT  AaHDH (305) , SUM1 ) 
*,(OATAAHOH(30b),NYP) 

* , (0ATAaHoh(307),WYA) 
*,(OATAaHuH(308),GHAY) 

*, (DAT  aaHDH ( 309) ,PR0D2) 
*,(DATAahDHv310),PHO03) 


w 
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• , ( OA  f AivRQR  (ill  ),SUM2) 

*,  (OAT  awRDR(312)  , PRODS ) 

*, (OAT  AWHDR ( 313) fOHCPF  J 
*, (OAT AWHOR(314) ,*ORC) 

• , (DATAnRDk( 315),G0RC) 

*,  (0ATA*RDR(316)  , GNYDRm) 

*, (UA  r AwRUK (31 n , OK  TKK) 

*, (DAT AWROR (3l8) , PROD  7) 

*, (DATArtROK(319),GTRR) 

*,  (DAU*RDR(320),DR1R1M) 

*,(OA(A«*ROR(321),SUM3) 

*, (OAT AWROK (322) , ORCMO) 

*, (DATAwRDR(323) ,DACMC) 

*,  (DATA*kUR(324)  , GR  XF  D ) 

*, (DATAwROK(325) ,GDR£) 

*, (DATAwRDR(327),nR) 

*, (OATArtY T ( 10) , hOLO ) 

*»  (OAT  AtlYT  ( 11  ) , IC) 

• , (DATAbY ( ( 13) » 1URNCOHU) 

• , (oat  Ayr t ( 14) , autubf) 

• » ( ANGLfc ( 1 ) , ALFA) 

*,  ( ANGLE  ( 2 ) , PH  I ) 

• » (ANGLE (3)  , T Ht  TA) 

DIMENSION 

•GP3T (4) , 

*ALPHMINT(4), 

•TALPHMAX (4,2)  , 

• TKPI T (4,2)  , 

*G  T R£  T (4), 

• GDSfclT (4)  , 

• KDAMLINT (4)  , 

•KDAMPAKT (4)  , 

•GR3T (4) , 

• GOAC  T (4)  , 

*GPT (4) , 

•GNYORMT (4) , 

• GRAY  T (4)  , 

•GDHC) (4), 

•GTRR1 (4) 

DATA  DEGRAD/57.29577951/ 

0 A I A GPSf/0. 4,0. 123,2.0,2.5/ 

DATA  ALPHMIM/-4. 0,0.0,  .6,  ,8/ 

DATA  T ALPHMAx/20.0,  15.0,  15.0,30.0, .6, .8, 1.5, 3.0/ 
DATA  T*P1T/15.Q,81.5,81.5,60.0,1.2,3.0,3.5,5.0/ 
DA i A GTRfcT/1 .0,0.3, 10.0,  12.0/ 

DATA  GOSa?/0.1,0.25, 1.5, 3.0/ 

DATA  KDAMLIN  T/0.08, 0,9324, 0.6,  1 .5/ 

DATA  KUAMPAHl/0.0636,0.008,0.6,1.5/ 

DATA  GRST/1.0,0.25,0.6, 1.5/ 

OATA  GOALT /0.45,2.5, 1 . 15, 3.0/ 

DATA  GPT/10.0, 1 .5, 10.0,45.0/ 

DATA  GNYORMT/O. 05,0.015,250.0,900.0/ 

DATA  gkaYI/10.0,3.2,2.0,3.0/ 

DA (A  GDRC  T / 1 200.0,3800.0,1.15,5.0/ 

DA  I A gJRRI/0,2, 1.0, 2. 0,5.0/ 

DATA  UAMAx/20.0/ 

DATA  PRHCOH/1.15/ 

DA  1 A GPRrtCT/0.3/ 

OATA  HP 1ADU/1.125/ 
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0 A I A NUOO/2.0/ 

DATA  NQHl/4,0/ 

DATA  UPCUIOFF/1 .0/ 

DATA  DP Jt  T 1/0.0/ 

DATA  ULDPt/0,0/ 

DATA  DFIJET1/0.0/ 

DATA  UdDB/0.25/ 

DATA  UAUH/0.2/ 

DATA  I G J E f 1 2/ 0 , 5/ 

DATA  TGJET21/5.0/ 

DATA  GHPAN IL/0,5/ 

DATA  GHPANTE/-0.5/ 

DATA  GPHrlC  TE/0.6/ 

DAIA  GPKHClL/0.2/ 

DATA  GUMhGP / 1.5/ 

DATA  GUA/1.0/ 

DATA  NZMAX/3.75/ 

DATA  GUNL/4,0/ 

DATA  NZM1N/-1 .0/ 

DATA  bSHHl/25.0/ 

DATA  GSdLU/10.0/ 

DATA  GPPAnT/2.0/ 

DATA  GPP/1.0/ 

DATA  GlPA/u.2/ 

DATA  G/32.174/ 

DATA  PPT A0b/1 . 125/ 

DATA  bDK 11/1.0/ 

DATA  YBDFlLT/l ,0/ 

DATA  DDF  AUd/0,5/ 

DATA  DbFdDd/l .0/ 

DATA  kDK^L IN/0,0131/ 

DATA  KOPMPAK  / 0 . 042/ 

DATA  DDF  HA/1. U/ 

DATA  DBFHB/2.0/ 

DATA  DFkADH/0.9/ 

DATA  OFftBOd/ 1,0/ 

DATA  PtdOd/O.2/ 

DATA  PE  ADB/ 0.1/ 

DATA  NB 1/1,0/ 

DATA  MW2/2.0/ 

DATA  NH 3/3.0/ 

DATA  hiR4/4.0/ 

DATA  GWCSA/-2.0/ 

DATA  DBFDCS1  /U.O/ 

DATA  uBFDC  /O.O/ 

DATA  01,025/0.0/ 

DATA  DSbtl/0,0/ 

DATA  HI  A>\(PhiN/0.0/ 

DATA  BCSLN/O.O/ 

DATA  UFFILh/O.O/ 

DATA  N/PNl .NZPN2/0, 0,0.0/ 
DATA  PBOU22N/0 . 0/ 

DATA  PK(J024i*/0 ,0/ 

DATA  t)fc  T AN/  0 , 0/ 

DATA  nYPnI  ,NYPrJ2/0. 0,0.0/ 
DATA  PWOD2N/0.0/ 

IF (ICO  GO  ID  1000 
3 CONTINUE 
NYsAY/32, 174 
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j NZsAZ/3*!.  1 7 a 

1 ALPHArALFA*OtGRAD 

CGPUSsl£9i.i6*XbAR+£46.0 

PbPB00Y*DEGRa[) 

UeUdOOY*OtGKAD 
k=kbUUY*UEGkAO 
SINALP*SlN( ALF A) 

S1NALF IbSINALF 

1 IF (SINALF .Ll .0.087155)  S I N ALF 1 aO . 08 7 1 55 

COSALFsCOS(ALFA) 

| S1NPHI«31N(PHI) 

||  * C0SPHI*C0SIHH1) 

‘ I COST  HE  T A»COS ( THE I AJ 

DETMPANsQ.O 

IF  (PiTCHPAN.Gt.O.SJOt IMP ANsl .0 
IF (PI TCHPAN.LE.-0.S)DETMPAN»-1 .0 
OtTMKMCsO.O 

IF (P1TCHSTK.GE.O.S)OETMRhC=1  .0 
IF  (PUCHSTK.Lfc  .-O.S)DETmkhC3-1 ,0 
DATMPAN*0.0 

IF  (Ri>LLPAN.GE.0.5)OATMPAN=1 .0 
IF (RULLPAN.LE .-0.5)OATMRANs-l . 0 
DATMkHC=0.0 

i IF(KOLLSlK.Gt.0.5)DAIMKHC=1.0 

IF  ( ROLLS  T Lt  . -0 . 5 ) 0 A T MRHCs- 1 .0 
okrso.o 

IF ( Y A* TRIM. Gt. 0.5) OR T*1 ,0 
IF (YArtTklM.LE.-O.S)OHTs-i.o 
UBFMANsO.O 

IF (FLAPCMO.GE.O.SJOdFMANsl  .0 
IF (FLAPCMO.LE.-O.S)DBFMANb-I .o 

* 

* 1 ME  FOLLUWlNG  SECTION  CALCULATES  ALL  TABLE  LOOKUP  FUNCTION  VALUES 

* 

UPSsFUNCT I UN (MACH, GPS T J 
GJt  f*0.0 

IF (UbAW.&t. TGJETldJGJETsTGJtT^l 
CALL  POINIF  (MACH, 4,  I A T , F K AC , T A L PHM A X ( 1 , 2 ) ) 

ALPHMAXST ALPHMAX l I A T , 1 ) ♦( T A LPHM A X ( I A T ♦ 1 , 1 ) - 1 ALPMMA X ( 1 AT , J ) ) *FR AC 
alPhm]N=FUNCT 10n(maCH, ALPHMINT ) 
call  PUINTF (MACH,4, IPT,FkAC, TKPIT ( 1,2)) 

KPJTbTKPIT ( I P T , 1)*(TKPIT ( IPT ♦ 1 , 1)»TKPIT (IPT,1))*FRAC 
G T RE  =F  UNC I I ON ( MACH , G t kt  T ) 

GOSBsF  UNC  T ION(MACH,GOSBT ) 

K 0 AMP A R3F UNCTION (MACH, NO AMP ART ) 

KOAML INbFUNCTIuN(MACH,KOAMLINT ) 

GRSsFUNC T ION(MACH,GRST  ) I 

GOACsFUNC T 1UN(MACH,GDAC1  ) ! 

GPsFUNCT  IUN((JBAR,GPT) 

GNYUkmbFUNCTION(TAS,GNYURmT)  I 

GRAYSFUNCT IUN(MACM,GRAYT) 

KDRCsFUNCT IUN(MACH,G0RC r ) 

GfkksFUNL  T IOn(MACm,GTHHT ) 

* 

* PITCH  FCS  PRUCESSING 

* 

* OEMSP 

OEMAnbO 

IF  (PRHCSOP.LT ,-PHHC OB) OEMANsPRHC SOP ♦PRHCUB 

155 


■■■i  m 


IF  (PKtlCSOP  ,GT  .PHHCOb)OfcMAN*PWHCSOP-PPHCDb 

ESHAPEs(0.36+0.Q484*AbS(D£MAN))*OEMAN 

IF  IE8HAPE.GT  ,23.Q)EShAPE=23.0 

IF  (ESHAPE.LT  ,-23. 0)E8HAP£s-23.0 

0£M3sESHaPE*GPS 

DEMSPaUtMS 

IF((MACH.GT.1.5) .AND. (PHHCSUP.GT . 19.5) ) OEMSP*0£MS*2 . 0 

*DE  1 P 

20  Dfc  TMGPPsOE  TMPHC*GPPMC  I 
IF(IC.OR.MUID)  GO  TO  29 

UE  T P*F 1 L T 1 (DETMGPP, .09,0. 0,-1 .0,DETHN) 

IF (DETP.GT . 1 .5)  D£  TH*1  ,5 
IF<UETP,LT.-1,5)  OETPs-1.5 
29  CONTINUE 
•DPJET 

OPJE  T*t)EMSP-U 

*UYCMD 

AbSDPJET  = AbS(OPJ£T ) 

NQsNQH  J 

IF ( AbSOPJET .LE.OPCOTOFF)  NU*NGLO 

CALL  HYSIEK(UAD0,ObDH,DPJETl , DP JE T , OUTDPJE T ) 

OPJET laDPJET 
UYCMD=NU*OuTOPJET 

IF ( (M6CH.LE. 1 .5) .OP, (OBAM.GT .20,01 JUYCMDsO.O 
* U Y G J E T 

50  UYGJETsGJET*UYCMD 
*PT ANPHI 

TANPMlsSlNPHI/COSPHl 
IF ( TANPrtI ,GI . I .0) T ANPHIsl  ,0 
IF (TANPHI .LT.-l .0) TANPHl*-i ,0 
60  PTPHI»R*T  A NPH I 

IF(hOLO)  GO  TO  62 

PT ANPHIsF ILll(H(PHl,.0196l, , 0 1 96 1 , - , 9608 , P T ANPH1 N) 

62  CONTINUE 
*BCSL 

IFIMACH.GT, 1 ,5J  HTAnPHIsO.O 
70  DJPTP=-OPJE1-PTANPHI-OETP 
I F ( HOLD ) GO  TO  72 

BCSLaFILT HO JR  TP,  1 .9615,-1 .8846, -.9231 ,bCSLN) 

72  -CONTINUE 
•DCSLLU 

IF(MOLD)  GO  TO  74 

OFF  I L *F 1 L T 1 (U, . 02988 , . 02988 ,-. 992 , UF F I L N ) 

74  CONTINUE 

OCSUSUF  F IL*GUMwGP 
DULOAs( ALPHA-ALPHMAX)*GQA 
NZ  AsNZ ♦ 1 ,8b49*U00T 
IF (HOLD)  GO  TO  76 

NZP»F ILT2(N/A, .06475, . 1 295 , . 064 75 , - 1 . 30 94 , . 5683 , NZPN 1 , NZPN2 ) 
76  CONTINUE 

OOLUNs (NZP-NZMAX ) *GONL 
OOLU=OuLOA 

IF (UULUN.GJ .OOLOA)  OOLUxOULON 
80  OCSLLOaOQLO+OCSO 
•OCSLHl 

DOh]As(ALPHA-ALPHMlN)*GQA 
OOhlNS (NZP-NZMINJ  *GONL 
OUHI =OOHl N 

IF (UUHI A.LE.OUHlNlOQHIsOOHl  A 


90  DCSLHl=OUHl*DCSu 
•MIDPICK 

IF  (DCSLLO. Gt. DCSLHI ) GO  TO  98 
IMDCSLHI.LT. BCSt)  GO  To  91 
GO  10  9 i 

91  IMOCSLMl.Gf  .DCSLLO)  GO  TO  97 
GO  TO  98 

93  IMBCSL.Gt. DCSLLO)  GO  TO  98. 

GO  fU  99 

97  MIDPICKsDCSLHI 
GO  TO  100 

98  MIDPICKsBCSL 
GO  TO  100 

99  MIOPICKxDCSLLO 

• DC  SL 

100  IF (MACH.GT . 1 , 5 ) GO  TO  105 
DCSLxMlDPICK 
GO  TO  110 

105  IF  (UBAK.GT .2.0)  GO  TO  108 
OCSLsO.0 

GO  TO  110 

106  OCSL=bCSL 

• DUCT 

110  CONTINUE 

GDUSKP1 T/SUHT (QBAH+4.0 ) 

IF (GDU.G£.6.0)GOUS6.0 
IF  (GUO. LE. 0.2)UDQ*0.2 
DUCT*DCSL*GDU 
•GTKEDOCT 

GTHEUUCT*GIWt*OUCT 

• DSB  X T KS 

DSBXThfsOSBPC*GDSd 
IF (DSB. Gt. 30.0)  GO  TO  115 
DSbXTWSsOSBXTk*GSBLO 
GO  T U 120 

115  OSdXTKSaOSBXTH*GSBHl 
•Dfc  TP 

120  OE  TP*0E  TMpAN*GPPANT 
DETDSBaDSBxrHS-OETP 

•ETPIM 

ETHlMlNsGIHfcOQCT-DETP-OYGJET 

if(ic.oh.hold)  go  to  i a i 

ETWIMeF IL1 1(ETHIMIn,.09,0.0,-1.0,ETH1MN) 

121  CONTINUE 

IF(ETPIM.&t.20.0)tTKIMsao.O 
IF  (fc TPIM.LE.-35.0)tT«IMa-35.0 

•ELFBK 

ELFBKINs(OfcTOSB/GT«E)*DE 
IF (1C .OK.HULD)  GO  TO  122 
TEMPsbO , 0/GT  HE 

ELFBKsF ILT 1 (ELFBK IN, 1 . 0/ ( TtMPt 1 . 0) , 1 . 0/ ( TfcMPf 1 . 0 ) , 

* ( ” T tMP* 1 ,0) /( T£mP*1 ,o) ,ElFBKN) 

1 22  CONTINUE 

•detpim 

IF  (MACH, LE.  12,0)  GO  TO  125 
DETKIMxt 1RIM 
GO  TO  130 
125  OETKIMstLFBK 
130  CONTINUE 
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•DECMO 

140  DECMU=DETklMtOUCT 

IFCDECmU.GT.20.0)  UECMD=20 . 0 
IF COELMO.L T .-35,0)  OECMO*-3S.O 

* 

* END  OF  HITCH  PROCESSING 

* 

* ROLL  FCS  PK0CESSIN6/PART  1 

. - a 

* DAMAN 

DAMANsO , 0 

IFCRKHCSOP.LT. -1,1b)  UAMANsKKHCSQPf 1 , IS  . 1 

IF CKKhCSOP.GT  . 1,15)  D AMANaRRHC SOP-1 . IS 

* PR0D21 

PRU02 l=DAMAN* (KDAMLlN+NOAMPAKoAbSCDAMAN) ) 

IF (PHD 021 ,GT .OAMAX )PK0U2ls0AMAX 
IF CPK0U21 .L  f ,-0AMAX)PK0D2ls-DAMAX 
1 F ( HOLD ) 00  TO  155 

PK0022  = F IL  T 1 (PK0021 , . 0909 , . 0 909 , - . 6 1 82 , PR0022N ) 

155  CONTINOE 
*GKS 

PR0D23sGKS*PKUD22 

IF (MACH.G  T . 1 .5)  PRUU24sOA?MKHC*.6 
IF( MACH. LE. 1.5)  PW0024sOATMKHC*.2 
IF C I C . OK. HOLD)  GO  TO  1 7 1 

PK0025=F 1LT 1 CPKUD24, . 0 4 , 0 , - 1 . 0 , PK0U25N ) 

171  CONTINUE 

IF (Pk0D25,GT.2.5)PK0D25s2,5 
IFCPK0D25.LT.-2.5)  PR0D25*-2.5 
SUM2l=PKOD25+PROD23 

PKOD2('3SUm^1*2,0 

IF(PKGD27.GT.10.0)PRQD27sl0.0 

IF(PK0027.LT,-10.0)PK0027a-10.0 

PR002b=PK0D27*SlNALFl 

* 1)  R P H I 

DKPHI=57,3*G*SINPHI*C0STHE1 A/TAS 
IF  (TUKNCOKO)GUTO  17*4 
ORPH 1=0 

* DKPkM 

174  DKPhm=R-ORPHI 

I F C HOLD  ) GO  TO  17  3 

PK0D29=F IL T 1 CDRPRM, , 9804 .9804 , 9608 , PW0D29N ) 

178  CONTINUE 

* GKh 

GRh=28,8/ ( OB AH ♦ 10.0) 

IFCGKH.LT. .05)  GRH=.05 
IFCHOLD)  GO  TO  185 

HETAF ILIsFILT 1 ( BE T A , . 999 , - . 999 , - . 998 , BE T AN ) 

185  CONTINUE 

BE  T Au=BE  T AF  1LT 

IF  C08AR.LE .2.0) BET AU*BET A 

IF  COBAR. GT .20.0) BE TAOxO.O 

SUM22  = C0SALF  *BETAU-PR0D29*GRH-2*DRPRM 

SUm23=PK0D26+SUM22 

* RSTAB 

RST ah=okpkm*cosalf-p*sinalf 

* PSTAtJ 

PST AB=0KPKM*S1NALF ♦P*COSALF 
PST  ABGsPSI  *t»*GPP 


• Vt  ' 


SUM24=SUM21»PSTABG 

0ACMaSUM24*Ut)AC 

I END  UF  ROLL, HART  J PROCESSING 

* YAW  FCS  PROCESSING 


ORMAN 

5f (KPTASUP.LT .-1  . 1 25)  L)HMANsRPTASOP  + 1.125 
IFtRPTASOP.GT. 1.125)  UKMANaHP I ASOP-1 . 125 

220Sdkm3S(kdkmun.kokmpak*absU)Rman))*oRman 

1MDRMS  .61.  22.5)  OKMS*22.5 

IFtORMS  .LT.  -22.5)  0RMS*-22.5 
OR  T MS 

ORTGDH*ORT*GDRTI 

IF (IC.OK.HULO)  GO  TO  221  mJTMSFN) 

ORTMSF*FILTl(l)RTGOK..0a,O,-1.0,DRTMSFN) 

221  CONTINUE  _ , . 

IF (DRTMSF.U ,-2.0)DRTMSFs-2.0 
IF (DKTMSF .GT.2.0)DR1MSF*2.0 

ORTMS*ORTMSF  tORMS 
SUM1 

PROD  1 a GNYOKM*DRTMS 
1 F ( TURNCORO ) GO  TO  225 
SUM 1 s-PROO 1 
GOTO230 

NYA«:NY+RD0l*l.Bb4R 

CONTINUE 
PKQD3 

PRU03=GRAY*PR002 
SUM2 

SUM2*P«00iYRST  AH 

GDKt 

GDRC*KORC/CO«AR+80) 

IF (GORC.GT , 1 5 • 0 ) GDRCs 1 5 • 0 
IF (GOkC.LT .1 ,2)G0RC«1.2 

234  PKU05=SUM2*GDRC 

IFCMACM  .LE.  1.5)  GOTO  235 
GORE (EARLY ) 

GORE  = -t>00.0/GbAk  . 

IF (GORE.GT .-1 . 0 ) GDRE*- 1 .0 
IF  (GOKE.LT  .-b.0)GDRE*-6.0 
PR 00  7 ■G0KE*SUM2  3 
0RCPF*YBi*FIL1*PR0D7 
IF (MACH  ,LE.  5.0)  GOTO  236 
PROD7=0.0 
ORCPFsO.O 
GOTO  236 

235  DKCPFsY«dFILT*PR0D5 

236  CONTINUE 
OKTRR 

orthk*ukcpf*gtrr 

IF(IC.OR.HOLD)  GO  TO  244 


225 


230 

232 
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DR  f H I M*P ILT 1 ( DR  TRH  » ,04f0»-i.0»DRTHIMN) 

244  CONTINUE 

IF (OKTRIm.GT ,9.0)DRTRIM«9. 0 
IF (0KTHIM.LT .-9.0) OKTNIMs-9.0 
GRXFO*lSO/UbAR 
IF (GKXFD.GT . 1 .0)  GRXFOal.O 
IF (GRXFO.LT. .1)  GRXFOs.l 
0ACMC=0 

IF ( (MACM.Lt . 1 .5) . ANO. ( TUHNCOKO) ) OACMC«GHXFO*OACM 
SUM3«0HCPF*DRTHIM-0ACMC 

* DRCMD 

IF(SUM3.GT.22.8)  8UMis22.b 
IF (SUM3.LT .-22.8)  SUM3*-22.8 
DRCMOSSUM3 
ORJfc  T*SUM23 

IF(MACH.GT. 1 .5)  GO  TO  262 

* OH  JE  T 

DRJE1 -0 

IF(MACH.GI.l.O)  0RJETa(.RCSA*SUM2 

* NR 


YA 

YB 

262  ABSDRJE I *AbS(URJET  ) 

IF (UBmH.LE.20)  GU  TO  263 
NRsNK  1 

IF ( AbSDHJET ,GT  . .5)  NH  = NR2 
IF(ABbOKJtT.GT  .1,0)  NHSNH3 
IF (AbSDRJE r .GT . 1 .5)  NR=NR4 
GO  TU  264 

263  NHsNR 1 

IF ( AHSDHJt T .GT . .5)  NHsNH2 

264  YADBs.i 
YbOBs.2 

IF (MACH.G  T . 1 ,5)  GO  TO  265 

YAObs.5 

YBOBs.b 

265  CALL  MYSTtR(YAOb/Y80B/UHJETl»ORJET»OUTOHJ) 
DRJET 1 sORJE  T 

UZCMUsNR*OUTOKJ 

ENO  OF  Y Aw  PROCESSING 
HULL  FCS  PHUttSSlNGrPART  2 
PE 

COT ALFsCOSALF/SINALF 
IF(CUTALF.GT,11.43)  CO  T ALF  a 1 1 . 4 3 
IF (CUT  ALP .LT .-1 1 ,43)  C OT ALF s- 1 1 . 4 3 
PE  = ORPRM*COT ALF-Bt  T AU*SINALF-P 
ABSPE® AbS ( PE ) 

YAHXFEEO 

CALL  HYSTER(PEAOB,PEBOB,OLDPE,PE,PEOUT) 
OLUPEsPt 

YAwXPEE0*PkUUT*Ab3PE 
PH004 1 s 0 

IF (UBAR.bT  .2)  PK0041*GP*Y AHXFEED 
OLSP 

IF (MACH.GT . 1 .5)  GO  TO  270 
GOA»200/ (uBAR*80) 
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IF  (GOA.L  I . . 1)  GuAs.l 
If-  (GDA.GT.(1.2/GUAC))  GDAa 1 , 2/GOAC 
DCSP*UDA*DACM 
GO  TO  275 
270  CONTINUE 

GO As  1 50/GB AH 
IT-  IGDA.LT  . , 1 ) GDAa.l 
IMGDA.GT.l)  GO  A3 1 
DCSPaGDA*PRUD41 
275  CONTINUt 

PR0D42a0CSP*GTRA 

OXCMD 

NP=2 

IF (ABSPt.GI . 1 ) NPa« 

UXCMOaO 

IF  ( (UBAR.LE.  10)  . ANO.  (MACH . G T . 1 . 5 ) ) UXCMD**P*PtOUT 
PR0D43 

SUM25*SIGN(  1 .0,UXCMU)-SIGN( 1 ,0,UZCMD) 

SUM25SaSIGN( 1 .0.SUM25) 

AtJS25sABS(SUM25S) 

SUM2b=SUM25-0L025 

0L0253SUM25S 

P«ODajsSIGN( 1 .0, SUM26) *ABS2S 
GTX,PR0D44 
GTXal .25 

IF (MACH. GT ,10)  GTXa50/ (UBARflO) 

I F ( G T X , L T . 1 . 25 ) GTXal. 25 
IF(GTX.GT.S)  GTxaS 
PR0044*PR0D43*GTX 
OATR 
OAT  RaO 

IF (MACH.Lt. 1 .5)  0ATR«PR0D42 

IF ( (MACH.GT .5) . AND. (UBAR.GT .2) ) DA (R»PKOD44 
OATP 

DATPaOA IMFAN*.5 

IF (MACH.GT , 1 .5)  DATPs-.5*0ATMPAN 
DATSUMI 

DA  TSUMaOATR+OA  TP 

IF (IC. UK. HOLD)  GO  TO  300 

DATSUMIaFRT  1 ( OA  T SUM,  . 04 , 0 , - 1 . 0 , 0 AT  SUMN  ) 

300  CONTINUt 
DACMD 

DACMU*OCSP+DATSUMl 
IF (OALMO.GT. 10)  OACMO=10 
IF (DACMO.LT .-10)  OACMOa-lO 


tNO  OF  ROLL  PROCfeSSING 


ACTUATOR  PROCESSING  FOLLOWS 


DtCDRR 

l)t  LOKks  ( ( Dt CMD-DACMl) ) - Ot R ) *20 , 0 
IF (DtCOKR.Gt. 20.0) DtCDRR=20.0 
IF (DtCORR.Lt.-20.0)DtC0RR*-20.0 

Dt  W 

430  DtR*.04*06CDRR*D£R 

IF (DtR.Gt  ,20.0)OERa20.0 
IF (DtR.LT ,-35.0)DERs-35.0 


* DECDRL 

433  DiCDHL=(  IDECMD  + DACMD)-DElWO.O 
IF (UtCURL.GE .20,0) DECDRL *£0,0 
IF(DECDRL.LE.-20.0)OECORL*-20.0 

* DEL 

440  DEL*.04*DECDRL*0EL 

IF  CDEL.GT  ,2u,0)DEL*20.0 
IF ( DEL. LT.- 55.0) DEL *-35.0 

* DE  AND  l)A 

DE*(DELtDER)/2.0 

UA*COEL-UfcH)/2.0 

* DSbC 

450  DSBPCsSBHP-uSBCl 

I F ( DSBPC .GT.6.1)  DSBPC*6.1 
IF (DSbPC.LT.-lO.Bb)  DSBPC*-10.86 
DSBC 1 *DSbC 1 tDSBPC 
USBCs.5*USBCl 

* DRCDRR 

[)PCDRR*( (DWEMD-DSbC) -DUMP) *10.0 
IF ( DRCDRR. GE. 1 0. 0) DRCDRR* 1 0.0 
IF(DHCURR,LT,-10.0)ORCOHH*-ly.O 

* URRP 

460  DRRP=.04*OHCDRRtDRRP 

IF (DRRP.GT.54.88)  0RRP*54.88 
IF  CUWKP.LT .-54.88)  DRRP=-54.88 

* ORCDRL 

463  ORCDRL* ( (DRCMD+DSBC)-DRLP)*10.0 
IF (URCORL.GE. 10.0) ORCDRL* 10.0 
IF(DHCDWL.L1 .-10.0)DHCDRL*-10.0 

* DRLP 

470  DRLP=.04*DKCDRL+DKLP 

IF (OPLP.GT .54,08)  DHLPS54.88 
IF (UKLP.Ll .-54.88)  DRLP=-54,80 

* UP 

QRs(URRP*URLP)*0.5 

* USB 

DSBsURLP-ORRP 

* DBF  PC 

IF  (DBFRC.Ui.3,0)OBFRC*3.0 
IF (DBFPC.LE.-l ,0)D8FRC=-1.0 

* DBF 

440  08F*UBF*0.04*DBFMAN 

IF (DBF .GE. 22. 5) DBF *22. 5 
IF (DBF .LT. -11 . 7)DBF*-1  1 .7 
Hi TURN 

« 

* I ME  FOLLURiNG  SECTION  INITIALIZES  ALL  INTEGRATOR  OUTPUTS  A NO 

* INTERMEDIATE  TRANSFER  FUNCTION  nOOES  TO  ZERO. 

* 

1000  DETR*ETR1M*ELFBK*0 

PR00P5*DR rMSF=ORlRIMsO 
OAT  SUM  IsQ  . 0 

* THE  THAnSFEH  FUNCTION  NODES  ARE  INITIALIZED  TU  ZERO  HERE. 

* 

Dt  f RNSO.O 

ETKlMNstLFBRNsO 

PHUDZ5NXUR1 mSFNsDRTKIMNso, 0 

DA  T Sumn*0 

GO  TO  5 
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... 


FUNCTION  F ILT 1 (X1N,GX1  ,GX2,  Gx 3, ANODE) 

F !LT l«XNOUtMXlN*GX  I 

XNUDE»XIN*GX2-F  1L  f l*GX3 

RE lUKN 

tNO 


FUNC  1 ION  FILT2(X1N,GX1,GX2,GX3,&X9,&X5,XNUDE1»XN0DE2) 

FILT<JsXIN*GX1  + XNOOE1 

XNODElaXlN*GX2-F 1LT2*GX9+XNUDE2 

XNODE2*XlN*GX3-F 1L [2*GXb 

RE TURN 

END 


H Y S 1 E K 


* 

* 


SUBROU  T I NE  HYSTER  <L1M1,L1M2,ARGT0,ARGT1, FUNCT) 
REAL  LIM1,LIM£ 


[HIS  SUBROUTINE  PERFORMS  [HE  HYS  [ EKESIS  FUNCTION  USING 
parameters  f rum  the  calling  program.  The  ARGUMENTS  ARE  AS  FOLLOWS: 


LIM1  AND  L I M2  ARE  THE  POSITIVE  BREAK  POINT  VALUES. 

(LIM1  LESS  THAN  LIM2).  IT  IS  ASSUMED  THAT  THE  FUNCIION 
IS  SYMETRICAL  ABOUT  THE  ORIGIN. 

AHGTO  AND  ARGT1  ARE  THE  INPUT  VALUES  FOR  THE  TwO  MOST 
RECENT  TIME  FRAMES  hITH  ARGT1  THE  MOST  RECENT.  THE  CALLING 
ROUTINE  MUST  CORRECTLY  UPDATE  THESE  VALUES  BEFORE  CALLING 
HYSTER. 

FUNCT  IS  THE  OUTPUT  VALUE  RE  TURNED  10  THE  CALLING  PROGRAM. 
FUNCT  IS  SET  TO  -1,  0,  OR  tl 


TEMPI  s ABS(ARGTl) 

IF  (TEMPI  - LIM1)  10,20,20 
10  FUNCT  a 0.0 
900  RETURN 

20  IF  ( TEMPI  - LIM2  ) 90,30,30 
30  IF  ( AkG  T 1 ) 60 , 5u , SO 
SO  FUNCT  a 1.0 
GO  TO  900 
60  FUNCT  a -1,0 
GO  TO  900 

90  IF  ( ARGT1)  80,70,70 
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* COMt  HtPt  If-  1 NPU I IS  BETrtttN  »L1W|  AND  -LIM2, 

* CHECK  PREVIOUS  INPUT  VALUE  A NU  St  I F UNC T ACCORDINGLY  OH 

* LtAVt  IT  UNCHANGED. 

80  IF  l APG10  ♦ LIM1  ) 900,900,10 

* COMfc  Mfcrft  IF  INPUT  IS  HE  UE£N  LIM1  AND  HM2. 

* CHELR  PPtVlUUS  INPUT  V A L Ufc 

70  IF  l AHGIO  - LIMI  ) 10,900,900 

tNU 


C 'FUNCTION'  IS  USED  TU  OtTtRMINE  THE  OUTPUT  UF  A FUNCTION  SCHEDULE  THAT 
C HAS  MINIMUM,  MAXIMUM,  AND  INTERMEDIATE  UNfcAH  VALUtS.  T hfc  FUNCTION 
C SCHEDULE  CAN  have  A MAXIMUM  UF  1 *0  BREAKPOINTS.  TaBLE(3)  IS  T Ht 

C SMALLEST  AKUUMENT  BREAKPOINT  VALUE  ANO  T AbLt ( 9)  IS  THE  LARGEST, 

C TABLt(l)  IS  THE  UU I PUT  VALUE  FUh  TmE  COHHfc SPONOInG  T ABLt ( 3 ) (SMALL) 

C ARGUMENT  A.vu  I ABLt  ( 2 ) IS  THE  COUPE  SPONU1NG  VALUE  FOP  THE  TARLEC9) 

C (LAKGt)  ARGUMENT, 

FUNCTION  FunC T ION( AKG, TABLE  ) 

D I ME  NS  1 UN  TABLt(9) 

IF ( AHG.GT , 1 AHLt ( 3) JGUTU3000 
FUNCT  1 U N * T ABLE ( J ) 

GOTO  3002 

3000  IF ( ARG.LT . I ABLt (9) )GOTO  3001 
F UNC T IONS f ABLt (2) 

GOTO  3002 

3001  FUNCT I UN* I AbLt ( 1 ) ♦ ( APG- T ABLE ( 3 ) ) * ( ( T AbLE ( 2 ) - T ABL t ( 1 ) ) / ( T AHLt ( 9 ) - 


*TAbLt (3) ) ) 
3002  RETURN 


ENU 
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